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CHAPTER I 


INTRODUCTION 


1.1 Object and Scope 
This thesis is a study of the effectiveness of the surrounding 


soil in restraining the lateral movement of shallow-buried arches subjected 
to blast overpressures from a nuclear explosion. An understanding of this 
effect is important in the safe and economical proportioning of such 
structures. 

The theory is developed for parabolic arches with circular or 
barrel arches included as a special case. It can be applied equally well 
to arches oriented with their longitudinal axes perpendicular to or parallel 
to the direction of travel of the blast wave. 

The arch is considered to be buried near the surface of a semi- 
infinite soil mass. The term ‘shallow-buried" as used herein refers to a 
configuration and amount of earth cover such that reflected pressures and 
dynamic pressures need not be considered in the loading of the arch and yet 
not so great as to provide significant free-field attenuation of air-blast- 
induced overpressures in the soil. This is a convenient and meaningful area 
of investigation. It is convenient because it eliminates the dynamic and 
reflected pressure loading components at one extreme and the free-field 
attenuation of overpressure in soil at the other. Consideration of either 
of these extremes would serve further to complicate the problem while con- 
tributing little to the significance of the results. It is meaningful in 
that an economical balance may be reached in practical installations by 
locating structures near the ground surface in order to minimize excavation 
requirements and facilitate access while at the same time using earth cover 


to eliminate dynamic and reflected pressure loading components from the blast 





loading which the structural system must withstand. In this study the depth 
of cover over the crown is limited to one-half of the arch span. This 
restricts the location of the structure to a region adjacent to the earth's 
surface in which there appears to be no measurable attenuation of free-field 
overpressure in the soil and in which the restraint offered to soil movement 
by the overpressure acting on the ground surface is significant. 

The two-hinged rather than the fixed arch is selected as repre- 
senting more closely the shallow-buried arch structures used in full-scale 
nuclear tests to date, (19295) While these structures were not actually pinned 
at the arch abutments, their behavior under blast loading appears to the 
writer to resemble more closely the pinned than the fixed case. Further, such 
techniques as separating the arch footings from the floor slab in order to 
reduce accelerations imposed on the floor slab, and keeping the arch footings 
small to encourage downward movement under blast loading in order to reduce 
direct compression stresses in the arch, tend to produce a more nearly pinned 
than fixed condition at the abutments. The greater structural flexibility 
of the two-hinged arch is an additional favorable feature in enabling under- 
ground structures to adjust toimposed loading systems. 

The principal variables studied in the numerical examples include: 

1. Depth of cover 
2. Overpressure loading 
3. Arch material 


4. Arch span 


1.2 Method of Approach 
The arch is replaced by a model consisting of four inextensible 


y 
bars connected at moment-resistant roxanne This model is developed 


(1) Numbers in parenthesis refer to the List of References. 
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herein in Section 3.1. The arch is thus reduced to a two-degree-of-freedom 
system. This greatly facilitates the dynamic analysis while at the same 
time permitting meaningful assumptions to be made regarding the loading and 
resistance parameters at the two haunches and at the crown. 

The two independent variables chosen to describe the system are 
the horizontal displacement components at the haunches. While any two in- 
dependent displacements could have been chosen, these are most convenient 
since they provide a direct picture of the primary mode of response of the 
model and are used directly in the determination of the state of stress within 
the soil mass at each side of the arch. 

The two independent modes of response of the model are referred to 
herein as the symmetric mode, in which the haunches move simultaneously in 
opposite directions, i.e. both inward or both outward, and the asymmetric 
mode, in which the haunches move simultaneously in the same direction, i.e. 
both toward the right or both toward the left. 

Solutions for the response of this system are obtained using the 


(4,5,6,7) 


Newmark B-methods. These are step-by-step numerical integration 
techniques in which the accelerations, velocities and displacements of the 
dynamic system are determined at each time interval from the time of arrival 
of the shock wave until the significant response of the system has occurred. 
These techniques are readily adaptable for use in digital computers. One of 
the University of Dllinois Structural Research Laboratory Codes was modified 
and used for the solution of the numerical examples presented in this thesis. 
The blast overpressure loading is concentrated as discussed in 

Section 2.4 and Chapter VI and is placed directly on the arch. Modification 
of the blast loading due to its transmission through the soil mass from the 


ground surface to the arch rib is handled separately as "soil resistance’. 


This is discussed further in the following sections. 





In general, the resistance of the arch-soil ‘structure’ is handled 
by: 

1. Considering the bending resistance of the arch as it deforms, 

2. Considering the frictional and shearing stresses developed 
within the soil as the arch moves into or away from it, 

3. Combining these arch and soil resistances vectorially as they 
act upon the model, and 

4. Considering the mass effects of the arch and soil. 

While these resistances are closely interrelated and act in combination during 
the response of the system, it is convenient and admissible to handle them 
separately in developing the equations of motion for the system. 

The soil resistance is developed by considering the states of stress 
which exist within the soil surrounding the arch as the arch-soil "structure" 
deforms laterally. 

Direct compression in the arch is handled separately. This is an 
approximation but a necessary one with this model, since to include it direct- 
ly would violate one of the assumptions upon which the model is based, i.e. 
the bars are inextensible. Note that it would also change the model from 
a two-degree-of-freedom system to a six-degree-of-freedom system. However, 
the response studies presented remain valid with respect to lateral movements 
as long as the arch section is sufficient to withstand direct compression 
and buckling. A separate "rib-shortening correction” (cf. Section 3.5) is 
applied to certain of the results to give an approximate picture of its 
Significance. This is discussed further in Chapter IX. 

The arch supports are considered to be unyielding. While this is 
a common assumption in structural problems, it is an approximation. Again, 


the assumption is necessary in order to limit the system to two degrees of 





freedom; however, this does offer an interesting possibility for further 
study as stated in Chapter IX. 

Throughout this thesis displacements, velocities and accelerations 
to the right and upward are considered positive, and conversely. JLoads are 
considered positive when they act on the structure to the right or upward. 
Resistances are considered as negative loads. 

The method of approach described above is approximate. The neces- 
sity of approximating actual conditions is discussed in Chapter II and is 
also pointed out as occasions arise in the development of the loading and 
resistance expressions for the system. The approach used gives calculated 
magnitudes of response which differ from values measured in field tests. 
However, as will be seen, the results do describe the behavior of the actual 
Structure and permit meaningful comparisons to be made of arch behavior 
under varying combinations of structure geometry and loading conditions 


(cf. Chapter IX). 


1.3 Notation 
Each symbol used in the text is fully explained when it is first 
introduced. <A summary of the most important notation is presented here for 


the convenience of the reader. 


a = vertical projection of a lower bar of the dynamic model 
(cf. Figure 2.1) 
A = cross-sectional area of the arch rib 


Ay, Aas A, = influence coefficients for total joint rotation (cf. Equa- 
tiong(s.2,i0)) 


b = horizontal projection of a lower bar of the dynamic model 
(cf. Figure 2.1) 

B = span of the arch (cf. Figure 2.1) 

: = vertical projection of an upper bar of the dynamic model 


(cf. Figure 2.1) 





J. 


J+) 


J 


J. 


eo) 36 


seismic velocity of the soil 


apparent vertical component of the velocity of the shock 
wave in the soil 


horizontal projection of an upper bar of the dynamic model 
(cf. Figure 2.1) 


depth of soil cover over the crown of the arch 


vertical dimensions of the soil prisms used in the mass 
calculations (cf. Figures 5.1 and 5.2) 


total depth of soil cover from the ground surface to the 
arch abutments 


modulus of elasticity of the arch material 
equivalent mass factor (cf. Eyuations (7.1.4)) 
acceleration of gravity, 52.2 ft/secé 
equivalent mass factor (cf. Equations (7.1.4)) 
thickness of the arch rib 

displacement 

velocity 

acceleration 

effective external load 

equivalent external load 

effective moment of inertia of the arch rib cross section 


as a subscript, unless otherwise specified, refers to the 
horizontal component at joint j 


as a subscript, unless otherwise specified, refers to the 
vertical component at joint j 


influence coefficients for the equivalent load on the dynamic 
model (cf. Equations (6.3.2)) 


influence coefficients for equivalent soil stiffness 
(cf. Equation (4.4.10)) 


observed value of the ratio between the horizontal and 
vertical soil pressures for the condition of earth pressure 
at rest 


hydrostatic pressure ratio for the “active” state of stress 
in the soil 





OA 


Kop 


bending stiffness of the arch rib as concentrated ata 
joint of the dynamic model 


ratio of the soil pressure in the direction of P' to the 
vertical soil pressure for the condition of earth pressure 
at rest 

ratio of the soil pressure in the direction of P! to the 
vertical soil pressure for the condition of e pressure 
at rest 


hydrostatic pressure ratio for the "passive’ state of stress 
in the soil 


ratio of the horizontal to the vertical overpressure induced 
in the soil by the shock wave in air 


soil stiffness during the development of an "active’ state 
of stress in the soil adjacent to haunch j 


soil stiffness during the development of a "passive' state 
of stress in the soil adjacent to haunch j 


length of bar j of the dynamic model (cf. Figure 2.1) 
concentrated mass at a joint of the dynamic model 
arch contribution to Ba 

soil contribution to we 

equivalent mass 

moment resistance of the arch-rib cross section 


maximum moment-resisting capacity of the arch rib 
cross section 


as a subscript, refers to time t = is 


as a subscript, refers to time t = tol 


overpressure in the shock wave in air 


peak value of the overpressure in the shock wave in air 


vertical component of the average overpressure acting on 
the arch 


uniform radial component of the average overpressure acting 
on the arch 


average overpressure acting at the ground surface on 
loading region j 





| 
POA 


tt 
POA 


t 
Pop 
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ar 
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average overpressure acting at the ground surface on the 
potential "active" failure wedge in the soil adjacent to 
haunch j 


average overpressure acting at the ground surface on the 
potential "passive" failure wedge in the soil adjacent to 
haunch j 


force on the side of the arch due to the weight of the soil 


force on the side of the arch due to the "surcharge" of the 
blast overpressure acting on the ground surface 


force acting on the side of the arch due to the weight of 
the soil and corresponding to a fully developed "active" 
state of stress in the soil mass 


force acting on the side of the arch due to the blast over- 
pressure and corresponding to a fully developed "active" 
state of stress in the soil mass 


component in the direction of P! of the force acting on the 
side of the arch due to the weight of the soil and correspond- 
ing to the condition of earth pressure at rest 


component in the direction of P! of the force acting on the 
side of the arch due to the blast overpressure and corre- 
sponding to the condition of earth pressure at rest 


component in the direction of P‘ of the force acting on 
the side of the arch due to the weight of the soil and 
corresponding to the condition of earth pressure at rest 


component in the direction of P\ of the force acting on 
the side of the arch due to the blast overpressure and 
corresponding to the condition of earth pressure at rest 


force acting on the side of the arch due to the weight of 
the soil and corresponding to a fully developed "passive" 
state of stress in the soil mass 


force acting on the side of the arch due to the blast over- 
pressure and corresponding to a fully developed "passive" 
state of stress in the soil mass 


influence coefficients for equivalent arch resistance 
(cf. Equation (3.4.7)) 


equivalent resistance due to deformation of the arch rib 


effective resistance at an "interior" haunch due to deforma- 
tion in the soil mass 


effective resistance at an "exterior" haunch due to deforma- 
tion in the soil mass 





equivalent resistance due to deformation in the soil mass 
radius of a circular arch 

maximum value of soil resistance (based upon a consideration 
of gravity forces) for a fully developed "active" state of 
stressin the soil 

maximum value of soil resistance (based upon a consideration 
of overpressure forces) for a fully developed “active” state 
of stress in the soil 

maximum value of soil resistance (based upon a consideration 
of gravity forces) for a fully developed "passive" state of 
stress inthe soil. 

maximum value of soil resistance (based upon a consideration 
of overpressure forces) for a fully developed "passive" state 
of stress in the soil 

shear in bar j 


time 


elapsed time following passage of the shock front at a point 
of interest 


duration of the positive phase of the shock wave 


time of arrival of the shock front at the windward boundary 
of loading region j 


time of arrival of the shock front at the leeward boundary 
of loading region j 


thrust in the arch rib 

natural period of vibration of the system 
velocity of the shock wave in air 

effective inertial force 

equivalent inertial force 

element of the constraint loadings 

one-half the angle of opening of a circular arch 
weight per unit volume of the arch material 


weight per unit volume of the soil 





Sh 

SL 
J 

5¢@ 
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At 


10 
friction angle corresponding to the friction developed 
between the surface of the arch and the soil 
"rib-shortening' correction 
change in length of bar j due to thrust in the arch rib 
rotation of bar j 
time interval used for numerical integration 


inclination with respect to the horizontal of the shock 
front in soil 


inclination of bar j with respect to the horizontal 
multiplying factor for the constraint loadings 

influence coefficient for e,uivalent mass (cf. E uation (5.2.6)) 
influence coefficient for equivalent mass (cf. Equation (%2.6)) 
rotation of the lower bar of the model at haunch 2 

rotation of the lower bar of the model corresponding to full 
development of an ‘active’ state of stress in the soil 
adjacent to haunch 2 

rotation of the lower bar of the model corresponding to 

full development of a "passive state of stress in the 

soil adjacent to haunch 1! 


influence coefficient for e,uivalent arch resistance 
(cf. Equation (7.2.3)) 


angle of internal friction of the soil 
rotation of joint j 
joint rotation required to produce yielding of the joint 


inclination of the lower bar of the model with respect to 
the vertical 


circular natural freyuency of the system 
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CHAPTER ITI 


GENERAL DISCUSSION OF THE VARIABLES 


2.1 Arch Resistance 

The inertial resistance of the arch is discussed separately in 
Section 2.3. The discussion here considers the arch resistance in terms of 
the forces required to pull a msssless, deformable structure into the deflec- 
tion configuration defined by the two independent displacements. 

For the purpose of the dynamic analysis the arch, which hss an 
infinite number of degrees of freedom as it responds to dynamic loading, is 
replaced by a mathematical model. This model consists of four rigid, inex- 
tensible bars connected together at three moment-resistant joints and 
connected to the arch abutments at two frictionless hinges as shown in 
Figure 2.1. This model has two degrees of freedom. Note that it permits 
assumptions to be made regarding conditions at the two haunches and at the 


(8) 


crown. It has been demonstrated that these are the critical sections of 
arch behavior, and thus the model provides a relatively simple though useful 
tool for the study of the dynamic response of the arch itself. 

In this thesis, an arch of unit width is considered. Also it is 
assumed that the arch section under consideration is sufficiently far removed 
from endwalls or intermediate stiffeners that they have no effect on its 
dynamic response. This, in effect, reduces the problem of one of two 
dimensions. 

The flexural resistance of the arch is determined herein by first 
finding the moment-resisting capacity of the arch cross section and then 
concentrating this moment resistance at the three moment-resistant joints 


of the model. The joint resistance is then represented by an idealized 


moment-rotation relationship. This is believed to be a reasonable 
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representation of the flexural resistance of the arch, since it will be shown 
(cf. Chapter IX) that the flexural resistance of the arch contributes rela- 
tively little of the total resistance of the soil-arch "structure’ to lateral 
deformation under blast loading. Thus, modification of the flexural resistance 
to allow for dead and live loads acting on the arch prior to the arrival of 
the shock wave or for thrust in the arch due to the blast loading does not 
appear warranted and no modification is made herein. 

Since the model used to represent the arch is composed of inex- 
tensible bars, it is necessary to handle the direct compression resistance 
of the arch separately. This is an approximation; however, its effect upon 
the lateral deformation response of the soil-arch “structure” is not considered 
serious as long as the arch section is sufficient to resist the maximum thrust 
imposed on the arch by the blast overpressure without general yielding or 
buckling inward. Well compacted soil surrounding the arch should restrain 
any tendency of the arch to buckle outward into the soil mass. It is important 
to note at this point also that the arch section must be sufficient to with- 
stand imposed dead and live loadings during construction and by normal service 
conditions prior to the nuclear explosion. Also worthy of note is the probable 
loading condition on the arch after the nuclear blast wave has passed and 
dynamic response has ceased. The importance of attentionto these conditions 
in the design of any buried structure is obvious. 

An approximate "rib-shortening correction is applied to certain 
of the numerical results obtained in this study. While this is not directly 
a part of a study of lateral deformation, it does give some indication of the 
total deformation which may be expected under the given conditions. It also 
serves to indicate the effect of direct compression in the arch rib on arch 


response. 
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2.2 Soil Resistance 

The present state of knowledge regarding the interaction of soil 
and underground structures contains many uncertainties. This is true for 
"static" behavior and is even more true for "dynamic" behavior such as that 
now under consideration. The development in this thesis is considered 
reasonable but is necessarily approximate, conse ,uently a high degree of 
refinement is not attempted. 

As with the arch discussion in the preceding section, the inertial 
effects of the soil are covered separately in Section 2.3. The discussion 
here considers the shear and frictional forces developed within the soil 
Surrounding the arch as the soil-arch “structure’ deforms under blast loading. 
It is assumed throughout that prior to the arrival of the blast wave the soil 
surrounding the arch is in the "at rest’ state described below. 

The concept of soil resistance used in this thesis derives from a 
consideration of the states of stress which are developed within the soil 
mass surrounding the arch as the arch undergoes lateral deformation. The 
confining action of the blast overpressure on the soil is also considered. 

The states of equilibrium within the soil mass which have been 


(9) 


used in developing this concept are: 

1. Active State of Plastic Equilibrium - the state of stress 
within a soil mass which has undergone a "stretching, i.e. the soil mass 
has been permitted to expand laterally, until the ratio of the horizontal 
to the vertical pressure within the soil mass has reached a minimum value, 
called the "coefficient of active earth pressure.” 


©. Passive State of Plastic Equilibrium - the state of stress 


within a soil mass which has undergone a lateral compression until the ratio 
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of horizontal to vertical pressure has reached a maximum value, called the 
"coefficient of passive earth pressure." 

3. Earth Pressure at Rest - this is a state of stress which is 
intermediate between the two states of plastic equilibrium described above 
for the soil mass and exists after the soil has been deposited by natural 
and artificial processes in the vicinity of the arch. It is a state of 
elastic equilibrium. Since this definition includes "artificial" processes 
for deposit of soil, it will be restricted herein to rejuire careful com- 
paction for optimum density of all fill materials used during construction 
and backfilling of the arch. 

The hydrostatic pressure ratios, or ratios of horizontal to 
vertical soil pressure, referred to in the foregoing definitions differ 
because of the differences in magnitude and direction of the shear and 
frictional forces in the soil mass. In general, these stresses are developed 
within the surrounding soil as the arch-soil "“structure’ deforms under blast 
loading. The extent of development of these stresses at any particular time 
depends upon: 

1. The state of stress within the soil mass, and 

2. The restraint imposed upon the deformation of the soil mass 


by the confining action of the blast overpressure. 


Since the blast loading is applied to the arch rib through the soil, these 
stresses will assist the arch in "resisting the loading. 

The degree of restraint imposed upon the deformation of the soil 
mass by the blast overpressure will, of course, vary with time as the blast 
overpressure varies. 

The primary variables determining the state of stress within the 


soil mass are the direction and amount of deformation to which the soil mass 
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has been subjected. These will, of course, vary with time, but, at any 
given time, will also vary from one region to another around the arch. In 
this thesis, three regions are considered, viz. the regions adjacent to each 
of the haunches and the region adjacent to the crown. 

As noted above, prior to any deformation it is assumed th’ t the 
"at rest’ stress state exists in che soil surrounding the arch. As a haunch 
moves inward, an ‘active state of plastic equilibrium" tends to develop. 
This does not develop fully immediately but requires a certain amount of 
deformation within the soil mass as discussed below. As a haunch moves out- 
ward into the soil mass, a ‘passive state of plastic e .uilibrium’ tends to 
develop. This, too, does not develop fully until sufficient deformation 
has occurred. 

The amount of deformation re ,uired for the full development of an 
active state of plastic equilibrium within the soil mass was studied by 


Geer) Terzaghi. (2° 


In a series of tests on a large-scale retaining wall 
Dr. Terzaghi demonstrated that the lateral force on a wall varies as the 
wall undergoes lateral movement. The relationship between the force and 


(11) The resistance- 


the movement is shown j,ualitatively in Figure 2.e. 
displacement functions developed in this thesis are essentially idealizations 
Siethis figure. 

While Figure 2.c is qualitative only, Dr. Terzaghi and others (19/11,22) 
have estimated the amount of wall movement re uired to reduce the earth 
pressure to the active state for retaining walls and specific soil conditions. 
These recommendations form the basis for the quantitative assumptions of 
deformation required to develop fully the active state of stress which are 
used in this thesis. There is little guidance in the literature regarding 


the amount of deformation required to develop fully the pussive st:te of 
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stress within the soil mass. Limited investigation 


has indicated, 
however, that it is greater than that required to develop the active 
state of stress. Accordingly, in this study it is assumed to be such that 
the idealized resistance-displacement function for the soil has a constant 
slope in the active and passive ‘elastic’ regions for a given value of over- 
pressure restraint (cf. Section 4.4). 

In determining numerical values of the hydrostatic pressure ratios 
corresponding to the limiting states of stress described above the classical 


(14) It is assumed that: 


Coulomb assumptions have been used. 

1. The soil slides along a solid pressure surface, i.e. the soil 
failure wedge is bounded on one side by the pressure surface and on the 
other side by the surface of rupture. 

e. The obliquity of the resultant pressure on the pressure surface 
has a value less than or equal to the friction angle of the soil. 

4. Sufficient movement has occurred within the soil mass to 
develop the state of stress being considered, i.e. active or passive. These 
assumptions yield directly the pressure acting on the pressure surface. A 
comparison between results obtained from the Coulomb Equation and those 
obtained by the Trial Wedge Method is given in Appendix A. The equation 
gives results which are considered sufficiently accurate for this study. 

An additional assumption re ,uired to define the soil failure wedge 
is the character of the surface of rupture within she soil mass. A plane 
surface of rupture is assumed throughout. 

The pressure surface is taken as the lower bar of the model on 
each side and is considered to continue upward in the same direction from 
the haunch to the ground surface. This facilitates the development of the 


soil resistance concept used in this study and, for the relatively shallow 


depths of burial considered, does not appear unreasonable. 
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The soil resistance at the crown is limited primarily to the 
inertial resistance of the soil at and adjacent to the crown. While there 
must obviously be shearing forces developed in this soil mass as the arch- 
soil “structure” deforms, the relatively small volume of this region for 
the shallow depths of burial considered, as compared with the volumes of 
the soil wedges influenced by lateral movements of the haunches, indicates 
that these shearing forces will be relatively small. Also, it is possible 
that they will act, at least in part, only to carry the load to another 
portion of the arch. The extension upward to the ground surface of the 
pressure surfaces at the sides of the arch as described in the preceding 
paragraph is an attempt to include that portion of these shearing forces 
which will be effective in reducing the load brought to the arch. No load 
reduction is made, i.e. no ‘deformation resistance’ is developed, for the 
remainder of the soil in this region. While this approach is somewhet 
arbitrary, it is not considered to have a significant influence on the 


response of the system as a whole. 


2.5 Mass Effects 

The inertial resistance of the arch-soil “structure” has an 
important effect on both the magnitude and time history of its dynamic 
response. As noted previously in Section 2.1, the model used in ths study 
for the dynamic analysis permits appreciable latitude in the choice of the 
masses which are considered to move with the system at any given time. This 
property of the model will be used to advantage in adapting the model to fit 
changing conditions during the response history of the actual structure. 

The mass of the arch is concentrated at the joints connecting the 


inextensible, massless bars, except for the mass of the lower portion of the 
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arch on each side. This portion of the mass of the arch is concentrated 
at the abutments, i.e. is not considered to have a significant effect upon 
the response of the system. This technique is commonly used in dynamic 
analyses and gives reasonable results. 

The mass of soil which is assumed to move with the arch is 
selected according to: 

1. The location of the joint, 

2. The direction of movement of the joint, and 

4. The depth of burial of the arch. 

Varying the mass assumption with the direction of movement of the 
joint has the effect of giving the dynamic system a variable mass. While 
this increases the complexity of the analysis somewhat, it may be handled 
readily by the numerical integration procedure, and it appears to give a 
better approximation to the actual field condition than the use of a constant 
mass assumption. logic indicates, for example, that the quantity of soil 
influenced by the inward movement of a haunch is different from the quantity 
influenced by the outward movement of the haunch into the soil. 

Limits are placed on the dimensions of the soil masses considered. 
In no case is a dimension considered which is greater than the corresponding 


(a5) 


span of the arch. 


2.4 Blast Pressure Loading 
Since the object of this study is to investigate the response of 


a shallow-buried arch to blast loading from a nuclear explosion, the blast 
overpressure is the primary weapon effect considered. In applying the 
results of this study to any field situation, other weapon effects obviously 


must be considered as well. Further, the fact that the limiting cases of 
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the numerical examples studied herein would not prove feasible in the field 
because of other-than-blast effects should not detract from their value in 
developing an understanding of structural response to the blast overpressure 
itself. 

The results of this study are considered applicable to structures 
located outside of the crater, i.e. outside both the zone of rupture and 


(16) It is doubtful that structures can be 


the zone of plastic deformation. 
built to withstand weapons effects levels at closer ranges. 

Direct ground shock is not considered. For structures located near 
the ground surface and outside the crater this omission appears reasonablé. 
Much of the energy of the direct ground shock is dissipated in disruptive 
effects in and immediately adjacent to the crater. Further, as the direct 
ground shock wave travels outward from the point of the explosion, its rise 
time increases so that it becomes even less effective in causing damage. 

Thus, the air-induced ground shock is considered the major struc- 
tural damaging effect. It is assumed that if the structure can withstand 
the air-induced ground shock it can resist both the direct ground shock and 
differential soil movements. This latter consideration may become more 
important with larger structures than those considered herein. However, no 
truly "large" underground arch structures have beentested to date, and the 
possible importance of differential soil movements in this case has not 
been established. (19) 

The shock wave in air is assumed to have a single front, i.e. the 
structure is assumed to be in the region of Mach reflection. This assumption 
is not necessary but is convenient when defining the load on the structure. 


The study may be extended readily to cover other types of blast waves if 


desired. 
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No attenuation of free-field overpressures with depth is assumed. 
For the relatively shallow depths considered in this study, this is entirely 
reasonable (cf. Section 1.1), $18) 

The decay of free-field overpressure with time in the soil mass is 
assumed to be the same as that of the shock wave in air. Again, for the 
relatively shallow depths considered in this study, this is a reasonable 
Eeeumption. ‘2°? 

The free-field horizontal-to-vertical overpressure ratio in the 
soil is assumed as 1/4. This is the lowest of the values normally recon- 
(15,17,19) 


mended and was chosen to accentuate any tendency of the arch to 
move laterally. A low value of this ratio acts, of course, to decrease the 
horizontal component of the pressure acting on the arch during passage of the 
blast wave and thus to decrease the restraint which the load itself offers 

to lateral movement. 

A plane shock front is assumed for the shock wave in soil. A study 
of the inclination of this shock front is presented in Appendix B. The rise 
time of the overpressure at this shock front is assumed instantaneous. The 
fact that studies and field measurements do indicate a finite and increasing 
rise time as the shock wave travels downward from the ground cunrece 1° 1s 
partially allowed for in developing the loading on the arch by considering 

(16) 


the average overpressure over a region of interest rather than the 


instantaneous overpressure at any particular point. 
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Figure 2.1. Model Considered in Analysis 
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CHAPTER III 


ARCH RESISTANCE 


In developing an expression for the resistance offered by the arch 
to lateral deformation under nuclear blast loading one must consider the manner 
in which the arch will respond and the relative importance of the various re- 
sistance parameters. In general, as the arch responds to a dynamic loading, 
the arch will mobilize both inertial and deformation resistance. This chapter 
is limited to a discussion of resistance due to the deformation of the arch 
rib. Inertial effects are discussed separately in Chapter V. 

As the arch moves laterally in response to a dynamic, unsymmetrical 
loading, significant deformations occur at the haunches and at the crown. It 


(8) 


has been shown that these are also the critical sections for bending moment 


and that the thrust is nearly uniform throughout the arch. 


3.1 Development of Model 
Since the arch itself is continuous, it has an infinite number of 
degrees of freedom as it responds to a dynamic loading. Thus, it is incon- 
venient to analyze the arch directly. In this study, for the purpose of the 
dynamic analysis, the arch is replaced by the four-bar model shown in 
Figure 2,1, (4) 
The model consists of four inextensible bars connected at three 
interior joints. At the interior joints are concentrated both the moment 
resistance of the arch rib and the mass of the arch. These joints also act as 
node points. The model permits the introduction of load and soil resistances 
into the analysis at each joint location. As mentioned in Section 2.1, the 


four-bar model is connected to the arch abutments by frictionless hinges. 


This is permissible for the two-hinged arches studied in this thesis, 
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If the bars were capable of deformation, the model would be a six- 
degree-of-freedom system. Since the bars are inextensible, one can write four 
"constraint relationships among the joint displacements, one "constraint" 
relationship for each bar. This reduces the number of degrees of freedom to 
two. These ‘constraint’ relationships are based upon the condition that the 
deformations of the model at the joints correspond to no change in the length 
of the members. These relationships are developed for the model in Section 3.2. 

In this study the constraint relations will be derived using four 
loading patterns which produce only axial stresses in the bars. They thus 
cause no movement of the model and, of course, do no work when applied to the 
model. The principle of virtual work and the fact that these force patterns 
Cause no movement are then used in deriving relationships among the displace- 
ments and among any general set of loading components. These relationships 
are then consistent with the two-degree-of-freedom systen. 

A limitation in the use of this model is that direct compression in 
the arch cannot be handled directly. This follows from the inextensible nature 
of the bars. This limitation would be serious if this study were of the total 
stability of the arch, since the arch rib must be able to withstand the thrust 
loading imposed upon it without failing by general yielding or buckling. 
However, if the arch section is sufficient to carry the imposed thrusts, it may 
then be studied from the standpoint of lateral deformation response without 
simultaneous consideration of direct compression load components. 

The prime reason for accepting this limitation in the model is, of 
course, to retain the convenience of the two-degree-of-freedom system. If 
the bars were extensible, the model would have not two but six degrees of 
freedom and the analysis would become correspondingly more involved. The nature 


of the behavioral assumptions upon which this study is based is approximate 
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enough that a more refined analysis is not considered justified at the present 
time. 

The independent displacement components chosen to define the con- 
figuration of the model at any instant are the horizontal displacements of 
the haunches. While other pairs of independent displacements could be used, 
those chosen are particularly convenient in that they: (1) indicate directly 
the primary mode of response of the system (cf. Article 2.1) and (2) simplify 


the expression of the soil resistance (cf. Chapter IV). 


3.2 Constraint Relationships 


The constraint relationships relate the various arch parameters and 
serve to reduce the continuous arch to a two-degree-of-freedom system. These 
relationships are developed herein by considering four independent loading 
patterns, each of which produces only axial forces in the bars of the model. 
Thus, since the bars are inextensible, no deformation of the arch can take 
place under any of these loading patterns. 

The four loading patterns chosen are shown in Figure 3.1. For con- 
venience these loading patterns will be referred to herein as the "constraint 
loadings". Other loading patterns could be developed by combining these; 
however, these possess certain orthogonal properties which are convenient in 
deriving the following relationships. The constraint loadings may be tabulated 


conveniently as shown in Table 5.1. 


5.2.1 Loading Relationships 
The load concentrations which will be applied to the model are 


shown in Figure 3.2. Relationships among these load components are now derived 
using the constraint loadings shown in Figure 5.1 and the fact that these con- 


straint loadings, when applied to the model, produce no movement of the model. 
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Now any general system of loads can be replaced by a set of loading 
patterns consisting of some proportion of each of the four constraint loadings 
plus a loading pattern which does produce movement of the model. This latter 
pattern is that portion of the original system of loads which is effective in 
producing dynamic response of the system. This concept is expressed in 


Eyuation ie. 2.1). 


i 
Hy = 7 Ay Os tH, (3.2.1) 
m=] 
where H = any system of loads on the model (cf. Figure 3.2) 
ra = multiplying factors for the constraint loadings 
An = elements of the constraint loadings 
H, = that portion of the loading system which is effective 
in producing motion of the model 
j = load component index; j = 1,2,... 6 
m = constraint loading index; m = 1,2,3,h. 


In developing the load relationships from Equation (3.2.1) the 
constraint loadings serve to reduce the number of equations from siX to two. 
This is accomplished algebraically by seeking four values of the multiplying 


factor, 0 such that four of the quantities A, are zero. Since the horizontal 


displacements at the haunches have been chosen as the independent variables, 


the corresponding load components H, and H, will be retained here. Thus, in 


the following derivation multiplying factors are sought such that a = Hy = 


tg = (O- 


Using either Figure 3.1 or Table 5.1 and Equation (3.2.1), one may 


write: 
Hy, - a X,, +ary =O 
d 
H - 2(=) A = 0 
e rae (5.2.2) 
= O 
i. ar, 
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Note that the constraint loadings were chosen so that Equations 


ma 


(3.2.2) are orthogonal, es 


J 
Thus, the multiplying factors, ae may be obtained from: 


a oe = 0, for m#n and j = 2,4,5,6. 


2 
hoe He = 0,4,5,6 (5.2.3) 
m 

d 

j 


Equation (3.2.3) gives directly: 


il 
7s a(H ) 
eer (H, + H-) 
2 2a 6 
(3.2.4) 
= oq) 
1 


M, = 7 Ba Hy > Hg) 
Substituting Equations (3.2.4) into Equation (3.2.1) and solving for 
Hy and H, yield the load relationships for the two-degree-of-freedom system. 
They express the "“eyuivalent" load on the system, i.e. that portion of the 
total load which is effective in causing movement of the structure in the 


horizontal direction at the haunches. 


— ie be b 1,d db 

He Pe oa ee ee _ 
eae 

— nk be 1,da b ve) 

He 2 Ee cate * eo oe ee 


These load relationships will be used later in developing the 


equations of motion for the model (cf. Chapter VII). 


9.2.2 Deflection Relationships 
In deriving relationships among the displacement components of the 


model use is again made of the constraint loadings of Figure 5.1 together with 





et 


the principle of virtual work. Since none of the loading patterns produces 
any displacement of the model, the product of the loadings in the pattern 
times the displacements in any consistent set of displacements of the struc- 
ture in the direction of the loadings in the pattern must be zero. This is 


expressed algebraically in Equation (4.2.6). 


& 
» hy 4, = 0 (3.2.6) 
j=l 
where or = displacement components 
oe = elements of the constraint loadings 
j = displacement component index; j = 1,2,...6 
m = constraint loading index; m = 1,2,3,4. 


Again referring to the choice of the horizontal displacement com- 





_ ponents of the haunches as the independent variables, the other displacement 
components may be written in terms of them using Equation (3.2.6). 
Solving Equation (3.2.6) explicitly for each of the deflection 


components yields: 


h, = (1 - be/aa)(h, + h;) 
hy = - (b/a)h, 
ld bd 3.2.7) 
By og Oe 
hy = (b/a)h, 
where a,b,c,d = generalized arch dimensions (cf. Figure 2.1) 

hy = horizontal displacement component of joint 1 

h,, = horizontal displacement component of joint 2 

te = horizontal displacement component of joint 3 


hy, = vertical displacement component of joint 1 
h. = vertical displacement component of joint 2 


he = vertical displacement component of joint 5 
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3.2.5 Joint Rotation Relationships 


It is necessary to know the relationships between the rotations of 
the joints and the independent displacement variables in order to determine 
the moment resistance developed in the joints of the model for any particular 
deflected shape. Since the horizontal displacement components of the haunches 
have already been chosen as the independent variables, they will be used as 
such in this derivation also. 

Referring to Figure 3.4, it can be seen that for small displacements 
the rotations of the bars may be expressed as: 


(h -h.) ese 6 


59. = a | Eee ee) (3.2.8) 
2 L 
J 
where Ee, = rotation of bar j 
“4? oe = horizontal displacement components at the ends of bar j 
L, = length of bar j 
a = inclination of bar j 


Substituting into Equation (3.2.8) and rearranging terms give the 
following expressions for the bar rotations of the model in terms of the in- 
dependent displacement coordinates and the geometry of the arch: 

66, = - hy/s 
Bene = =(1/c + b/ad)h, - =(1/c - b/ad )h., 


2 
(3.2.9) 
ee 1 
56, a Aili Z b/ad )h, + Awe + b/ad )h, 
6), = - h,/a 


The total rotation at each of the joints is simply the algebraic 
difference of the rotations of the bars connecting at the joint. This may be 


expressed as: 


ae a 3 ele 
Ur 50,1 7 28, (3 ) 





eg 


where i) 


J 
b6. 
J 


total rotation of joint j 


rotation of bar j. 


Substituting Equations (3.2.9) into Equation (3.2.10) and intro- 
ducing simplified notation give the desired joint rotation relationships. 
iy Wageae > 5s 
Be Ash + ALB, (3,2. 3 
UL = h, - Ah 
oe Sk 5 


where A, = =(2/a + 1/c + b/ad) 


eS 
{ 


2. 
| 


A, = i/c 


A, = =(1/c - b/ad) 


35.5 Moment-Rotation Relationships for Joints 

As the arch displaces laterally under load, bending moments are 
developed in the arch which act to resist the displacement. Likewise in the 
model, as a change occurs in the angie between the bars meeting at a joint, 

a resisting moment is developed in the joint. In Section 3.2.3 expressions 
are derived relating joint rotations and model displacements. It is the pur- 
pose of this section to derive relationships between these joint rotations and 
the corresponding moment resistance developed by the joints. 

If bending in the arch is assumed to be uniform over a region such 
as that shown in Figure 3.4, the moment-rotation relationship for the joint 
may be derived as follows: 

Let M, = moment resistance of the arch rib cross section 


J 


E. = modulus of elasticity of the material from which 
J the arch rib is made 


L. = moment of inertia of the arch rib cross section 


L. = length of bar (assumed of equal length on either 
J side of joint ) 


¥, = joint rotation 
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Then, from Figure 3.4: 


¥ M, L 


cs epee ws 
2 EL 2 
Jd 
Hence: 
EI 
M - ody 
J a J 
For an arch of constant cross section, my and may be considered 
constant. ne may be chosen as a constant when setting up the model. For these 
conditions it is convenient to simplify the notation as: 
M, = ¥ ie | 
; Ss, ; (3.3.1) 
where Ke r= EI /L, = bending stiffness of the arch rib as concentrated 
5 J at joint j of the model. 


Strictly speaking, Equation (3.3.1) is valid only in the elastic 
range. Since it is intended to permit plastic deformation of the arch rib, 
the moment-rotation relationship is extended as shown in idealized form in 
Figure 3.5. This follows directly from Equation (3.3.1) and the requirement 
that the arch rib possess sufficient ductility to continue to rotate without 
rupture or decreasing moment resistance once the maximum moment resistance is 
developed. This can be achieved readily in the design of the arch rib. 

The limits of this moment-rotation diagram for any particular case 
will thus be determined by the strength, ductility and geometry of the arch 
rib cross section as well as by the geometry of the arch itself. As mentioned 
previously (cf. Section 2.1), no modification of bending resistance is made 
for thrust in the arch caused by the blast loading, or for dead and live 
loads existing on the arch prior to the arrival of the shock wave. In general, 
these effects will tend to reduce the stiffness of the arch. Again, while 
the arch must be capable of withstanding these loads, their net effect on the 


total resistance of the arch-soil "“structure"to lateral deformation is 
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relatively small and such modification is not considered warranted for this 


study. 


3.4 Development of Equivalent Arch Resistance 


In order that the moment resistance of the joints may be inserted 
conveniently into the equations of motion for the model, relationships are 
needed expressing the resistance of the arch to lateral deformation. These 
are derived in this section first in general terms and then specialized for 
the choice of independent displacement variables used in this thesis, viz. 
the horizontal displacements of the haunches. 

In the sign convention used for this derivation clockwise shears are 
considered positive. Bending moments are considered positive when they tend 
to produce compression in the outer fibers of the arch rib. Bar rotations 
are measured at the left ends of bars, counterclockwise from a horizontal 
exis directed toward the right. As stated in Section 1.2, loads and displace- 
ments are taken as positive horizontally to the right and vertically upward. 

Referring to Figure 3.6 and considering the equilibrium of either 


one of the bars, it can be seen that: 


M. - M 
g - wsdl (3.509 
3 L. 
J 
where = shear in bar j 
u. = resisting moment in joint j 
oe = length of bar j 


Summing horizontal forces acting on the free body shown in 


Figure 3.7 gives: 


Q, = S, sin (6, 7 60, ) - Ss ) (3.4.2) 


, : sin(6,.) + 86 


j+2 jt+1 
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where Q, = horizontal component of arch resistance at joint j 
2 = shear in bar j 
oe = inclination of bar j 
eo rotation of bar j under nuclear blast loading. 


Summing vertical forces acting on the free body shown in Figure 3.7 
| gives: 


Q 545 = - S, cos (6. 4 66) fe Sa cos (6.4 rs 86.4) (3.4.3) 


where Q = vertical component of the arch resistance of joint j. 


ny 

In deriving the equivalent arch resistance at the haunches, it is 
convenient to express the inclination of the shearing forces in the bars in 
terms of the arch geometry. This is readily accomplished using the familiar 
trigonometric identity stated in Equation (3.4.4). 


sin(6@ 56.) = sin 8, cos 56, + cos @. sin 56. 4.4 
(ORD EEN) COI 3 j j Tad 


Using the notation of Figure 2.1], and limiting consideration to small rotations, 


one can write: 


sin(6, a 60, ) = [a - (2)ny] 

sin(@, + 66,) = z [ce + =(a/c A b/a)h, - =(a/c - b/a)h, } aed 
sin(6, “ 665) = pt - =(a/c - b/a)h, + =(a/c an b/a)hs] 

pin(@), + 66) ) = i, [-a - (b/a)h, ] 


It is now possible to write expressions for the equivalent arch 
| resistance in the horizontal direction at the haunches by substituting 


Equations (3.4.1) and (3.4.5) into Equation (3.4.2). ‘Thus: 








where 


5D 


ao [eG j- [Agee debs as | 
2 
(4.46) 
-M 
5 = [Ape] HE HE Bos] ecd [= 


a and a, = horizontal components of arch resistance at joints 1 
and 3 respectively 


M > M,) ms = resisting moments developed in joints 1, 2 and 3 
respectively 


a,b,c,d = generalized arch dimensions as given in Figure 2.1] 


h,, h, = horizontal components of displacement at joints l 


i? 5 
and 5 respectively. 


By limiting consideration to equal bar lengths, rearranging terms 


and simplifying notation, one can write the equivalent arch resistances as: 


where 


5.5 


M,(q, + a,b) - abs) - Mota, + dqby - ahs) 


“My(q, + a)h) - ahs) + MQ(a, + ah) - aghs) 


“A 
2 


(3.4.7) 


a, (1/L°)(a + c) 

dy = (A/eL*)(a/e - d/a) 

a = (c/E") 

qq = (1/2L°)(a/e + b/a) 
Rib Shortening 


The results of the dynamic analysis developed elsewhere in this 


thesis may be modified using the "rib-shortening’ correction derived in this 


section in order to give an indication of the significance of the effect of 


direct compression in the arch rib on the total displacement of the arch. The 


necessity of considering direct compression effects separately was discussed 


earlier (cf. Sections 2.1 and 3.1). Essentially, the model cannot accept axial 
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changes in the length of the bars and retain its characteristics as a two- 
degree-of-freedom system. While this is a definite limitation, it is not 
serious in a study of lateral deformations as long as the arch rib has suf- 
ficient capacity to withstand the thrust in the arch without general yielding 
or buckling. 

The approach used in this section is to consider the response of 
the arch to a uniform compression loading. This response is then superimposed 
on the lateral response in order to approximate the total response of the arch. 
This is an approximation in that these effects are in fact interrelated and 
cannot be entirely separated, although if one considers the relative compres- 
sibilities of the arch and the soil, direct compression effects are not as 
important as might otherwise be expected (cf. Chapter IX). 

Note that the soil provides a high degree of damping to the response 
of the arch in this mode. This effect cannot be neglected as it tends to pro- 
duce an over-critically damped response. Further, the extremely short natural 
period of vibration of the arch-soil "structure" in this mode makes the rise 
time of the shock wave in soil significant. This also tends to reduce the 
magnitude of the response. Therefore, for the purpose of the approximate 
correction derived in this section, a dynamic magnification factor of unity 
is assumed for this motion. 

The arch loading does contain a significant "uniform compression’ 
component because of the rapid envelopement of the arch by the shock wave 
and the relatively slow decay of overpressure with time behind the shock 
front. Therefore, it appears reasonable to consider this component of the 
loading in computing the thrust in the arch rib. In this study, the uniform 
compression loading on the arch is taken as a proportion of the average vertical 


overpressure acting on the arch (cf. Chapter VI). This proportion depends 
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principally upon the ratio of the horizontal to the vertical overpressure 

in the soil, although the fact that the crown is subjected to a higher value of 
overpressure than the average radial overpressure on the arch must also be 
considered. 

Two additional assumptions are made in order to simplify the compu- 
tation of the rib shortening itself. First, the average thrust in the arch is 
assumed to be equal to that in a circular ring under uniform compression. 
Second, the arch is assumed to maintain its same general shape as it deforms 
under the uniform compressive loading. This is approximate, but then the 
nature of the correction derived in this section does not warrant a rigorous 
solution for rib shortening in the arch. As will be seen, these assumptions 
lead to simple expressions which give reasonable results and are therefore 
justified on this basis. 

Consider the arch to be loaded as shown in Figure 3.8. The thrust 
in the arch, T, may then be determined by equating the sum of the vertical 
forces acting on the arch to zero. 


OTsin B = ORsin B pt 


= i 
T = P, R 
Tsp! B (3.5.2) 
a a 6 e 
as Ty " 
where p. = uniform compression component of the average overpressure 


loading the arch (cf. Chapter VI). 
The change in length of any one of the bars of the model caused by 


the thrust is given by Equation (3.5.2). 


OL, = TL ,/AE 
OL, = = p! BL ,/ AE (3.522) 
where OL, = change in length of bar j due to thrust 
L, = length of bar j 
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A cross-sectional area of the arch rib (assumed constant) 


E 


modulus of elasticity of the arch material. 

The model is shown in its undeformed and deformed configurations in 
Figure 3.9. The arch abutments are considered immovable (cf. Section 1.2); 
therefore, joint O remains stationary. For the model to maintain the same 
general shape as the bars shorten, let bar 2 rotate about point 0' such that 
joint 1 moves in a direction perpendicular to the original direction of bar é. 
This is shown in Figure 3.10. Now, from symmetry joint 2 must move vertically 
downward. Therefore, joint 2 may be considered first to perform the same 
motion as joint 1 and then to rotate and shorten until it is vertically below 
its original position. This is illustrated in Figure 3.11. 

The displacement coordinates of joint 1, i.e. the "rib-shortening” 
correction for joint 1, may be written directly from Figure 3.10. 


by 8L, /sin(6, -6,) 


8L, /sin @.cos @.-cos @.sin @ 


i 2 1 2 
= SL, (L)*/ad-be (3.5 98m 
bh, = - 5)sin 6, 
=e (OL, )Le/ad-be (3.5009 


oh), = + 0, COB 6, 


+ (OL, )La/ad-be (3.5.59 


The displacement coordinates for the crown may be written directly from 
Figure 3.11. 


bh. = 0 (3.5.69 


eae oL,,/sin 0, 


(SL, )L/c (3.5.7) 
The "rib-shortening’ corrections for the model in terms of the 


arch geometry and the uniform overpressure which is effective on the arch 





DT 


may now be obtained by substituting Equation (3.5.2) into Equations (3.5.4) 


through (3.5.7) and considering symmetry. 


2 
dh See BL c 
1 cra AE(ad-bc ) 
dh, = 0 
oh, = - Phy 
l1- seal 

: = — t 
ony = + 3 Pe | 6-5-8 

; oe 
| 1 —: | BL | 
ES ae 
Ohy = oh), 


Note that if the overpressure, Py? is expressed in units of force/(length)*, 
the cross-sectional area of the arch rib, A, must be expressed in units of 
Sreneth )-/iength. 

An alternate approach to the derivation of the "rib-shortening™ 
corrections is to allow the haunches and the crown to move such that, if the 
elevation of the abutments is used as a reference, the ratio of the elevation 
of the haunches to the elevation of the crown remains constant. Expressions 
were derived on this basis and found to give results, for the numerical 
examples studied herein (cf. Chapter VIII), which do not differ significantly 
from those given by Equations (4.5.8). Since the alternate expressions are 


somewhat more involved, they were discarded in favor of Equations (3.5.8). 
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CONSTRAINT LOADINGS 
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Figure 3.1. Constraint Loadings 
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Figure 3.2. Load Concentrations on Model 





Figure 5.5. Diagram Used in Determining 
Joint Rotation Relationships 
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Figure 3.4. Diagram Used in Determining Moment-Rotation 
Relationships for Joints 
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Figure 3.5. Idealized Moment-Rotation Diagram for Joints 
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Figure 3.6. Free-Body Diagram of Typical Bars and Joint 
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Figure 3.7. Free-Body Diagram of Typical Joint 
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Figure 4.6. Uniform Radial Load on Arch 
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Figure 3.9. Deformation of Model Due to Rib Shortening 
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Figure 3.10. Diagram Used in Determining Rib- 
Shortening Correction at Haunch 





Figure 3.11. Diagram Used in Determining Rib- 
Shortening Correction at Crown 
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CHAPTER IV 


SOIL RESISTANCE 


In the case of an underground arch the soil contributes significant- 
ly to the total resistance to nuclear blast loading of the arch-soil “structure”. 
As with the arch, this resistance can be conveniently considered in two parts, 
viz. inertial resistance, or that due to the mass of the soil which moves with 
the arch, and deformation resistance, or that due to the presence of shear and 
frictional forces within that portion of the soil which deforms with the arch 
under blast loading. Deformation resistance is discussed here. Inertial or 
mass effects are considered separately in Chapter V. 

As discussed in Section 2.2, the concept developed in this study 
regarding soil resistance due to frictional and shearing forces in the soil is 
based upon consideration of the state of stress existing in the soil at any 
particular instant during the response of the structure. Specifically, in the 
derivations which follow, the limiting stress conditions of the "Active State 
of Plastic Equilibrium’ and the "Passive State of Plastic Equilibrium," as 
defined in Section ©.2, are considered as well as the intermediate condition 
of "Earth Pressure at Rest." In general, differences in the gravity forces 
acting on the arch for different states of stress in the soil, together with 
the difference between the blast loading exerted on the ground surface and 
that transmitted to the arch by the soil, are considered as due to "soil 
resistance." A linear interpolation is used between the "at rest" and the 
"plastic equilibrium’ states, and an unchanging soil resistance assumed for 
deformations greater than those required for "plastic equilibrium." This in 
effect assumes an elasto-plastic resistance function for the soll which is 


consistent with the present state of knowledge of soil-structure behavior. 
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Since the model developed for the dynamic analysis (cf. Section br 
contains three interior hinges, three regions, one adjacent to each hinge, are 
considered separately from the standpoint of soil deformation and resistance. 
Conditions are different at the haunches from those at the crown. Further, 
since conditions at a haunch are different depending upon whether the haunch 
has moved inward into the arch or outward into the soil, these cases are treated 
separately. For convenience, a haunch which has moved inward into the arch with 
respect to its original position prior to the arrival of the shock wave, or with 
respect to its original position as modified by any "permanent set’ in the soil 
adjacent to the arch, is referred to herein as an "interior" haunch. A haunch 
which has moved outward into the soil with respect to the same reference posi- 
tion is referred to as an “exterior” haunch. This classification is particular- 
ly useful when studying "rebound," or response of the system after the first 


maximum response has occurred (cf. Section 4.4). 


4.1 Interior Haunch 

As a haunch moves inward into the arch, the soil tends to follow it 
under the influence of the overpressure acting on the ground surface. This 
tends to create an "active state of stress within the soil surrounding the 
haunch. Further, since the blast loading is applied to the arch through the 
soil, the shearing and frictional stresses developed within the soil mass will 
assist the structure in resisting the load. 

To illustrate how this "soil resistance’ is developed, consider the 
loading "sensed" by the arch as the arch deforms. As mentioned in Section 2.2, 
the pressure surface, or surface on which the soil pressure acts, at the side 
of the arch is assumed to coincide with the lower bar of the model and in 
addition to extend upward in the same direction from the haunch to the ground 


surface. In Figure 4.1 P' represents the force on the side of the arch due to 
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the weight of the soil. It is assumed to act at the lower third point of the 
pressure surface. P' represents the force on the side of the arch due to the 
"surcharge of the blast overpressure acting on the ground surface. It is 
assumed to act at the midheight of the pressure surface. Since classical 
earth pressure theory (Coulomb assumptions, cf. Section 2.2)is being used to 
define the limiting stress conditions for the calculation of hydrostatic pres- 
sure ratios for the soil stresses, P' and P’ are assumed to act at an angle to 
the direction of the pressure surface. This angle, 5, is the friction angle 
corresponding to the friction developed between the surface of the arch and 
the soil. 

Now consider the hinges of the model to be "locked" so that when the 
blast loading is applied, no deformation can take place. If no movement takes 
place, the state of stress in the soil cannot change, and it therefore remains 
in the "at rest" condition. For this condition of "no deformation’ one can 


write expressions for the forces acting on the sides of the arch as: 


rer e 
oa ™ 3 ato Son (4.12.1) 
where POA = component in the direction of Pi of the force acting on the 
side of the arch due to the weight of the soil and corre- 
sponding to the condition of earth pressure at rest 
Us = weight of soil per unit volume 
qd, = total depth from the ground surface to the arch abutment 


Kon = ratio of the soil pressure in the direction of Pi to the 
vertical soil pressure for the condition of earth pressure 


at rest, 
: a es 
and ; Poa = Pe jt Kon (4.1.2 
where Pon = component in the direction of P, of the force acting on the 


side of the arch due to the blast overpressure and correspond- 
ing to the condition of earth pressure at rest 








Pas = average overpressure acting at the ground surface on the 
J potential “active failure wedge in the soil (defined 
below). Note that since p., acts downward, it is considered 


negative in accordance with the sign convention used herein 
(cf. Section 1.2). 


Now consider the hinges of the model to be "unlocked" and the 
model permitted to deform until a fully developed "Active State of Plastic 
Equilibrium’ is reached. At this point frictional and shearing forces exist 
within the soil mass and the expressions for the forces acting on the side of 
the arch are: 
ak (4.1.3) 


where Pa = force acting on the side of the arch due to the weight of the 
soil and corresponding to a fully developed “active state of 
stress within the soil mass. 


K, = hydrostatic pressure ratio for the ‘active’ state of stress 


‘i in the soil, 
and: P" =+p .aK (4.1.4) 
A sj o A 
where eA = force acting on the side of the arch due to blast overpressure 


and corresponding to a fully developed “active state of stress 
within the soil mass. 


The results of the preceding discussion are presented graphically 
in Figure 4.2. The "no deformation’ load components, Equations (4.1.1) and 
(4.1.2), are plotted along the vertical axis, i.e. corresponding to no rota- 
tion of the pressure surface. The "fully active" load components, Eyuations 
(4.1.3) and (4.1.4), are plotted at a deformation Ee corresponding to that 
amount of rotation of the pressure surface required to reach the "Active 
State of Plastic Equilibrium’ within the soil mass. A discussion of the 
amount of deformation required is contained in Sections 2.2 and 4.4. 

Because of the forces developed within the soil during deformation, 
Pi’ + P, is less than P/, + one This difference is defined herein as the 


A A OA 
"soil resistance” corresponding to an inward movement of a haunch sufficient 
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to develop fully an "active' state of stress in the adjacent soil. An 
expression for the soil resistance corresponding to this "fully active" condi- 
tion can now be derived by taking the difference of each of the load components 
considered, i.e.: 


Let: R, and Ra = maximum values of soil resistance corresponding to a 
fully developed "active' state of stress in the soil. 


Then, using Equations (4.1.1) and (4.1.3): 


SN oy ea 
= 57685 Non - 3 7040 Ka 
Ri = = 7d (Ky, - K,) (4.1.5) 
and, using Equations (4.1.2) and (4.1.4): 
nue onan 
om Py 8Ko, ¥ Pe jFoka 
Ra = - B, 54, (Ko, - K,) (Ho Gn 


The nature of the variation of the soil resistance between the 
"at rest and the "fully active’ states can only be inferred from field 
observations of test structures, laboratory model studies and a general con- 
Sideration of the nature of stresses in masses of soil. No exact relationship 
is known to the writer. One possible relationship is shown qualitatively in 
Figure 4.3. This corresponds, in general, to Figure 2.2. For the purposes 
of this study the relationships between deformation and soil resistance will 
be approximated using linear functions. This is discussed in detail in 
Section 4.4, 

While Figures 4.2 and 4.3 are helpful in visualizing the development 
of the soil resistance concept used in this study, they are somewhat misleading 
in that the two components of soil resistance as expressed by Equations (4.1.5) 
and (4.1.6) cannot be added directly to one another for insertion in the 


' equations of motion for the model because they act at different points. In 
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order to combine these quantities into what will be termed herein the 
"effective soil resistance at the haunch, the lower bar is considered as a 
beam simply supported at the haunch and at the abutment and "loaded" with the 
two soil resistances Rj, and Ra The "reaction" at the haunch is then taken 
as the soil resistance which is "effective’ at the haunch. This concept is 
illustrated in Figure 4.4. While this neglects moments at the haunch, it 
is sufficiently accursete for this study. 
Also shown in Figure 4.4 is the soil wedge corresponding to the 
development of the ‘active" state of stress in the soil. This wedge is 
bounded by the ground surface, the pressure surface and a plane passing upward 
through the arch abutment at an angle of (+ z) with respect to the horizontal, 
where ? is the angle of internal friction of the soil. This concept follows 
directly from the use of classical earth pressure theory to define the "active" 
state of stress in the soil plus the assumption of a plane failure surface for 
the determination of Ka As with the arch itself, this wedge is considered 
herein to have a unit thickness in the direction of the longitudinal arch axis. 
The overpressure term, Poy in Equations (4.1.2), (4.1.4) and (4.1.6) is taken 
as the average overpressure acting at the ground surface on this failure wedge. 
Keeping in mind that Rj, and Ry as given by Equations (4.1.5) and 


(4.1.6) are maximum values, the maximum effective soil resistance for an 


interior haunch may be obtained directly from Figure 4.4 as: 


0. = = cos(w + &) (RA (a,/3) ne RA (a,/2)) 
Q5 j = 7 (QR) je aha cos (w + 8) (hole 
A _ bd, | ba, 
and : Qs 543 2 So sin(w + 8) RA (3) “ R, (2) | 
a 
> Oo ' rT) 
Qe 543 Se (2R) NS 5R,) sin(w + 5) (4.1.83 





pp: 
A 
where ae = maximum value of the horizontal component of the soil 
resistance which may be considered effective at an 
interior haunch 
A 
5 543 = maximum value of the vertical component of the soil 


resistance which may be considered effective at an interior 
haunch 


a,b,c,d = generalized arch dimensions (cf. Figure 2.1) 


Ww = inclination of the lower bar of the model with respect to 
the vertical 


6 = friction angle corresponding to the friction developed 
between the surface of the arch and the soil. 


d 


5 depth from the ground surface to the arch abutment 


Ryka = maximum values of soil resistance components as defined by 
Equations (4.1.5) and (4.1.6) 


Substituting Equations (4.1.5) and 4.1.6) into Equations (4.1.7) 


and (4.1.8): 


2 
d 
Q.5 = Ze (Kog - K,M(7,4, - 3B, Jeos(w + 6) (4.1.9) 
2 
© “0 = § h 
Qe 543 = - Zz (Ky, = K, (7,4, - 5p, ;sin(u 0) (4.1.10) 


These equations will be used in Section 4.4 in developing expressions for the 
total soil resistance of the system. 

One further point to be considered here involves the hydrostatic 
pressure ratios for this case. Methods for determining Ky are compared in 
Appendix A. Values for Kon cannot be calculated using classical earth pressure 


(10) In the numerical 


theory, although values have been measured in tests. 
solutions of this study an observed value of the ratio of horizontal to 
vertical soil pressures corresponding to the "at rest" condition is assumed. 
It is then necessary to "resolve" this ratio from that corresponding to hori- 
zontal pressures to that corresponding to pressures parallel to the line of 


action of Pe A diagram of the relationships involved is given in Figure 4.5. 
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Referring to this figure, one may write: 


al 
raat a® K 
= pe ee r x? 2 
_ From which: 
v 2 eal 
Kon =\/1 + KG COs | ( + 5) - arctan . (4.1.11) 
where kK, = observed value of the ratio between the horizontal and 
vertical soil pressures for the condition of earth pressure 
at rest. 


4.2 Exterior Haunch 

The soil resistance developed at an exterior haunch may be derived 
in a manner analogous to that used in the previous section for the interior 
haunch. The essential difference is that movement of the arch outward against 
the soil tends to produce a “Passive State of Plastic Equilibrium" (cf. Sec- 
tion 2.2) within the soil. As in the previous case, the "soil resistance” 
will be taken as the differences in the loads transmitted to the arch by the 
soil in the "locked" and "unlocked" conditions. The forces acting on the arch 
for this case are shown in Figure 4.6. 

If the arch is now considered to be "locked" so that no deformation 
can take place, hence no frictional forces develop in the soil, and the blast 


loading applied, one can write: 


r ii e 
Pop = 3 7545 Kop (4.2.1) 
where Pop = component in the direction of Ps of the force acting on the 
side of the arch due to the weight of the soil and correspond- 
ing to the condition of earth pressure at rest, 
Kop = ratio of the soil pressure in the direction of Pl to the 


vertical soil pressure for the condition of earth pressure 
at rest, 





5) 
and: Ee Ps s¢4oKop (hk senen 


where Pp... = component in the direction of P, of the force acting on the 


side of the arch due to the blast overpressure and corre- 
sponding to the condition of earth pressure at rest, 


p', = average Overpressure acting at the ground surface on the 


potential "passive" soil failure wedge (defined below). 
Note that since p'. acts downward, it is considered negative 


in accordance with’ the sign convention used herein (cf. 
Section 1.2). 


If the structure is now "unlocked" and permitted to deform until the 
"Passive State of Plastic Equilibrium’ is reached, frictional forces will 


develop within the soil mass and the forces acting directly on the side of the 


arch will be: 
ee e ; 
where Po = force acting on the side of the arch due to the weight of the 


soil and corresponding to a fully developed "passive" state 
of stress within the soil, 


hydrostatic pressure ratio for the “passive” state of stress in 


ws 


in the soil, 
. ass 5 
and: Eee Ps aK, (ho.4) 
where P = force acting on the side of the arch due to the blast 


Overpressure and corresponding to a fully developed "passive" 
state of stress within the soil. 


These load components are plotted in Figure 4.7. The "locked" 
condition is represented by TF = 0, i.e. no rotation of the pressure surface. 
The “unlocked” condition is represented by Shy? i.e. sufficient rotation of 
the pressure surface to develop fully the "Passive State of Plastic Equilibrium’ 
within the soil (cf. Section 2.2). 

As in the case of the interior haunch, the differences between the 
"locked" and "unlocked" conditions are taken as the "soil resistance." 

Referring to Figure 4.7, these differences may be expressed by subtracting 


Eyuations (4.2.1) and (4.2.2) from Equations (4.%.3) and (4.2.4). 
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Thus, if: Ry and R, = maximum values of soil resistance corresponding to a 
fully developed "passive" state of stress in the soil, 


Then: Rp = Po - Pop 
_ L 
=u 7 5doKp “3 Haake 
a. 72, 1 - X,) (4.2.5) 


- a i . ft 
And: Rp - Ps Pop 


Be one a ca O Kop 
Ry = - Pe d(K, - (4.2.6) 
Since the load components P' and P" are assumed to act at different 
points on the pressure surface (cf. Figure 4.6), the corresponding soil re- 
sistances R, and R, must be combined in a similar manner to that used in the 
case of the interior haunch. Again considering the lower bar of the model as 
a simply supported beam acted upon by the soil resistances Rj and Rp» the 


effective soil resistance components at the haunch are (cf. Figure 4.8): 


‘ da d 
os(d - w) E (ee (2) 
ge (se + 5R,) cos(5 - w) (h.2.99 


© | 


a 
Sy 


. : bd, | bas 
And: Og 7 7p tinld - 6) [Rp (<2) + RY (5) | 
a 


- er + 3Ry) sin(S - w) Glue o.8) 


where a = maximum value of the horizontal component of the soil 
resistance which may be considered effective at an 
exterior haunch, 


= maximum value of the vertical component of the soil 
resistance which may be considered effective at an 
exterior haunch, 


8343 


a,b,c,d = generalized arch dimensions (cf. Figure 2.1), 





92 


W = inclination of the lower bar of the model with respect 
to the vertical, 
© = friction angle corresponding to the friction developed 


between the surface of the arch and the soil. 
Substituting Equations (4.2.5) and (4.2.6) into Equations (4.2.7) 


and (4.2.8): 


2 
a 
a - rae Oe - Kop) (7.4, - 5p, ,)cos (8 - w) (4.2.9) 
of 
a = - z (K, - Kp)(7,4, - 3pi,)sin(6 - w) (4.2.10) 


Note that the soil failure wedge now considered (cf. Figure 4.8) is 
bounded by the ground surface, the pressure surface and a plane through the 
arch abutment inclined at an angle of (F - z) with respect to the horizontal, 
where ? is the angle of internal friction of the soil. A unit thickness of 
wedge parallel to the horizontal arch axis is considered in this thesis. The 
assumptions leading to this wedge are the same as those leading to the defini- 
tion of the "active" soil failure wedge described in Section 4.1, i.e. the 
use of classical earth pressure theory (Coulomb's assumptions) and the assump- 
tion of a plane failure surface in the determination of the hydrostatic pressure 
ratio corresponding to full development of the "passive" state of stress within 
the soil (cf. Appendix A). 

Equations (4.2.9) and (4.2.10) are used in Section 4.4 in deriving 
expressions for the total soil resistance of the system. In order to use 
these equations, one must be able to calculate the hydrostatic pressure ratios 
for the soil pressures corresponding to the ‘at rest" and "passive" states of 
stress. Methods of finding Kp are compared in Appendix A. Kop is found by 
"resolving" observed hydrostatic pressure ratios relating horizontal to 
vertical soil pressures to a ratio relating soil pressures acting parallel to 


the line of action of P! 


p to vertical soil pressures. The technique used is 
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analogous to that described in Section 4.1 for finding Kon: 


relationships involved is given in Figure 4.9. Referring to this figure, one 


A diagram of the 


may write: 
Ay 2 
oa i a6 eS 
18 2 OS eee [arctan ay (5 - ») 
i Pian ay Se 
a 7545, 1+ K, 
a it 
From which: Kop Wa +k cos [arctan eas (5 - ») | (4.2.10) 
O 
where k = observed value of the ratio between the horizontal and vertical 


s0il pressures for the condition of earth pressure at rest. 


4.3 Crown 

As the soil-arch "structure" deforms in response to a blast loading, 
shear and frictional forces are developed in the soil near the crown. The 
definition of these forces is not an easy matter. In evaluating the relative 
importance of these forces, consider that for the shallow depths of cover 
included in this study (cf. Section 1.1) the volume of soil at and adjacent to 
the crown is relatively much smaller than the volume of the soil wedges defined 
in the preceding sections, affected by lateral displacements of the sides of 
the arch. Thus, the maximum resultant value of these soil forces will 
necessarily be relatively smaller at the crown than at the sides of the arch 
because of the smaller volume of soil involved. Further, some of these 
resisting forces serve only to distribute the applied load away from the 
crown to other regions of the arch, thus contributing no net resistance to the 
system. A portion of these forces can, however, be assumed logically to act 
to transmit some of the applied load away from the arch. This latter portion 
of the soil forces does assist the arch in resisting an applied load but is 


not readily defined explicitly. 
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An attempt to include at least part of this "effective" portion of 
the soil resistance near the crown in the total soil resistance of the system 
is made by continuing upward the pressure surfaces at the sides of the arch, 
in the same direction as the lower bars of the model, until they intersect 
the ground surface. This is not a rigorous approach, but it is believed to 
be sufficiently accurate for the present. Since the forces involved are 
relatively small, it would seem to follow that their influence on the response 
of the system as a whole should be relatively small as well. 

The effect of the inertial resistance of the soil near the crown 
is discussed in Chapter V. 

As a result of the foregoing discussion, the effective soil resistance 
to be placed on the model at the crown is zero. Thus, the following expressions 


may be written: 


Q.5 = 0 (4.3.1) 
where Q50 = horizontal component of the effective soil resistance at 
the crown, 
and : 
2,5 7 0 (4.35.2) 
where Be = vertical component of the effective soil resistance at the 


crown. 


4.4 Development of Equivalent Soil Resistance 
Since the model is a two-degree-of-freedom system (cf. Section 3.1), 


the effective soil resistances derived in the preceding sections must be 
combined such that their net, or "equivalent," effect on the model is defined 
in the direction of the independent displacement coordinates. This is done 
in this section for the set of independent coordinates already selected, viz. 


the horizontal displacements of the haunches. Thus, the object here is to 
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define the resultant resistance offered by the soil to horizontal movements 
of the haunches. Note that this is not the same as the soil resistance 
components derived in Sections 4.1 and 4.2 since now the entire system is 
considered simultaneously. 

The technique used herein is simply to apply the effective soil 
resistance components to the model as a system of loads. Then, by applying 
the "Load Relationships’ derived in Section 3.2.1 (Equations (3.2.5)), the 
equivalent soil resistances in the horizontal direction at the haunches are 
determined. This may be expressed algebraically by substituting the effective 


soil resistance components, for the load components, aa? in Equations 


Q, 3? 
(3.2.5). Thus: 


Q64 F 851 _ me 7 aa) 250 - (2) 854 7 ee = Sa 
(4.4.1) 
Qy5 = 5 (= GR) Got Qs - 5 (= B) Me + CE) Me 
where Q.4 = equivalent soil resistance at the windward haunch, 
Q. = equivalent soil resistance at the leeward haunch. 


Since the expressions for the components of the effective soil 
resistance at the joints were derived in the preceding sectinns in terms of 
the limiting or maximum values, it is convenient here to write maximum values 
for the equivalent soil resistance. Further, since the haunches may be 
displaced either inward or outward, it will be convenient to have two sets of 
equations, one for interior haunches and one for exterior haunches. These 
equations may then be used in various combinations for dynamic analysis as 
the configuration of the system dictates. 

For the case of both haunches displaced inward, the maximum equiva- 
lent soil resistance may be found by substituting Equations (4.1.9), (4.1.10), 


(4.3.1) and (4.5.2) into Equations (4.4.1). This gives: 








Note that 


Thus, one 


where 


AS 
Q 


S 


{> 


65 = 


the 


may 


| > 


{> 


p 


Sj 


BJ 


SJ 


Be, 


qc 


1° A (Ko, - K, (7,4, - 3B.) | cos + 5) + (2)sin(w i‘ 5) | 


a b 
Zo (Koy ~ K)(7g4, ~ 38,5) | eos(w + 8) + (Boin(w + 8)] 


preceding expressions differ only in the overpressure tern. 


write, in general: 


qc 
O 


= : b 
Z (Koy ~ KM(7gd, ~ 38,5) [cos(w + 0) + B)ein(u + 8) | 
(eh. oF 
maximum value of the horizontal component of the eyuivalent 
soil resistance at haunch j when the haunch is displaced 


inward, 


average overpressure acting at the ground surface on the 
"active soil failure wedge adjacent to joint j (cf. Figure 4.4). 


Similar expressions may be written in the case where both haunches 


are displaced outward. Substituting E,uations (4.2.9), (4.2.10), (4.3.1) and 


(4.3..) into Equations (4.4.1): 


A ae i . 

C1 = 5 (K, - Kop )(7,.4, = 3P.4) | cost - 5) + (~)sin(w 2 °) | 
A ac 

Q! 2 Saas ; D 

aa = Z_ (Kp ~ Kop) (7.4, - 3g) cos (u = Bb) + (~)sin(w 7 8) | 


Note that these 


in general: 
a 
| 

wee 
fd 
q 

where 8 J 
4 

Pej 


expressions also differ only in the overpressure term. Thus, 


a5 7 b 
Z2 (Ky - Kyp)(748, ~ 3544) [eos(w - 8) + (Boin(w - 8)] 


(4.4.3) 


= maximum value of the horizontal component of the equivalent 


soil resistance at haunch j when the haunch is displaced 
outward, 


average overpressure acting at the ground surface on the 
“passive” soil failure wedge adjacent to joint j 
ter. Figure 4.8). 
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Note that both Equation (4.4.2) and Equation (4.4.4) express the 
equivalent soil resistance at a haunch solely in terms of the components of 
the effective soil resistance acting at that haunch. This permits the haunches 
to be considered separately and reduces the number of cases which must be con- 
sidered in the remainder of this discussion to two, viz. the interior haunch 
and the exterior haunch. The resulting expressions can then be applied to any 
conceivable configuration of the systen. 

An important aspect of the soil resistance developed as the system 
responds to a dynamic load involves the amounts of deformation required to 
develop fully the "active" and "passive" states of plastic equilibrium. In 


(10,11,12,13) 55 these 


Section 2.c reference was made to a number of estimates 
parameters. There is general agreement that in the case of retaining walls 
backfilled with cohesionless soil, a rotation of the wall about its base 
corresponding to a movement away from the backfill at the top on the order 
of 0.1% of the wall height is sufficient to reduce the state of stress in the 
soil behind the wall from the “at rest’ to the "fully active" state. In the 
case of the development of the "passive" state of stress, quantitative esti- 
Mates are generally lacking, although it appears certain that the required 
movements are greater than for the "active" case (cf. Figure 2.2). A 
technique for estimating this required movement is outlined below. Note that 
the pressure surfaces assumed in deriving the effective soil resistances 
(cf. Section 4.1 and 4.2) are similar to the pressure surfaces behind retaining 
walls. Therefore, the significant deformations of the model for the develop- 
ment of soil resistance are the rotations of the lower bars, since these define 
the positions of the pressure surfaces. 

In this study the idealized soil resistance-displacement function is 
obtained by assuming a linear buildup of soil resistance with increasing dis- 


placement until a maximum value is reached. The soil resistance is then 
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assumed to remain at this maximum value for further displacements in the same 
direction. This results in a so-called "elasto-plastic" resistance diagram 
as illustrated in Figure 4.10. 

Thus, the soil resistance-displacement function for an interior 
haunch may be defined by using Equation (4.4.2) to find the maximum value of 
the soil resistance, and the assumption that this value of soil resistance 
corresponds to an inward rotation of the lower bar, op of 0.001. This is 
illustrated in Figure 4.11. The range of possible values for the soil resist- 
ance corresponding to each value of displacement results from the overpressure 
term, Poy? in Equation (4.4.2). 

The maximum value of soil resistance for an exterior haunch is given 
by Equation (4.4.3). Since quantitative estimates of the corresponding 
deformations are lacking, the assumption is made that the slope of the idealized 
soil resistance-displacement function in the elastic range is constant for any 
given value of overpressure acting on the soil failure wedges. For the specific 
examples studied herein, this deformation corresponds to a rotation of the 
pressure surface on the order of 0.005, which appears to be in reasonable agree- 


ment with existing theory and experimental evidence. (1013) 


Note, however, 
that this value has not yet been determined definitely. An experimental in- 
vestigation of the deformation corresponding to full development of the 
“passive” state of stress in the soil is certainly needed (cf. Chapter IX). A 
complete, idealized soil resistance-displacement diagram is shown in 
Figure 4.12. 

Expressions for the "yield" deformations of the soil may be derived 


as follows: 


Let: at = rotation of the lower bar of the model as the 
structure responds to a blast loading, 
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! rotation of the lower bar of the model corresponding 
Y to full development of the "active" state of stress 
within the soil adjacent to the haunch, 


g 


rotation of the lower bar of the model corresponding 
to full development of the "passive’ state of stress 
within the soil adjacent to the haunch. 


t 
5 by 
Then, from Figure 4.12: 


a 
Eg (Sa) (4.4.4) 
Ly Q 
8) 


Substituting Equations (4.4.2) and (4.4.3) into Equation (4.4.4): 


(ae - 
. | z~ (K, = Kop) (7,4, = 5P 55) cos(w - &) + (2)ein(w - 5)] | 
B= t,, 


%5 \( 5) 6) + (2) 6 
ba (Koa = Ky 7 545 = oP sj cos + aa = )sin(w + ) 


(4.4.5) 


is very small except for 


t —_ 
Ly 


Now, the difference between Py j and Ps, 
the brief interval during which the shock front traverses the structure and 
the adjacent soil wedges. This is due to the relatively small dimensions of 
the structure as compared with the spatial dimensions and slow decay with time 
of the nuclear blast wave. Further, the large differences which do exist 


between Ps and p! during the shock front transit, when pressure buildup 


J 
occurs, are not significant here because during this interval the structure 
does not have time to respond sufficiently to develop significant soil 


resistances. Thus, Equation (4.4.5) may be reduced to: 


ro So 7 aes! cos(w - 5) + (b/a)sin(w - 6 
* by 7 * ty Fax : K, | | $28 w+) + (b/a)sin(w + o (4.4.6) 


Now, note that: 


cos(w - ©) + (b/a) sin(w - 5); _ cos (w - 5) + tan w sin(w - 6 a 
son} + &) + (b/a) sin(w + 5)} ~~ cos : + 5 + tan w ae a ~ 


(4.4.7) 
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Therefore: 


Ka - Kop 


ee ice 


(4.4.8) 
ty “Kon ~ Ky 


' — 
Alternatively, the idealized soil resistance-displacement function, 
Figure 4.12, may be defined using either the maximum soil resistances or the 
"yield" displacements for the soil masses as derived above, and the slope of 
the “elastic line, or “soil stiffness." For an interior haunch the "soil 


stiffness" may be expressed as: 


& 
= Q, 
K .= ot (4.4.9) 
sj” Ey 
where K = "stiffness" of the soil during the development of the 


"active state of stress. 
Substituting Equation (4.4.2) into Equation (4.4.9): 
ae 
pot te (Koy = Ky)(7—8, - 37g3) | cos(w + 8) + B)sin(w + 8)| 


or: Ko = k' -k Pay (4.4.10) 


q@? 
where k' = oor - K,) | cost a eee (2)sin(w " ») 
(<4) k! 


(ee 


k" 


While for an exterior haunch, the'soil resistance’ may be expressed as: 


Q! 
= ng (4.4.11) 
Ly 
where K'. = "stiffness" of the soil during development of the "passive" 


SJ 


state of stress. 





64 


Substituting Equations (4.4.3), (4.4.7) and (4.4.8) into 


Equation (4.4.10) gives: 


ac 


Baca ee Se), ', - Faby) eae) ss = 2) 


8 J é cP ~ NOP ) cos(w - 5) + (b/ajsin(w - & 
ly Koa ~ K, ae a ~ + pease - f 


ors webs ap Pe a ee 
Rag = k k Pe 


; (4.4. 1m 


1 


Note that the development of the two expressions for "soil stiffness, ' 


K. and Kt (Equations (4.4.10) and (4.4.12)), is consistent with the assumption 


of a constant slope for the "elastic'’ line of the soil resistance-displacement 


function for a given value of average overpressure, i.e. for Pej = Ps ; 
Equations (4.4.10) and (4.4.12) express the soil stiffness in terms 


of a " 2. 5 os, " diagram (cf. Figure 4.12). It is sometimes more convenient 
to deal directly in terms of the horizontal displacements of the haunches 


rather than with the rotations of the lower bars. Now, from Equations (4.2.9): 


80, = 80, = - 6) = - h Ja. (4. 4. a 


Therefore, the soil stiffnesses may be expressed as: 


(=) K = = (k'-k" p.,) (4.4.24) 


sj sj 


and: (=) K' 


1 a | ee | 
te = (k' - k PL) (4.4.15) 


At this point one should consider the resistance offered by the 
soil as the structure "rebounds, or reverses its direction of motion, during 
response to a blast load. Keeping in mind the fact that the overpressure 
tends to force the soil into intimate contact with the arch, the variation in 


soil resistance during rebound is assumed to be as shown in Figure 4.13. 
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It is now possible to examine more closely the classification of 
haunches according to their displacement into "interior" and "exterior’ 
haunches as was done at the beginning of this chapter. Referring to Figure 
4,15, the original position of the haunch is represented by the origin of the 
coordinate system. The effect of the “permanent set” is to transfer the 
origin of the coordinate system to allow for any "plastic" deformation which 
may take place in the soil during the response. This is, of course, an 
_ approximation of the behavior of the soil, but it is consistent with the 
behavioral assumptions made previously in this chapter. Note that a reversal 
in the direction of haunch movement does not cause an immediate reversal in 
the direction of the soil resistance. There is some evidence to indicate ‘2° 15) 
that the reversal of the direction of the soil resistance will occur more 
rapidly with reversal of the direction of movement than Figure 4.13 indicates; 
however, further investigation of this behavior was not included in this 
study (cf. Chapter IX). 
Consistent with the sign conventions used elsewhere in this thesis, 


soil resistances acting to the left are positive, and conversely. 
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Figure 4.1. Forces Acting on Pressure Surface for "Interior" Haunch 
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Figure 4.3. Soil Resistance-Displacement Diagram for 
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Figure 4.4. Effective Soil Resistance for "Interior" Haunch 
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Figure 4.6. Forces Acting on Pressure Surface for "Exterior" Haunch 
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Figure 4.8. Effective Soil Resistance for "Exterior" Haunch 
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Figure 4.10. 
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Figure 4.11. Partial, Idealized, Soil Resistance-Displacement 
Diagram Showing Effect of Overpressure 
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Figure 4.12. Complete, Idealized, Soil Resistance-Displacement 
Diagram 
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Figure 4.15. Assumed Variation in Soil Resistance During Rebound 





1) 


CHAPTER V 


MASS EFFECTS 


The mass of a system subjected to dynamic loading directly affects 
the inertial resistance of the system and thus has an important effect on 
both the magnitude and the time history of the response of the system. It 
is the purpose of this chapter to list and discuss the assumptions made 
regarding the mass of the arch-soil “structure” and to derive the relationships 
necessary to incorporate inertial resistance into the equations of motion of 
the system (cf. Chapter VII). 

In general, the approach used here is to consider the masses of arch 
and soil which may logically be assumed to move as the system responds. The 
inertial resistances of these masses are then calculated and applied to the 
model as a system of loads. The resultant effect of these inertial forces 
is taken as the equivalent inertial resistance of the model. As in the case 
of the equivalent deformation resistances determined earlier, the equivalent 
inertial resistance is determined in agreement with the choice of independent 


displacements for the model. 


5-1 Mass Concentrations on the Model 

Since the model has three interior joints, masses may be placed at 
each joint in such a manner as to approximate in the most realistic possible 
fashion the inertial resistance of the actual system in the region adjacent 
to that joint. The model itself imposes no inherent restrictions on the choice 
of the masses which may be placed at the joints. The mass of the arch may be 
handled quite readily by concentrating at each joint the mass of the arch rib 
for a one-half bar length to either side of the joint. Note that this results 


in the concentration of the mass of the lower half of each of the lower bars 
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at the arch abutments. Since the abutments are assumed immovable (cf. Sec- 
tion 1.2), this portion of the mass of the arch is thus assumed to have no 
significant effect upon the response of the system. As pointed out in 
Section 2.3, such assumptions are not uncommon in dynamic analyses and give 
reasonable results. 


The mass of soil which may logically be assumed to move with the 





arch is more difficult to define. Im general, the mass of soil moving with 
a joint depends upon: 

1. The location of the joint, e.g. the mass of soil moving with 
the crown will be different from the mass of soil moving with one of the 
haunches. 

e. The direction of movement of the joint, e.g. the mass of soil 
affected by the inward movement of a haunch is different from that affected 
by its outward movement. 

5. The depth of burial of the arch, e.g. the mass of soil moving 
with the crown will tend to increase as the depth of cover increases. 

Further discussion of the latter two considerations is warranted 
here. With regard to the direction of movement of the joint, note that this 
approach is somewhat different from that used in Chapter IV in developing the 
soil resistance due to deformation. The effect of both approaches is the 
same during initial response. However, once rebound occurs, the "deformation’ 
resistance of the soil is assumed to decrease linearly to zero before building 
up in the opposite direction (cf. Section 4.4). On the other hand, by con- 
sidering only the direction of movement when choosing the effective soil mass, 
the start of rebound causes an immediate change in the magnitude of the mass 
of soil considered to move with the system. This apparent contradiction is 


not as serious as it might first appear, since the start of rebound corresponds 
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to a maximum displacement, hence zero velocity, of one of the two independent 
motions of the system, so that there is actually no momentum when the mass 
assumption is changed. The expressions derived below permit the mass assump- 
tions to be changed independently for each degree of freedom of the systen, 
since in general the independent motions will not reverse their direction of 
movement, i.e. “rebound,' at the same time. 

With regard to the depth of burial of the arch, it has been suggested 
in the case of rectangular, buried structures that a depth of cover greater 
than the span length should not be used in estimating the mass of soil moving 


tithe structure, (2? 


This is a logical restriction and is incorporated 
into the derivations in this chapter. 

In deriving the mass concentrations on the model consider first that 
these concentrated masses will affect the inertial forces developed at each 
of the joints as the system responds to a blast loading. Second, the magnitude 
and direction of these inertial forces is variable but can be determined at 
any time during the response. Now, instead of dealing with the resultant 
inertial forces at the joints, consider instead the horizontal and vertical 
components of these forces. The mass concentrations on the model which are 
derived in the remainder of this section correspond to these horizontal and 
vertical components. This approach is more convenient than dealing with the 
masses corresponding to the resultant inertial forces at the joints. Note 
that this concept is not based upon horizontal and vertical "components" of 
the masses themselves. 

The masses at each joint corresponding to the horizontal and vertical 


components of the inertial forces at the joint may be expressed as: 


no Wag + ae (5.1.99) 





where 


m 


ms. 
BJ 


ne. 
ad 
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mass concentration at joint j corresponding to the 
horizontal, or the vertical, component of the inertial 
force at joint j. 

soil contribution to im 


arch contribution to ee 


As mentioned previously, the mass contribution of the arch at each 


joint is the mass of the arch rib for half a bar length to either side of the 


joint. 


where 


This may be expressed as: 


tl 


ad 5 ahs 


S22 (L. | +L,) (5.1.2) 


mass contribution of the arch at joint j 


weight per unit volume of the material from which the 
arch is made 


thickness of the arch rib 
acceleration due to gravity 
length of the bar to the left of joint j 


length of the bar to the right of joint j. 


The masses of soil associated with the horizontal components of the 


inertial forces at the joints are shown in Figure 5.1. While both of the 


haunches are free to move either inward or outward, Figure 5.1 is presented, 


for convenience, as though the left-hand haunch, joint 1, were moving inward 


and the right-hand haunch, joint 3, moving outward. Note that in order for 


the dimensions of the soil blocks not to exceed the "corresponding span” of 


the model, restrictions have been placed on two of the dimensions as follows: 


day = (S+¢+4,) < (b+ a) (5.1.3) 
Cc 
do = (5 +4,) < 20 (5.1.4) 


where the dimensions are defined in Figure 5.1. 
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Referring to Figure 5.1, one may write: 








Y 
= 2 b Sa) Ore 
(moa in * Bg 4a) Ce) ee 
atc \;d 2 rn 
+ (Gey ear ne) vp oles cot(y + 3) | (5.1.5) 
where Cee Pe = mass contribution of the soil at joint 1 corresponding 


to the horizontal component of the inertial force acting 
at joint 1 when the joint is moving inward 


ae ¢ 
Mo = a une - i) (5.1.6) 


where m = mass contribution of the soil at joint 2 corresponding to 
the horizontal component of the inertial force acting at 
the crown for both inward and outward movement of the 
crown 


~ 


5 a, ;b+d 
(0s out i se | dy (BD x (dy 7 aCe 
Toe ate )(S) + ae, cot(T - 3) (5.1.7) 


where (m,.,) = mass contribution of the soil at joint 3 corresponding 
to the horizontal component of the inertial force at 


joint 3 when the joint is moving outward. 


out 


Note also that: 
(moan _ (m1 )in (5.1.8) 


where (moan = mass contribution of the soil at joint 35 corresponding 


to the horizontal component of the inertial force at 
joint 3 when the joint is moving inward 


(m5) (m,5) (5.1.9) 


out = out 


where (m,,) mass contribution of the soil at joint 1] corresponding 
to the horizontal component of the inertial force at 
joint 1 when the joint is moving outward. 


out . 


The masses Of soil associated with the vertical components of the 
inertial forces at the joints are shown in Figure 5.2. Note that the 


restrictions stated by Equations (5.1.3) and (5.1.4) also apply here. 


a 





Referring now 


where moh 
where m 
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where M6 

Now 
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to the figure: 


It 


un 
may = E (4a) + Cy - NPE) + Cy - SEB] (5.2.20) 


mass contribution of the soil at joint 1 corresponding to 
the vertical component of the inertial force at joint 1 
for both inward and outward movement of the joint 


yd 
ose ee - c/4) (5.1.11) 


mass contribution of the soil at joint 2 corresponding to 
the vertical component of the inertial force at joint 2 
for both inward and outward movement of the joint 


7 
so | dy (0/2) + (a, - a) (BES eae ae Ge 2) (5.1. diag 


mass contribution of the soil at joint 5 corresponding to 
the vertical component of the inertial force at joint 4 
for both inward and outward movement of the joint. 


substituting Equations (5.1.2) and (5.1.5) through (5.1.12) 


into Equation (5.1.1): 


: 
(men = BE | ayy B) + (gy - SPE + Cay - 8D) 


(m, ) 


(m 


(m, 


my 


2 1D, | 
+ qd 4 cot(y A 2) | + 7 nt/g 


Y 
out 2g ae (3) + (gq - BPE) + (a, - 832)6) 


din 


} 
o i ry | 
f 7 cot(y - 2) | - 7 ho L/g 


7 a 
yeas _¢ 
-  g (d_ 0 i) e 7 sh t/e 


(m,), 


out (my) out 


7 
_ Js faa) + (gy - SP) + (ay - $29G)| + 2,,1/e 


2g 
— (5.1/2 
7 2 


= my 
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Equations (5.1.13) are expressions for the effective mass con- 
centrations at the joints of the model. They are used in the next section in 
deriving expressions for the equivalent masses of the two-degree-of-freedom 


system. 


5.2 Development of Equivalent Mass 
It is now necessary to combine the inertial resistances of the 


arch-soil "structure" in such a manner as to obtain their resultant effect on 
the response of the model. Again, the expressions which follow are developed 
considering the horizontal displacement components of the haunches as the 
independent variables. 

The technique followed here consists of using the mass concentration 
expressions derived in the preceding section to write the horizontal and 
vertical components of the inertial forces acting on the model at each of the 
joints. From the fact that in the response of any one of the masses the con- 
ponent of its acceleration in a given direction is proportional to the correspond- 
ing displacement component, the accelerations are related to one another using 
the "Displacement Relationships’ derived for the model in Section 3.2.2. The 
inertial forces are then applied to the model and their resultant effect in 
the horizontal direction at the haunches found using the "Load Relationships” 
derived in Section 35.2.1. 

Thus, using Equations (3.2.7) and the principle stated above, one 


may relate the acceleration components of the masses, ae as: 


h = + = (3 - be/ad)(h, S h;) 


C 
hy = (b/a)h, 
ho = + = (a/e - b/a)(h, - h,) (5.2.1) 
hg = + (o/a)n, 
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where h, = horizontal component of the acceleration at joint 1 
h, = horizontal component of the acceleration at joint 2 
h, = horizontal component of the acceleration at joint 3 


hy, = vertical component of the acceleration at joint 1 
h. = vertical component of the acceleration at joint 2 

hg = vertical component of the acceleration at joint 35. 

Note that the acceleration components chosen for reference, i.e. hy and hs 

correspond to the horizontal displacement components of the haunches, which 

have already been chosen as the independent displacement coordinates for the 
dynamic analysis. 

Using Equations (5.1.13) and (5.2.1) it is now possible to write 
the horizontal and vertical components of the inertial forces acting on each 
joint of the model. Since there is nothing to be gained by rewriting here 
the mass concentration expressions given in Equations (5.1.13), the notation 
will be simplified by merely indicating the proper mass concentration, oe to 
use for each inertial force component. Thus, the inertial force components, 


oy? are: 


= mh, 


=> 


z= (m/2)(1 ~ ve/ad)(ny + by) 

Cains ; (5.2.2) 
2 = - m) (b/a)b, 

Z. = (m,/2)(a/e : b/a)(b, : h,) 


= mg (b/a)h, 


of 
! 


By applying the inertial force components, a to the model as a 
system of loads, the equivalent inertial force in the horizontal direction 


at the haunches, 2, may be found, using Ey uations (3.2.5), ess 
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Z= +b - By, - Garde - by, 
(re. 57 
2, = = (1 - 7E Jz, + 2. - = Ge *)2, - (2 7% 
Now let Ds on = equivalent mass of the system in the horizontal 
direction at haunch j. Then: 
Z, - Ms bh, (5.2.9 
Substituting Equations (5.2.2) and (5.2.4) into Equations 
2.5): 
By Py = m,h, + (m,/4)(1 - ve/ad)*(h, + h; ) 
+ m,(b/a)“h, 4 (m,/4) (a/c - b/a)*(h, - h,) 
my hs = (my/4)(2 - ve/ad)*(in, + h,) + mh, 
e ae aD oc 
- (m,/4) (a/c - b/a) (h, - h,) a mg (b/a)“h, 
Rearranging terms gives: 
in = Ee : (m,/4 )(2 - ne eal) + m (b/a)° + (m,/4)(E - S| 
h 
-3{- -(m,f4)(2 + ve/ea)® + (mg/4)(a/e ~ v/a)? | 
al 
. (5.2.5) 
By ay * | (mg/¥)(2 ~ be/ad)® +m, + (mg/4)(S ~ 2)? + mg (o/a)?| 
2 [aya - ve/ad)* + (m,/4)(a/c - v/a)? | 
Note that either (man or (m5), ut may be used for a particular joint in 


Equations (5.2.5) as the motion of the system dictates. 
Since the mass effects at the haunches are not separated in 
Equations (5.2.5), it is not possible to treat them separately here as was 


done in the derivation of expressions for equivalent soil resistance in 
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Section 4.4. This, together with the fact that each of Equations (5.2.5) 
can have either of two possible values, one for inward and one for outward 
movement of the corresponding haunch, makes it necessary at this point to 
consider the possible motions of the model (cf. Section 1.2). There are four 
of these motions as illustrated in Figure 5.3. For convenience, they are 
designated as follows: 

#1 = Asymmetric-Right. This motion is characterized by the 
simultaneous movement of both haunches to the right. 

#2 = Symmetric-In. This motion is characterized by the simultaneous 
inward movement of both haunches. 

#3 = Symmetric-Out. This motion is characterized by the simultaneous 
outward movement of both haunches. 

#4 = Asymmetric-Left. This motion is characterized by the simul- 


taneous movement of both haunches to the left. 


Note that in the case of each of the four motions defined above the crown may 
move either upward or downward depending upon the relative velocities and dis- 
placements of the haunches. 


The notation of Equations (5.2.5) may now be simplified as: 


7 
m, = ey 
1m” “I,m b, l,m 
7 (5.2569 
h 
= 1 
"3 ym = "3 ,m E 7 "3 ym 
vhere Hot mt (m,/%)(1-be/ aad )* + m,(b/a)> + (m,/4)(a/e-b/a)° 
My im = 7 Kitg/4)(1-be/aa)® + (m,/4)(a/e - b/a)® 
Hs = (mz/4)(1-ve/aa)” + m, + (mz/4)(a/e-v/a)° + me(b/a)® 
ae = = (m,/')(1-be/aa)° + (m,/%)(a/e-b/a)* 


m = index denoting the motion of the system as defined above. 
Meme. 54: 
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1. Asymmetric-Right 2. Symmetric-In 
eae SS 
_ 
3. Symmetric-Out 4. Asymmetric-Left 


Figure 5.3. Motions of Model 





CHAPTER VI 


BLAST PRESSURE LOADING 


It is the purpose of this chapter to define the loading on the 
dynamic model in terms of the geometry of the arch-soil "structure’ and the 
shock wave (blast overpressure) resulting from the explosion of a nuclear 
weapon. The majority of the blast parameters are taken from references (16) 
and (17). The shape of the overpressure decay curve behind the shock front 
is taken from reference (20). ‘The assumptions made in selecting these 
parameters are discussed in Section 2.4. 

In defining the load on a buried structure the shape and char- 
acteristics of the blast wave in soil are of obvious importance. Studies of 
the forces transmitted through the soil cover are being made by several groups 
at the present time. A summary of much of the significant work which had been 
done on this problem, together with a convenient numerical technique for solv- 
ing the one-dimensional problem of stress wave transmission through soil, is 
contained in reference (18). From these studies it appears that there is 
little or no significant attenuation in the vertical, free-field overpressures 
within the relatively shallow depths considered in this thesis. The results 
of a study of the configuration of the shock front in soil are presented 
herein as Appendix B. 

For derivations in this chapter the blast loading is considered to 
be transmitted by the soil directly to the arch without modification. This 
is consistent with the approach used in Chapters IV and V in deriving expres- 
sions for the deformation and inertial resistances of the soil surrounding the 
arch. The assumption is made that the ground surface is flat in the vicinity 
of the arch so that dynamic pressures and reflection phenomena due to surface 


irregularities need not be considered. 
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While the model and the methods of dynamic analysis used in this 
thesis are suitable for studying blast overpressure loadings due to shock 
waves approaching the structure from any direction, the expressions developed 
in this chapter are derived on the basis of a shock wave approaching from the 
side of the arch, i.e. with the direction of travel of the shock wave normal 
to the longitudinal axis of the arch. This geometry produces maximum asymmetry 
in the arch loading and thus should lead to maximum lateral response of the 
arch. 

The assumption that nuclear blast waves may be approximated using 
the idealized curves of references (16), (17) and (20) is basic to all of the 
derivations in this chapter. A convenient, but as will be seen not altogether 
essential, assumption is made that the seismic velocity of the soil is practical- 
ly constant between ground zero and the point of interest, i.e. the location of 


the structure. 


6.1 Definition of Loading Regions 


The arch itself may receive load at all points on its periphery 
during passage of the shock wave. However, in defining the load on the replace- 
ment model used for the dynamic analysis load is brought in only at the joints. 
These joint loadings are simply concentrations of the distributed loading 
sensed by the arch itself. The load concentrations are obtained by defining 
regions in the vicinity of each joint which will tend to bring load to that 
joint. Then the average blast overpressures acting on these regions are con- 
centrated and applied to the model at the joints. It is more convenient here 
to deal with the horizontal and vertical components of the load at each joint 
than with the resultant forces themselves, just as was the case with the 


inertial forces considered in Chapter V. 





In the following sections regions tending to bring load to each 
HOint, i.e. "loading regions,’ are defined for various shock wave-soil 
combinations. In general, for shallow-buried structures such as are considered 
in this study, the configuration of the shock front in the soil may be con- 
sidered as a plane perpendicular to the direction of travel of the shock front. 
It is oriented, in the plane of the arch, at an angle with respect to the 
shock front in air which is dependent upon the relative values of the seismic 
velocity of the soil and the velocity of the shock wave in air at the point 
of interest (cf. Appendix B). Since in this study the abutments are assumed 
immovable (cf. Section te loading regions tending to bring load primarily 
to the abutments are neglected. 
The three cases considered below are for: 
(1) A vertical shock front in the soil, 
(2) A shock front in the soil which is inclined at an angle 
less than the inclination of the upper bars of the model, and 
(4) A shock front in the soil which is inclined at an angle 
greater than the inclination of the upper bars and less than the inclination 
of the lower bars of the model. 
A series of problems (cf. Chapter VIII and Appendices C and D), 
which was solved using various angles of inclination for the shock front and 
holding the remaining variables constant, indicated that the lateral response 
of the model was not particularly sensitive to variations in this angle. Thus, 
the assumption of a constant seismic velocity in the soil between ground zero 
and the location of the structure, while convenient in defining the inclination 
of the shock front in the soil, is not too significant here. Therefore, the 
case of a shock front in soil inclined at an angle between the vertical and 


the inclination of the lower bars of the model is omitted. Also omitted is 
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the case of a horizontal shock front in soil. While this may be readily 
handled by the techniques developed herein, it is outside the scope of this 


study (cf. Section 2.4). 


6.1.1 Case I - Vertical Shock Front in Soil 

From the results of the study presented in Appendix B it can be seen 
that a nearly vertical shock front in the soil results when the velocity of 
the shock front in air at the point of interest and the seismic velocity of 
the soil are nearly equal. This case is illustrated in Figure 6.1. 

Since the configuration of the shock front in soil has been assumed 
as a plane perpendicular to the plane of the arch, it will, in general, strike 
the arch at some point of tangency or at an abutment. The problem of the 
shock wave enveloping the arch itself can be handled quite readily (cf. refer- 
ence (21)). However, for the purposes of this study it is more convenient and 
not significantly less accurate to consider instead the shock wave as enveloping 
the replacement model. 

In interpreting Figure 6.1 consider the shock front in air as moving 
along the ground surface from left to right while the shock front in soil, as 
an extension of the shock front in air, moves through the soil. Prior to the 
arrival of the shock front at point Oa', the shock wave is assumed to have no 
influence on the structure. As the shock front moves from 0a' to Ob, the 
shock wave engulfs the "passive" soil failure wedge corresponding to outward 
displacement of the windward haunch, joint 1 (cf. Section 4.2). The average 
overpressure acting at the ground surface between points Oa' and Ob is thus 
Ba As the shock front moves from point Oa to point Ob, the shock wave engulfs 
the "active" soil failure wedge corresponding to inward displacement of the 
windward haunch, joint 1 (cf. Section 4.1). The average overpressure acting 


between points Oa and Ob is thus Dae 
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As the shock front moves from point O to point la, the shock wave 
is assumed to load the left abutment, joint 0; therefore, this region is not 
considered in the arch loading. The shock front then moves from the midpoint 
of bar 1, point la, to the midpoint of bar 2, point 1b, and in so doing is 
assumed to load joint 1. The average overpressure acting between points la 
and 1b is designated P)- 

As the shock front moves on from point lb, the shock wave starts to 
load joint 2; therefore, for convenience, point lb is also labeled point 2a. 
Now the shock wave may be assumed to load joint 2 as the shock front moves 
from point 2a to point 2b. The average overpressure acting between points 2a 
and 2b is designated Pp: 

It may be noted here that the horizontal component of the loading 
between points 2a and 2 tends to act to the right at joint 2, while the hori- 
zontal component of the loading between points 2 and 2b tends to act to the 
left. For this reason it is convenient to designate the average overpressure 
acting between points 2a and 2 as ae and that acting between points 2 and 2b 
as Dp. 

As the shock front moves on from point 2b, it starts to load joint 3. 
Accordingly, it is convenient to label point cb also as point Ja. Now, the 
average overpressure acting between points 3a and 3b is assumed to load joint 3 
and is designated Ps. 

The average overpressure acting between points 3b and 4 is assumed 
to load the right abutment, joint 4, and is therefore neglected. 

As the shock front moves from point 4a to point 4b, the shock wave 
engulfs the "active" soil failure wedge corresponding to the inward displace- 
ment of the leeward haunch, joint 3. The average overpressure acting between 


points 4a and 4b is thus Pe: 
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As the shock front moves from point 4a to point 4b', the shock wave 
engulfs the "passive" soil failure wedge corresponding to the outward dis- 
placement of the leeward haunch, joint 5. The average overpressure acting 
between points 4a and 4b’ is thus Pass: 

For purposes of calculating the "rib-shortening" correction derived 
in Section 5.5, the average overpressure acting on the entire arch is needed. 
This is taken as the average overpressure acting between points O and 4 and 


Lf " 


is designated Pe For convenience in this case point 0 is also labeled aa 


"ab tt 


and point 4 is also labeled 
Note in the foregoing description that no actual “loading” of the 
arch takes place until the shock front reaches point la. It is assumed that 
no movement of the arch occurs prior to this time. Therefore, the time of 
arrival of the shock front at point la is taken as zero time for this case. 
This assumption is consistent with the derivation of the "deformation" resistence 
of the soil in Chapter IV. It is not unreasonable when one considers the high 
velocity of the shock wave, the relatively small dimensions of the arch-soil 
"structure, the inertia of the soil loaded prior to this time and the 
relatively small displacements of the structure during the early part of the 
loading history. 
In the derivation of expressions for the load concentrations on 
the model (cf. Section 6.2) it is convenient to use the times at which the 
shock front arrives at the windward and leeward boundaries of the various 
loading regions. Zero time for this case has already been chosen at point la. 
This is expressed algebraically by Equation (6.1.1). 


t,, = 0 (6.1.99) 


The remaining time parameters are given by Equations (G2); expressed in 
terms of the velocity of the shock front in air, U, and the geometry of the 


arch-soil "structure" as shown in Figure 6.1. 
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16 = - (b + 2d cot ae P ly/ou 
Oa O } ~ 9 
= = - (b + 2a cot E + Spe 
Oa re) y 2 
bee (b - aba /a)/2u 
ty = (> + a)/o2u 
aes 7 ays 
t, = (b + 2a)/2U 
(6.1.2) 
to = (> + 3a)/2uU 
a a op 
tap = (b+ 2a)/U 
ty fob 426 2bd_/a)/2U 
= (+ Fl EP yoy 
tis = -b + 2B + 2 oct as 2 | a 
ake: rl ee Ven 
thas = -b + 2B+ 2 cot ra 5 | 2 
to. 27 b/2U 
Cs (-b + 2B)/2u 


time of arrival of the shock front at the windward 
boundary of loading region j 


where, in general: ae 


time of arrival of the shock front at the leeward 
boundary of loading region Jj. 


as 


6.1.2 Case II - Shock Front in Soil Inclined at Angle Less Than 
Inclination of Upper Bars of Model 


The results of the study presented in Appendix B indicate that for 
seismic velocities less than the velocity of the shock wave in air at the 
point of interest, the shock front in soil will be inclined at an angle to 


the shock wave in air. The expressions derived in this section are for the 








2 


case where the inclination of the shock wave in soil with respect to the 
horizontal, 1, is less than the angle of inclination of the upper bars of 


the model, i.e.: 
1) < arctan 5 (6.1.3) 


This case is illustrated in Figure 6.2. Note that in this case 
the shock front in soil first strikes the model at the crown. Therefore, for 
this case the time of arrival of the shock wave at the crown is taken as zero 
time. 


t = 0 (6.1.4) 


The considerations in selecting this time as the reference are the same as 
for the preceding case (cf. Section 6.1.1). 

The loading regions shown in Figure 6.2 are defined in an analogous 
manner to those for the vertical shock front in soil and therefore are not 
described in detail here. However, there are several differences worthy of 
note, viz: 

1. The "active" and "passive’ soil failure wedges are not fully 
enveloped by the shock wave until the shock front in soil reaches the corre- 
sponding abutment. 

2. On the windward side of the model the shock front strikes the 
midpoint of the upper bar before reaching the midpoint of the lower bar. 
However, the midpoints of the bars are not relabeled. This introduces an 
inconsistency in the notation, but not a serious one, since the order in which 
the boundaries of a "loading region” are reached by the shock front is obvious 
in any set of computations. 

3. The crown, joint 2, is enveloped by the shock wave from the top 
downward instead of from left to right as was the case in the preceding section. 


Here also the notation is not changed. 
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4, The average overpressure acting on the arch for this case is 


taken as the average overpressure acting between points 2 and 4. This 


corresponds to the envelopement of the model by the shock wave as the shock 


front travels from the crown to the leeward abutment. For convenience, 


points 2 and 4 are also labeled "aa" and "ab." 


The time parameters defining the times of arrival of the shock 


front at the boundaries of the loading regions are given in Equations (6.1.5). 


These may be written readily by referring to Figure 6.2. 
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Equations (6.1.5) are written in terms of the apparent vertical 
component of the velocity of the shock front in soil, Ci. If the inclination 
of the shock front in soil, , with respect to the horizontal is obtained from 


Appendix B, Equation (B.3), C! may be obtained from Equation (6.1.6). 


Ci = U tan 9 (6.1.6) 


6.1.3 Case III - Shock Front in Soil Inclined at Angle Intermediate 
Between Inclinations of Upper and Lower Bars of Model 


The expressions derived in this section are for the case where the 
inclination of the shock front in the soil is greater than the inclination 
of the upper bars of the model and less than the inclination of the lower 
bars of the model, i.e.: 


£ 


es arctan = (6.1.79) 


arctan = 
This is illustrated in Figure 6.3. Note that the shock front strikes the 
windward haunch first. Therefore, the time of arrival of the shock front at 
the windward haunch is taken as zero time for this case. 


eitaatyn() (6.1.8) 


The considerations for the selection of ty as the reference time are the 
same as those listed in Section 6.1.1 for the selection of the reference time 
for the vertical shock front in soil. 

The loading regions shown in Figure 6.3 are defined in the same 
fashion as those in the preceding sections. However, several differences do 
exist for this case, viz: 

1. The windward haunch is assumed to be loaded as the shock front 
moves from the haunch to the midpoint of either the upper or lower adjoining 
bar depending upon which occurs later. In other words, the shock front may 


arrive at point la either before or after it strikes point lb. The order of 


arrival depends upon the relative geometries of the shock front and the model. 
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2. The soil failure wedges on the leeward side of the model are 
assumed to be loaded as the shock front moves from point 4a to either the 
leeward abutment or the intersection of the failure plane and the ground 
surface, whichever occurs later. The order of arrival depends upon the 
relative geometries of the shock front and the model just as in the case of 
the windward haunch discussed in the preceding paragraph. 

5. The average overpressure acting on the arch in this case is taken 
as the average overpressure acting between points 1 and 4. For convenience 
these points are also labeled for this case as “aa” and "ab". This corresponds 
to tne envelopment of the model by the shock wave as the shock front travels 
from the windward haunch to the leeward abutment. 

The time parameters defining the times of arrival of the shock front 
at the boundaries of the loading regions are given by Equations (6.1.9). These 


may be written readily by referring to Figure 6.3. 


to,1 = - (c + 4,)/c! - (b + a cot E : 3) iv 

cas = (c + a/c! - (bd + d cot E Zi 2 vo 

to, = a/C! - b/U 
t= 0/2! - b/ 2U 
22 c/2C! + a/2U a 
ea Ib 

Sy m= c/C) + d/l 

cae c/ 2c! + 3a/2U 

3, = Yop 

ts, = e/ec' + (ld + b)/2U 


* Use either t. or t.,., whichever is larger. 
la Lb 








the = - (c + d./c: + (-b +B- bd /a)/U 


+ 


me), = - (ec d/c:  (-b + B+ d cot E + 3 /\/v 


(ree = afC! + (-b + B)/U 


"(ty +), - (ce + a/c! + (-b+ B+ d cot E - 3 )i/v 
(typ (ti) 


an js 7 


cen (tino 
Equations (6.1.9) are written in terms of the apparent vertical 
component of the velocity of the shock front in soil, Ci. If the inclination 


of the shock front in soil, , with respect to the horizontal is obtained from 


Appendix B, Equation (B.4), then C’ may be obtained from Equation (6.1.6). 


6.2 Load Concentrations on Model 

The load concentrations cn the model are derived in this section by 
considering the average overpressure acting at the ground surface on the load- 
ing regions defined in Section 6.1. The approach used is illustrated in 
Figure 6.4. The overpressure is assumed to have no effect upon the region 
until time ae the time of arrival of the shock front at the windward boundary, 
ja, of the region, j. The average overpressure acting on the region, P,» is 
then assumed to build up linearly as the shock front traverses the region from 
the windward boundary, ja, to the leeward boundary, jb. At times greater than 


the time of arrival of the shock front at the leeward boundary, eae the 


* 
Use either (ty 1 oF (ty dos whichever is larger. 


iis, : 
Use either Opp or (ty dos whichever is larger. 
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average overpressure acting on the region is assumed to have a value equal to 
the instantaneous value of the overpressure at the center of the region. This 
approach is essentially the same as is used in reference (16) in defining the 
loading of flat surfaces oriented parallel to the direction of travel of the 
shock wave. Various techniques involving the integration of the overpressure- 
time pulse followed by an averaging of the impulse over the region of interest 
were studied and found to offer little or no advantage over the simpler approach 
described above. 

In general, the average overpressure acting on a loading region is 


expressed by Equations (6.2.1) 


fat p.(t) = 0 
uae P ,( ) 
t,, -t iG, ci. 
ja - jb Jj 2 api da 
ab = ‘iat’ 
ty Sts5 (t, +t p )+tgi p ,(t)= DA ae ) (6.2.9 
a) — 
z (t. tt» )ttg < t p ,(t)= O 
where p(t) = average overpressure acting on region j at any time, t 


(cf. Figure 6.4) 


p(t) = variation of overpressure behind the shock front with 
time (cf. reference (20)) 


5 a = time 
19 = time of arrival of the shock front at the windward 
Ja 
boundary of region j 
vay = time of arrival of the shock front at the leeward 


boundary of region, j. 

It has been pointed out earlier, in the cases of the inclined shock fronts, 
that because the reference points on the model are labeled from left to right, 
inconsistencies in the notation exist when the shock front approaches a region 


from a different direction. This occurs, for example, in Case II where the 
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shock front strikes point lb before arriving at point la (cf. Figure 6.2). 
This type of situation may be handled readily using Equations (6.2.1) by con- 


sidering the earlier time as t. and the later time as t irrespective of 


Ja Jb 
the designation of the corresponding points on the model. 

Expressions for the concentrated loads acting at the joints of the 
model may now be written in terms of the average overpressures acting on the 
corresponding loading regions, Equations (6.2.1). Here the assumption is 
used that the overpressure acting at the ground surface is transmitted vertical- 


(18) 


ly downward by the soil as a vertical overpressure in the soil. Further, 


it is assumed that the horizontal overpressure in the soil, BP bears a con- 


stant relationship, Kor? to the vertical overpressure, p,015+17,19) 


ae = P,/Py (6.2.2) 


Note that the following expressions, Equations (6.2.3), are for the horizontal 
and vertical components of the load concentrations at the joints rather than 
for the resultant loads themselves. This is simply for convenience in deriva- 
tion and for use later in determining the equivalent loads on the system 


(cf. Section 6.5). Referring to Figures 3.2 and 6.1, 6.2 or 6.4, one may 


write: 
at — 
H) = - 5 Kp (atc) 
a — — 
us rus 3 Ky (Pog Pap) 
on ib re 
HL, te K Pz (atc) (6.2.3) 
H), = + = p, (dtd) 
i. = + Pod 
= pe patora) 
6 anes 


Note that since the overpressure acts downward, it is considered negative in 


accordance with the sign convention used herein (cf. Section 1.2). 
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6.3 Development of Equivalent Load 


The concentrated loads acting at the joints of the model must be 
combined in order to obtain their resultant effect upon the two-degree-of- 
freedom system. The “Load Relationships’ derived in Section 3.2.1, 
Equations (3.2.5), express the resultant loads in terms of the equivalent 
horizontal loads at the haunches. This is in agreement with the choice of 
the horizontal displacement components of the haunches as the independent 
variables for the dynamic analysis. 

While Equations (6.2.3) and (4.2.5) may be used successively to 
Obtain the equivalent load on the system, it is more convenient to combine 
them and express the equivalent loads in terms of the average overpressures 
acting on the various loading regions. Equations (6.3.1) are obtained by sub- 


stituting Equations (6.2.3) into Equations (3.2.5). 


Mlk 


H=- 5 Py | Kor (ate) (o/a)(o+a) | - 5; (Po, Pop eK, (1-be/ad) 


ppa(a/e - b/a) 


NM |e 


(6.3.1) 


(Po_"Pop )eK,,(1 - be/ad) 


wr! 
mh 


- 2 piala/e v/a) + £55)[K, (are) + (o/a)(b+a) | 


NM 


The notation in the preceding equations may be simplified as shown 


in Equations (6.3.2). 


Hy = - Jp, - Jp(Py, - Pop) + JsPp 
a ms a a o (6.3.2) 
ney - J,(Py, - Pop) ~ SS ADE 
where , = 5K, (atc) + = (b/a)(d+a) 
1 
J, = 5 cK, (1-be/ad) 
J, = = a(a/c-b/a) 





ote that the dimensions of Jy» J, and J, must be in agreement with the 


BH , 


nits in which the overpressure is expressed. 
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Loading Regions for Case II 


Figure 6.2. 
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Figure 6.4. Generalized Overpressure-Time Diagram 
for Typical Loading Region 
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CHAPTER VII 


EQUATIONS OF MOTION 


Now that expressions for the equivalent forces, masses and re- 
sistances of the two-degree-of-freedom system have been derived, it is possible 
to develop the equations of motion for the system in terms of these quantities. 
While a simplified notation is used throughout this chapter, the expressions 
may be expanded readily by inserting the proper expressions from Chapters III 
through VI. As an aid in visualizing the two-degree-of-freedom system, the 
equivalent masses, resistances and loads are considered to exist and act at 
the haunches, i.e. at joints 1 and 3 of the model. This is only a physical 
analogy, but it is convenient in comparing the expressions derived in this 


chapter with those of the preceding chapters. 


7-1 Development of Equations 


The forces acting on one of the two equivalent masses are shown in 
Figure 7.1 in general terms. A general expression for the equations of motion 
may be obtained by applying d'Alembert's Principle and equating the sum of the 
horizontal forces acting on the mass to zero. The result is given in 


Equation (7.1.1). 


ms + Q, S Qe 5 = H (7.1.1) 
where My = equivalent mass of j 
h, = acceleration of mass ; 
Q, = equivalent deformation resistance of the arch at j} 
Qs 3 = equivalent deformation resistance of the soil at j 
H, = equivalent load acting on mass j. 


Since each of the equivalent masses is actually a function of both 


of the independent accelerations (cf. Equations (5.2.6)), the two equations 
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of motion as expressed by Equations (7.1.1) are not in a convenient form for 
integration. However, since Equations (5.2.6) are linear in the acceleration 
terms, the equations of motion may be rewritten explicitly for the independent 
acceleration components as shown in the following derivation. 

Equations (7.1.2) are obtained by substituting Equations (5.2.6) 


into Equations (7.1.1). 


eee aes eS ye etal 
Boro eee 5. 5.” S85 


Solving Equations (7.1.2) simultaneously for h, and h, gives Equations (7.1.3). 


(7.1525 








1 3 

, H m = oe a m 
a ee ieee =o) + (Heedmeci=) 

ug ee a See 1 % sl My m3 m7 "1,m"3,m as Ss 83 

(7.135 
a ee, ee lym (i, -0,-3,) 
h, = H. -Q, -Q + SOs 
i Q sl Bey an i 9 Rs, 83 


7m’ 3,m>1,m°5,m 
The notation of Equations (7.1.3) may be simplified as shown in Equations 
et). 
hy 7 Fy (Hy -% -9, 1) ui Gy (Hs a,-0.5) 
h, = Ps (Hy, -& -9,y) + Gs (Hs ~ Q,-Q.5) 


Lt 


where F ee 


(7208 


ae 1,m3,m 1,m°3,m 
G = de 
1,m 


lu “Vv v 
l,m’ 3,m l,m 5,m 


V 


ss 


oa “l,m 3,m 1,m3,m 
cue am 
2m "1,m'3,m 1,m'5,m 


m = index denoting the motion of the system as defined in 
Section 5.2. “m= 1,2)3,4. 


er 
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Equations (7.1.4) may now be integrated sinultaneously in order to 
obtain the velocities and displacements of the system at any time under a 


given set of dynamic loads. 


7.2 Natural Periods of Vibration 

The natural periods of vibration of any system are highly signifi- 
cant parameters in predicting the response of that system to dynamic loads. 
The natural periods of the model which is used in this study for the dynamic 
analysis may be found quite readily from the equations of motion derived in 
the preceding section. 

Before proceeding with the derivation of expressions for the natural 
periods of vibration of the model, however, it is convenient to express the 
resistance of the arch, Equations (3.4.7), in a slightly different forn. 
Equations (3.4.7) express the equivalent resistance of the arch in terms of 
the resisting moments developed at the hinges of the model. While this is con- 
venient for use in an elasto-plastic analysis, it is inconvenient for the 
purpose of finding the natural periods of vibration corresponding to the 
initial, elastic response of the system. 

The elastic resistance of the model may be expressed in terms of 
the independent displacement coordinates by substituting Equations (3.2.11) 
and (3.3.1) into Equations (3.4.7). 


Q, = K,(A,by -Ash, )(a,ta,hy -a,hs) - Ky(-Aghy +Aghs )(a+a gh) -a,bs ) 
(7.27 


SE K,(-A,h) +Aph, (a, +a,h) -aph; } + KA(-Ajb) +Agh, )(q,+a,hy ~aghs } 
Equations (7.2.1) may be simplified by limiting consideration to small dis- 
placements such that the products of the displacements are small with respect 
to the remaining terms and may be neglected. The simplified expressions are 


given as Equations (7.2.2). 
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Q = (Aya, +459, Kb - (A, +839, kph. 
R, aan (A,a+A59, Kh) - (A)9, +459, )Kgh,, 


The notation in Eyuations (7.2.2) may be simplified as shown in Equations 


(7:24 


7.2.3). 


(7.2.38 
85, = - T,h, + th 


vhere t, = (Aja, +4,0,)K, 
(K,/2L°)(5+a/c+b/d+2c/a+be/ad ) 
5 7 aerials: 
(K,/2L")(3 + a/c - v/a - be/ad) 


~ 
| 


The discussion will now proceed with the derivation of expressions 
for the natural periods of vibration of the system. One may write the equations 
of motion for the system in free vibration, Equations (7.259 directly from 
the general equations of motion for the system, Equations (7.1.1), by elim- 
inating the loading terms. 


rite) @ a eae 0 (7.2.4) 


If a sinusoidally varying time function is assumed for the motions of the 
system, the accelerations, hy, and displacements, hy) of the system are related 


as shown in Equations (7.2.5). 


ee e 

h =e. W h woe 

yor eh, (7.2.5) 
where W = circular natural frequency of the system for the particular 


motion being studied (cf. Section 5.2). 
Equations Gage) may be rewritten in expanded form by inserting Equations 


(5.2.6) and (7.2.5). The result is given as Equations (7.2.6). 
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2 2 ~ 
ia tem mm 37 et 8g 
, (7.2.6) 


Q a 
nan hy . HS mm h, rs &, + 863 =A) 


Substituting Equations (4.4.14) and (7.2.3) into Equations (7.2.6) gives: 


2 2 = 
4) mn PL - “ayn Bs - Thy + T3hs - (K,/a)by = 0 
(7.2.7) 
-V we h. +H we h, + Tih, - Tih, - (K /a)h = Q 
mmm 1 Dm ms > i maak io s4 is 


Note that Equation (4.4.15) must be used in lieu of Equation (4.4.14) for 
either of Equations (7.2.7) when the motion being considered is such that the 
haunch is displaced outward against the soil rather than inward into the arch. 
This qualification also applies to the remainder of this derivation. Rewriting 
and rearranging Equations (7.2.7) gives: 

(4) - t) - K,,/a)hy + (-¥) wt + 1; )h, 


(7.2.8) 
2 2 5 
(mv, Yh + tay rn (u, wo = 1 K.2/a)h, = 


ym m Ee 


il 
© 


| 
© 


Equations (7.2.8) yield solutions different from zero only if their 
determinant vanishes. This permits Equation (7.2.9), the frequency equation, 


to be written in terms of the circular natural frequencies of the system. 


h pe 
| 1 ,m"5,m : “am”, “m ; nym , K,3/8) 
a 1 2 
+ Hy Ty + K_,/a) - an + ere a (7.2.9) 


+ | (+ ty Ky/aory - Ky/a) - (5)*) = 0 


The circular natural frequency for any one of the four motions of 
the system (cf. Section 5.2) may be obtained by substituting the proper values 
of Ba a? Yam? and Ce (cf. Equation (5.2.6)) and the proper expressions 


for soil stiffness, i.e. either K. or Ki (cf. Equations (4.4.14) and (4.4.15)), 


J 
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into Equation (7.2.9) and solving for wes Note that since both Ky and Ke 
are functions of the overpressure restraining the soil, the circular natural 
frequency, hence the natural period of vibration, of the system will vary as 
the shock wave passes over the structure. Im general, solution of Equation 
(7.2.9) yields two values of we, The proper value may be selected by sub- 
stituting the values, in turn, into Equations (7.2.8) and noting the relation- 
ship between the displacements hy and n° The proper value is indicated when 


the relationship between h, and hz agrees with the system motion being 


i 
studied. 

Once the proper circular natural frequency has been found, it may 
be substituted into Equation (7.2.10) and the natural period of vibration of 


the systen, TW obtained for the particular system motion and overpressure 


conditions being studied. 


On 
Boo (7.2. 108 
m 


While there is no real advantage in solving Equation (7.2.9) in 
general terms, a simplified expression for the circular natural frequency can 
be written for those special cases in which the stiffness and mass parameters 
are symmetrical. In this study this condition exists, for example, for the 
“Symmetric-In" and "Symmetric-Out" motions with no overpressure restraining 
movement of the soil. Equation (7.2.11) is obtained by setting La es. 


H eta a on and K., = K 5 = K, in Equation (7.2.9) and solving for the 


value of _ corresponding to symmetrical response. 


T. To —eene 
Ar as Wil (7.2.10 


m J _ te oS 
Note that K'/a may be used in place of K /a in Equation (7.2.11) if the motion 


being studied so indicates. 





Figure 7.1. Forces Acting on Typical Equivalent Mass 
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CHAPTER VIII 


NUMERICAL STUDIES 


This chapter covers the basis of the numerical studies made in 
connection with this thesis. The range of variables included in the studies 
is given together with a description of the blast wave properties used. The 
numerical examples are described and, finally, a set of sample calculations 
is included. The results of these studies are presented in tabular and in 


graphical form in Appendices C and D. 


8.1 Range of Variables Studied 
The quantitative investigations made in connection with this thesis 


are designed to illustrate the behavior of shallow-buried, two-hinged arches 
subjected to nuclear blast loading. Two different arches are used, one a 
semicircular, reinforced concrete arch, 16'S" in diameter, and the other a 
corrugated metal arch in the form of a segment of a circle with a rise of 
11' 8 1/2" and a span of 25' 11/4". The concrete arch has a uniform thick- 
ness of 8" and is reinforced with No. 4 bars running circumferentially on 
10" centers at mid-thickness. Intermediate grade reinforcing steel is used 
with an assumed yield strength of 50 ksi under dynamic load. The elastic 
modulus of the steel is taken as 50,000 ksi. The average strength of the 
concrete as measured by the standard 6" x 12" test cylinder is taken as 

3.5 ksi at 28 days. Substituting this value for the concrete cylinder 
strength into the appropriate expressions of reference (22) yields a corre- 
sponding elastic modulus for the concrete of 3,820 ksi. The corrugated 
metal arch is 10 gauge with an elastic modulus of 30,000 ksi and an assumed 


yield strength of 30 ksi under dynamic loading. 





Bo iy, 


The arches are studied at depths of cover over the crown varying 
from zero to one-half of the arch span. Intermediate depths of cover of 
one-eighth and one-quarter of the span are used. 

The soil cover is taken as a granular, cohesionless material with 
a unit weight of 120 pcf, an angle of internal friction of 30° and a hydro- 
static pressure ratio of horizontal to vertical soil pressures in the "at 
rest’ state of stress of 0.5. The friction angle between the soil and the 
arch is assumed as two-thirds of the angle of internal friction of the soil, 
or 20°. The seismic velocity of the soil is assumed as 3,000 feet per second 
for the majority of the numerical work; however, a series of problems is 
included varying the seismic velocity from 1,000 feet per second to 3,000 
feet per second in order to study the significance of this parameter. 

The blast loading is assumed to result from a contact surface 
burst of a one-megaton nuclear weapon. Overpressures ranging from 50 psi to 
1,000 psi are considered. The shock wave properties are discussed in 
Section 8.2. 

Approximate calculations based upon reference (16) indicate that 
the overpressures selected occur outside of the crater and therefore are in 
agreement with the loading assumptions discussed in Section 2.4. Structures 
located at the higher overpressures will probably be inside the fireball and 
will also be subjected to heavy doses of nuclear radiation. Again, as dis- 
cussed in Section 2.4, these effects must be considered in any actual design; 
however, they will have little or no effect upon the structural response of 
the arch. 

The majority of the numerical solutions are computed assuming a 
ratio of horizontal to vertical overpressure in the soil of one-fourth. In 
addition a series of problems is included in which this ratio is varied from 


one-fourth to one. 
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8.2 Shock Wave Properties 

The shock wave parameters used in this study are taken from ref- 
erences (16), (17) and (20). A summary of the more significant parameters 
is presented in Table 8.1. There is some disagreement among authorities (17+ 20) 
on the duration of the positive phase of the shock wave in the higher over- 
pressure regions. The differences in the values given are not great, however, 
when one considers the actual variation in all shock wave parameters to be 
expected in the field. Further, small differences in positive phase duration 
are not significant in this study because of the rapid response of the arch 
(cf. Appendix D). However, the rapid response of the arch does make the 
system quite sensitive to the initial shape of the blast wave. Therefore, 
while the positive phase duration for the 1000-psi shock wave may be taken 
as the same as that for a 100-psi shock wave without significantly affecting 
the results, the initial shape of the two shock waves must be defined 
carefully. This latter aspect of the problem is discussed in detail in the 
following paragraphs. 

The results of a study of the configuration of the shock front in 
soil are presented in Appendix B. 

The decay with time of the overpressure behind the shock front is 
discussed in reference (20). Later studies of shock wave characteristics 
are available but are not used here because of their security classification. 
Reference (20) gives the following approximate expressions for the variation 
with time of the overpressure at a point behind the shock front during the 


passage of the positive phase of the shock wave. 


= ’ 
Ore fia £: p(t')/p.. = Mine op aft /t4) (222) 
= a 7d sO d 
aa 
45 psi < isere < 700 psi: Q= ay + Isc, (t'7t,) (8.272) 
where t' = elapsed time at the reference point after the shock front 


has passed 


p(t') = overpressure at the reference point at time t' 





i bs 
De. = peak value of the overpressure as the shock front passes 
the reference point 


t, = duration of the positive phase of the shock wave at the 
reference point 


e = natural logarithmic base = 2.71828... 


al _ a : 0.0332(p.o/P.)°s Pso eS» 
7 *s0 
= P )({0.88 + 0.072 P_)); < 150 psi 
ie . ae - ee ii > - 7 
50 oO SO O SO 


peel + 0.1845(p_/P,) - 104/((p../P,) + 10) 


P = ambient atmospheric pressure (taken as 14.7 psi for 
numerical work in this study). 


The overpressure-time relationship for 50 psi is obtained by sub- 
stituting p.. = 50 psi and t, = 1.60 seconds (cf. Table 8.1) into Equations 


(8.2.1) and (8.2.2). The resulting expression is given as Equation (8.2.3). 


p(t’)/50 = (1-t'/160) 7 (0. 705+3.83/(1+t') ]t/1.60 (8.2.3) 


The computer program used in this study for the solution of the numerical 
examples (cf. Section 8.3) approximates the "actual" overpressure-time 
curves with a series of straight-line segments. These line segments may be 
obtained for the 50-psi shock wave by plotting Equation (8.2.3) and selecting 
a series of straight-line segments which approximate the curve to a reason- 
able degree over the time interval of interest. The coordinates of the end- 
points of the lines selected are given in Table 8.2. 

The overpressure-time relationship for 100 psi may be obtained by 
substituting p_. = 100 psi and t, = 1.48 seconds (cf. Table 8.1) into 
Equations (8.2.1) and (8.2.2). The resulting expression is given as 


Equation (8.2.4). 


32 t 
p(t')/100 = (1-t'/1.48) By ROSE 1+35.78(t'/1. 1 (8.2.4) 


Proceeding as in the case of the 50-psi shock wave gives the lines defined 
by the overpressure-time coordinates listed in Table 8.4 as an approximation 


to the overpressure-time curve. 
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Since the 1,000-psi shock wave falls outside the range in which 
Equations (8.2.1) and (8.2.2) apply, the overpressure-time relationship for 
this shock wave may be obtained by interpolating in the curves of references 
(20) and (23). This was done and a series of straight-line segments selected 
which approximate the "actual" overpressure-time curve reasonably well. ‘The 
line segments selected are defined by the overpressure-time coordinates 
listed in Table 8.4. 


8.3 Description of Numerical Examples 

The properties of the arches selected for the numerical examples 
are given in Section 8.1. Table 8.5 indicates the variations in parameters 
applied to each of these arches. 

The lower value of overpressure to which the reinforced concrete 
arch is subjected is well below the level which should cause buckling or 
general yielding of the arch due to thrust in the arch rib (cf. Section 8.4). 
The higher value of 1,000 psi is selected to magnify any tendencies of the 
arch to respond laterally. It is not intended to be below buckling or 
general yielding limitations. As pointed out in Chapter IX, however, the 
arch rib contributes relatively little to the lateral resistance of the 


arch-soil "structure" as compared with the soil contribution. Therefore, 


the 1,000-psi results are indicative of the response of an arch with a 
sufficiently large cross section that buckling and general yielding are no 
donger critical. 

The 50-psi loading on the corrugated metal arch is reasonably 
close to the overpressure causing general yielding due to thrust in the arch 
rib for the assumed soil and loading conditions. In fact, the case of zero 
depth of cover over the crown could very well be ruled out of any practical 
consideration on this basis alone; however, it is included in this study for 
comparison purposes. The 100-psi and 1,000-psi loadings are used to give a 


comparison with the corresponding loadings on the reinforced concrete arch. 





A 


Further, the observation made above in the case of the reinforced concrete 
arch, regarding the relative contributions of the arch and the soil to the 
lateral resistance of the system, is even more applicable here. Therefore, 

it may be stated that if the arch is strengthened sufficiently to withstand 
the thrusts developed in the arch rib, the response results of the 100-psi 

and 1,000-psi studies should be indicative of the response of the strengthened 
structure. 

Integration of the equations of motion for each of the numerical 
examples was done on the Illiac, the high-speed digital computer at the 
University of Illinois. Structural Research Laboratory Code 46624) was 
adapted to perform this integration. 

Tabular comparisons of the results obtained from the examples 
listed in Table 8.5 are presented in Appendix C. Response curves and con- 
figuration diagrams for the arch at various critical times are presented 


in Appendix D. 


8.4 Sample Calculations 


Calculations are presented here for the reinforced concrete arch 
with a depth of cover over the crown equal to one-fourth of the arch span. 
The arch and soil cover are shown in Figure 8.1. The arch is subjected to a 
100-psi shock wave from the explosion of a one megaton nuclear weapon as a 
contact surface burst. The seismic velocity of the soil is assumed as 3,000 
feet per second and the horizontal overpressure induced in the soil by the 
shock wave is assumed to be one-fourth of the vertical overpressure. In 
general, the calculations are based on a one-foot width parallel to the longi- 


tudinal axis of the arch (cf. Section 2.1). 


en 
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8.4.1 Calculation of Parameters and Coefficients 


The horizontal displacements of the haunches, hy and h 


are 


bid 


chosen as the independent displacements. The remaining displacements are 


related to them by Equations (3.2.7). 


Ds 


ne 


i 


tl 


The total 


(1-be/ad )(h, +h, ) 


NMJ Nr 


(2-(2.440)°/ (5.893 )°){hy+ hy] 


© 


4145 (+ hs) (8.4.1) 


- (b/a)h, 
: (2. 440/5.895 )h, 


- 0.41741 h, (8.%.23 


5 (a/e-b/a)(h, -h, ) 


% (5.893/2.4N0 - 2.4H0/5.895 (ny -h,) 


1..0006(h, -h, ) (8.4.3) 
(b/a)h, 

(2. 440/5.893 )h, 

0.4242 hs (8.4.4) 


rotations of the joints are related to the independent 


displacements by Equations (3.2.11). 


at 


Nl N{r 


(2/at+1/c+b/ad) 
(2/5.893 + 1/2.440 + 2.440/(5.893)*] 
. 4098 ere 


i/e 
1/2. 440 
0.4098 et71 


Oo 
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“SS : (1/c-b/ad) 
= = (1/2.4u0 - 2,440/(5.893)*] 
= 0.1698 ft > 
¥) = Ajh)-A,h, 
v= 0. 4098 h, - 0.1698 h, (8.4.5) 
V5 = shh) tots 
V5 = - 0.4098 h, + 0.4098 hs (8.4.6) 
V5 = Ash) -Ajh, 
Vs = 0.1698 h, - 0. 4098 h, (8.4298 


The general form of the moment-rotation relationship for the model 
hinges is shown in Figure 8.2. The dimensions on this figure may be found 
using Equation (3.3.1). The material properties for the arch are listed in 
Section 8.1. Figure 8.3 is a sketch of the cross section of the arch rib. 


Now, in the derivation which follows, neglect direct stress in the arch rib 


and let: 

he = thickness of the concrete arch rib 

ds = effective depth of the concrete arch rib 

bo = width of the arch rib parallel to the longitudinal 
axis of the arch 

j = ratio of the distance, jd_, between the resultants of the 
compressive and tensile stresses to the effective depth, qd. 

k = ratio of the distance, kd _, between the extreme fiber 
and the neutral axis to the effective depth, qd. 

I. = gross mement of inertia of the arch rib cross section 

I, = transformed moment of inertia of the arch rib cross section 
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average moment of inertia of the arch rib cross section 
area Of reinforcing steel in a one-foot width of arch 


28-day strength of a standard 6" x 12" concrete test 
cylinder 


yield strength of the reinforcing steel under dynamic 
loading 


modulus of elasticity of the reinforcing steel 
modulus of elasticity of the concrete in the arch 


modular ratio for steel and concrete. 


From reference (22): 


Then: 


n 


5 +10, 000/f! 


5 + 10,000/3,500 
7.86 


E./n 


30 x 10°/7.86 


3.82 x 10° psi 


l 


y = Asfay dd, 


0.24(50)(7/8)(4)(1/12) 


4.50 ft-kips 


b(n, )?/12 


12(8)?/12 


pid ine 


5 2 2 
db. (xd. ) /3 + nAd/(1-k) 


22 § (5/8)(4)}> + 7.86(0.24)(4)°(1-3/8)° 


Coe ine 
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1 
Ea ™ 5 (I+ ie 


Caley 


268.6 ee 


K,= EI /L, 


(4.62-x 10°) (268.6)/(6.378)(12) 


" 


1117.2 ft-kips/radian 


<— 
Ht 


: M/K, 


3.50/1117.2 


0.003135 radians 
The relationships shown in Figure &.2 are expressed algebraically in 


Equations (8.4.8) 


OS V3 5 Yy! "3 Beg 
(8.4.8) 
vy Se M, = M 


Note that the value of ¥, must be adjusted for any "permanent set” in the 
hinge as the structure responds. Also, Equations (8.4.8) are valid for 
rotation in either direction. For convenience, the inequalities at the 
left of Equations (8.4.8) are shown for the positive direction of rotation 
only. 


The equivalent arch resistance may be found using Equations (3.4.7). 


q, = (1/L*)(atc) 
= (1/(6.378)°][5.893 + 2.440] 
= 0.2049 et 

a, = (1/2L°)(a/e-b/a) 


[1/2(6.378)*][5.893/2. 440 - 2.440/5.893] 
2 


0O.02K6 ft 








ee 


c/L° 


2.4h0/(6.378)* 
ab 


2 
i) i) 


0.0600 ft 


(1/2L°)(a/c+b/a) 


(1/2(6.378)*][5.893/2. hho + 2.4h0/5. 893] 
2 


Ya 


Uy) 


0.0348 ft 


M, (q+) -a,h5 ) - M3 (4, +438) -a,h;) 
- M3 (q.+a,hy -ayb,) + M3(a,+a,h) -agh, ) 


M, (0. 2049 + 0.02h6 ae 0.02h6 h; ) 


PP 


- M,,(0.0600 + 0.0348 h, - 0.0246 h,) 


- M, (0.2089 + 0,026 h, - 0.02H6 h,) 


% 


(8.4.9) 
% M,,(0. 0600 + 0.02h6 h, - 0.0348 hs) 


From the results of the study presented in Appendix A one may 
conclude that the hydrostatic pressure ratios for soil pressures corre- 
sponding to the fully active and the fully passive states of stress in the 
soil may be calculated to a reasonable degree of approximation using the 


Coulomb formulas (cf. reference GIEQ DF 


2 
ne 60s ee (8.4.10) 


A ie Vin ae Se 
4 sin( 04+? )sin(P-i 
cos’ cos( 5+) E +) 2 SPRY SS 


where K hydrostatic pressure ratio corresponding to the "Active 


State of Plastic Equilibrium" in the soil (cf. Section 2.2) 


? = angle of internal friction of the soil 


W = angle of inclination with respect to the vertical of 
the lower bars of the model 


6 = effective friction angle between the arch and the soil 


i = slope of the ground surface with respect to the horizontal 
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cos“(30° - 22°29.6' 
cos” (22°29.6' )cos(20°+22°29.6") E A sin(20-+30 )sin 


V cos(20°+220°29.6' )cos(22°29.6') 


Ky = 0.510 
2 
K, a Be Se) con) (8.4.11) 
. Soa ae Va 
2 . gin(P+5 )sin(P+1) | 
cos" cos(w-5) E a) eoulues lcoamee 
where = hydrostatic pressure ratio corresponding to the "Passive 





State of Plastic Equilibrium" in the soil (cf. Section 2.2). 


a sok 0°420°% 29.6! 
i 


Oo ,,0 O 
cos“ (22°29.6' )cos(22°29.6'-20°) E -/- a ‘ — sae = | 
cos(22°29.6'-20° )cos(22° 29.6') 
= 5.453 
The soil pressure ratios corresponding to the directions of Pa and Po for an 


"at rest" state of stress in the soil (cf. Section 2.2) may be found using 


Equations (4.1.11) and (4.2.11). 





Ko, = Vv 1aKé cos [ (45 )-arctan(1/k, )] 
= \/1+(0.5)* cos [(22°29.6'+20°)-arctan(1/0.5)] 
= 1.051* 
fo ae : 
Kop = V1+k- cos [arctan(1/k, ) - (w-8)j 


(etna) cos farctan(1/0.5) - (22°29,6' -20° )} 


ess: 
The deformation corresponding to full development of the "active" 
state of stress in the soil is taken as an inward rotation of the lower bar 


An error of 1.5% in K a iS introduced at this point. It was found to cause 
a change in the calculated response of less than 1%. 








12h 


of the model, S py? of 0.001 radian (cf. Section 4.4). The deformation corre- 
sponding to full development of the "“passive' state of stress in the soil is 
then taken as an outward rotation of the lower bar of the model, Ey? as 
given by Equation (4.4.8). 

(CE - “OP ) 
ore 

3.4353 - 0.543) 
TsO5n. =10.510 


Ely = & 


0.001( 


it 


0.005342 radian 


While Shy is not known to this degree of precision, it is convenient to 
retain a sufficient number of significant figures that calculation of the 
slope of the "elastic" line in the soil resistance-deformation relationship 
(cf. Figure 4.12) will give the same result for the build-up of both the 
“active and the "passive" states of stress in the soil at a given over- 
pressure level as assumed in Section 4.4. 

Expressions for the soil "stiffness" may now be obtained from 


Equations (4.4.10) and (4.4.12). 


y) 
Ya 
pe Za (Kog-K,) (cos (048 )+(b/a)sin(m45 )] 
Ly 


> 
0. 120 2, OO 1.0 1-0. 10 0 1 Oo 
= “B93 )(0 001 [cos(22°29.6'+20° ) 


+ (2420 )sin(20°29.6'+20°)| 


4647 kips/radian 


nm 
N 


(3/7 a. Mk! 
(3(3647)/(0.120)(12. 500) J (0. 144] 


1050 kips/radian-psi 
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Ke = k' =oh Da 

7 = 4647 - 1050 Ps j (8.4,.728 
Ke =k! - K'PS 3 
Kj = 3647 - 1050 Pe (8.4.13) 


Alternatively, the soil stiffness may be expressed using Equations (4.4.14) 
and (4.4.15). 


K, ,/a = (k'-k" Pp, ,/a 
= (3647-1050 p, ,)/5.895 
K,,/a = 618.9 - 178.2 Ps j (8.4.14) 
Ki /a = (k'-k" p ,)/a 
Ki /a = 618.9 - 178.2 Ps (8.4.15) 


The soil resistance-deformation diagram for this problem is shown in 
Figure 8.4. 

The mass concentrations corresponding to the horizontal and vertical 
components of the inertial forces acting at the joints of the model are given 
by Equations (5.1.13). The geometry of the arch-soil "structure’ is shown 
in Figure 8.1. The geometry of the soil masses assumed to move with the 
arch is shown in Figures 5.1 and 5.2. The limiting vertical dimensions of 
the soil masses at the haunches are determined using Equations (5.1.3) and 
2.4). 


d= ((a/2)ecta,] < (b+) 


? 
( (22528) + 2.440 + 4.167] < (2.440+5.893) 
9.554 > 8.333 
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*, use : ql =o ene 
do = [(c/2)+a ]< 22 


? 
[(2.440/2) + 4.167] < 2(5.893) 
5 0t = ll. OO 


°, use: do Be SOU Ete 


la @(a,- 8 2)(PES) + (a- SEE) +ae, cot(y + 2) | 


i 
E 


y ne 

0.120 2. 440 5895 /2.uh0 + 5.895 
B(32.2) 8.533 Bat (0959 2 ) 
+ (8.555 - 2822 + 2.) (2-895) 5 (8.333 )Pcot 60°| 





0.150(8 
2\) 


6.378 
i 2 


2, 


0.1782 k-sec“/ft 


7 
(owe = Z| yy BI(dag- SPE) + (agg SENG) + a2, cot(G - 2] 


+ 7 noe 


il 





2(32.2) 
+ (8.533 - 2928 + BNO) (3-095) « (8,533) cot 30° 


| | 895 )/2.4N0 + 5.893 
‘a  8.393¢2 iO) 4 (8,333- 2+825)(2 HO + 9.895) 


, 215008 
eld 


6.378 
T2(52.2 


0.5276 k-sec“/ft 








esl 





7. ; hl 

oa Ba one - 4) + 
_ 0.120(5.893) 2.440,  0.150(8)(6.378 
: 32.2 (5.387 - |) + “3505500 


= 0.12%7 k-sec“/ft 


(min = (,),,, 


= 0.1782 k-sec"/ft 


7 (7a) ) ot 
= 0.3276 k-sec“/ft 
y 
my = 52 [4 @)r(a,g- S12) + (a,,- SV) | + 


= BPs aeSy | 8.35322) + (8.535 - 22898) BHO + 5.895) 


a 
5.893 + 2.440, 5.893 
+ (8,333 - 282 + 2.) (2.8 )| + Bagh. 378) : 


PyGht@55) sce) ft 


ae i 


5 2 


= 0.1247 eee 


mon 


0.1035 k-sec“/ft 


The equivalent masses for the two-degree-of-freedom system can 
now be found using the mass concentrations calculated above and Equations 
(5.2.6). The calculations for these parameters can be shortened consider- 
ably by noting that (m5. = (ms )ey (ne (mut? m, = me and 


Vv = VY ; 
l,m 3,m 
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7 me be\2, by2. 5 pa bY * 
"1,al_, = 1.2 7 (my in® oat : oa? +m), (2) : tO - a) 
ot 
= 0.1782 + ati . axel rer OF 1035 (= ey 
(5.893 ) 
0.124 8 2.440 ; 
+ SES (Senko ~ 52855) 
Eres on eenee-/ Ft 
yal, 5, 7 Crows * BO - Braye BEB! 
2 
= 0.3276 + Od etT E - (2.0) - 0.1055(2 0)" 
(5.893) 
0.1247 ,5.8 2.440 ,° 
+ SSE (S85 > S895) 
= 0.4916 k-sec“/ft 
me a by 
Vv = = (1 - ee st (==) 
1,m| ee) i + oo c.6hOUa 
0.12h (2.440)" | | 0.124 
nom a -& (e.mo)" | , es a (ae - 5 a 
7" 893) 
= 0.1034 k-sec“/ft 
mn =e 
3,m| ee 1,3 1,m| a wl 


0.4916 k-sec"/ft 


id = H 
3,m| = 2,4 1,n|,, = 1,2 


= 0,400 k-sec“/tt 


“Lal, 


"3 ym : 
, mn = a ot = 1,2,5,% 


= 0.1034 k-sec@/ft 





The equivalent masses for the 


may now be written. 


Asymmetric-Right: ce 


+5 
oymmetric-In: My 5 
Ms sp 
Symmetric-Out: My 
os 


Asymmetric-Left: ce 


eo 


four motions of the system (cf. Figure 5.3) 


i 


O 


0 


The equivalent loads on 


by Equations (6.3.2). 


i 


.3402-0.1034 h,/hy 
—— (8.4.16) 

.4916-0.1034 h,/h, 

.3422-0.1054 h,/hy 
a (8.4.17) 

.3422-0.1034 h./h 

Eas 

-4916-0.1034 b,/h, 
a (8.4.18) 

.4916-0.1034 h, /h 

5 

4916-0. 1034 h/by 
(8.4.19) 


.3422-0.1034 h,/h, 


the two-degree-of-freedom system are given 


J = 5K (atc) + =(b/a)(b+a) 


Ae 


= 0.144 BG (5-695 + 2.4W0) + HEE)(2, uo + 5.695) | 


= 0.3985 k/psi 


1 be 
J, = cK (1 - = 


= (2. 4W0)(F)(O.144)[2 = (2.480)°/ (5.895 )*] 


= 0.0182 k/psi 


J, = = a(a/c - b/a) 


0.8491 k/psi 


2=825. (5.893/2. Who - 2.440/5.893 )(0.144) 


an) ee 
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=- J,P, - J,(py, - Pa) + JzP5 


Fal 


- Jn(Po, - Pop) - Ue ee 


OF! 


br] 
i) 


- 0.3985 Py . 0.0182(p,, Poy ) + 0.8491 a | | 
8.4.20 


Si 
i) 


- 0.0182 (p,, - Py) - 0.8491 p, + 0.3985 Ps, 


Since Equations (8.4.12) through (8.4.15) for equivalent soil stiff- 
nesses and Equations (8.4.20) for equivalent loads are written in terms of 
the average overpressures acting on the loading regions defined in Section 6.1, 
it is necessary to find these overpressures using Equations (6.2.1) and the 
proper set of time parameters from Section 6.1. Table 8.1 gives the velocity 
of the shock front at the 100-psi overpressure level as 2,950 feet per 
second. This is approximately equal to the assumed value of 3,000 feet per 
second for the seismic velocity of the soil. From the results of the study 
presented in Appendix B, one may assume, for the depth considered in this 
example, that the shock front in soil is essentially vertical and a continua- 
tion of the shock front in air. This corresponds to Case I, Section 6.1.1. 


The time parameters for this case are given by Equations (6.1.1) and (6.1.2). 


(b+24 cot i 3 j)/e0 


0a 
= = [2.440+2(12.500)cot 30°}/2(2950) 
= - 0.0078 seconds 

= al 

rot = 3 (tog: , top) 
z = (-0.0078+0.0013) 


= - 0.0032 seconds 









2 (b-2ba _/ a)/2u 
- [2.440-2(2.440)(12.500)/5.893 ]/2( 2950) 






+ 0.0013 seconds 










- (b+2d cot(T 3 iS )/2u 


(2.44042(12.500)cot 60°}/2(2950) 






0.0029 seconds 







ct | 
ul 





(+ att 


O ae 





( -0.0029+0.0013 ) 


oO 
Mir Nir 






0.0008 seconds 






ct 
uN 


(t bo tat 


ab 1p) 


(0+0.0014) 


+ Dik Ir 


0.0007 seconds 
ce. 3 (b+a )/2U 
= (2.440+5.895 )/2( 2950) 


= + 0.0014 seconds 


Van Satay 


0.0014 seconds 


i) 
+ 


(t e 


2a + top) 


(+0.0014+0.0034) 


— NS ne 


0.0024 seconds 











2b 


ca 


ct 


eb 


4a 


ct 


5b 
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(b+3a )/2U 
[2. 4h0+3 (5.893 )}/2( 2950) 


+ 0.0034 seconds 


(to. i t,) 


ble 


5 (+0.00144+0. 0024) 


+ 0.0019 seconds 


(b+2d )/2U 


(2. 4h0+2(5.893 ) | /2( 2950) 


0.0024 seconds 


+ 


(t,, a i) 


(+0.0024+0.0034) 


+ Nik le 


0.0029 seconds 


Yop 


+0.0034 seconds 


(t,o ie ts) 


KH nOlH 


5 (+0.0034+0. 0048 ) 


+ 0.0041 seconds 
(b+2a )/U 


(2.440+2(5.893 )]/2950 


+ 0.0048 seconds 


(~b+2B-2bd /a)/2U 
[ -2.440+2(16.667 )-2(2. 440 )(12.500)/5.893 ]/2( 2950) 


+ 0.0035 seconds 




















1 
‘ 5 (ty, bd ti) 
= = (+0.003540.0077) 


+ 0.0056 seconds 





= (~b+2B+2d cot (f - x) )/2u 
[-2. 440+2(16.667 )+2(12.500)cot 60° ]/2(2950) 


+ 0.0077 seconds 


(ty, e type) 


=5 (40.0035+0.0126) 


= + 0.0080 seconds 


tip! = (-b+2B+2d cot (ft - £})/2u 
= [-2.440+2(16.667)+2(12.500)cot 30°]/2(2950) 


= + 0.0126 seconds 


t = - b/2U 
= - 2,440/2(2950) 


= - 0.0004 seconds 


7 
Pepe a (ten. * tay? 
2 = (-0.0004+0.0052) 
= + 0.0024 seconds 
toy = (-b+2B)/2U 


= (-2.440+2(16.667) ]/2(2950) 


+ 0.0052 seconds 











134 


The solution of Equations (Se2.)) may be presented conveniently 
in tabular form as shown in Table 8.6. Since the integration of the equa- 
tions of motion is to be performed on the Illiac, each of the overpressure- 
time functions is approximated by a series of straight-line segments. 
Therefore, the overpressure-time coordinates given in Table 8.3 are used for 


the general overpressure-time function, p(t'), with the exception that the 


maximum value of each overpressure, i.e. Sore, = (ea), is taken 
directly from the p(t') vs. t' curve itself. This is easily done and improves 
the approximation slightly. Thus, Table 8.6 contains the coordinates of the 
end points of the series of straight-line segments used to approximate the 
average overpressure acting on each of the loading regions. 
The equations of motion may be written using Equations (7.1.4). 
Fo ae) 


m i ey ne 
1, 1,m"3,m ier 5 


ae 0.4916 


dey: (0.3422)(0.4916)-(0.1034)* 
2 
= 4.1213 ft/k-sec 
Y 
2 
Gm 5 


UM Hd = V 
mS ole dt 


0.1054 
eG 2 
) (0.3422 )(0.4916)-(0.1034) 


0.6565 ft/k-see" 
v 
ey ee) eee 


: Uh -V V 
29M HD a3 m7’ 1,m"3,m 


0.1034 
234 = 2 
, (0.3422)(0.4916 )-(0.1034) 


0. 6565une/k-see” 





diye 


u 
G = ae eee 
oe “i,m 3,m-1,m‘3,m 


7 0.3420 

G. ie = ° 
’ (0.3422 )(0. 4916 )-(0.1034) 

2 


2.1727 ft/k-sec 


The remainder of the coefficients are calculated in the same manner. 
Inserting the coefficients into the equations of motion for the four motions 


of the system (cf. Figure 5.3) gives: 


Asymmetric-Right: 


h) = ple (ie Ode) ns 0.6565(H, -2,-@.,) 
F a = _ _ (8.4.21) 
h, = 0.6565(H, -Q, -@.,) + 2.1727 (H,-@,-Q.5) 

Symmetric-In: 
hy = 5-2159(H, -@,-Q,,) + 0.9717(H, -@,-@.;) 
c _ _ _ _ (8.4.22) 
hs = 0.9717(H,-@,-Q.,) + 5-2159(H, -Q,-@.,) 

Symnetric-Out: 
h, = 2.1283 (H, -@, -@., ) + 0.4477 (H,-@,-.,) 
> i es - = (8.4.23) 
h, = 0. W477(H, -Q, -@_,) nu 2.1285 (H, -@, -@.) 

Asymmetric-Left: 
h, = DWE TUTE RS Oe ) a 0.6565(H, -@, -Q.,) 
. 7 _ _ (8.4.24) 
h, = 0.6565(H,-@,-@,,) + 5.1213(H,-@,-@,,) 


8.4.2 Investigation of General Yielding and Buckling of Arch Rib 


Before proceeding with the integration of the equations of motion, 
it is advisable to investigate the effects of the thrust in the arch rib. 


The average overpressure acting on the arch is given in Table 8.6 as Pa 


= eee 
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The thrust in the arch, T, caused by this average overpressure is given 
approximately by Equation Gaon which is repeated here for convenience as 


Equation (8.4.25). 


-2 lt: 
t= 5p, 5B (8.4.25) 
where Ps = effective uniform radial overpressure acting on the arch 
B = span of the arch 


Since the value of or being used in this problem is one-fourth, the radial 
overpressure acting on the arch will vary from 25% of P, at the abutments to 
100% of >. at the crown. Thus, because the arch is semi-circular, this gives 
an average radial overpressure on the arch of 62.5% of the vertical overpressure. 
However, this neglects the thrust induced in the arch as the sides are forced 
Outward against the soil during lateral response. Therefore, some higher 
proportion of the vertical overpressure is appropriate. For the approximate 
calculation here, pt is taken as three-fourths of Pa 

A measure of the capacity of the arch rib cross section to carry 
thrust loads may be obtained from Equation (8.4.26). 


T = 0.85f" Pet foe (8.4.26) 


J y 


where = thrust corresponding to general yielding of the arch rib 
fl 28-day strength of a standard 6"x12" concrete test cylinder 
A. = area of concrete in the arch rib cross section 
A = area of steel in the arch rib cross section 
ge dynamic yield strength of the reinforcing steel 


Equation (8.4.26) neglects bending moment in the arch rib and is therefore 
approximate; however, because the mode of response being considered here is 
essentially a breathing mode, and further because of the high degree of soil 


restraint, bending in the arch rib is not of prime importance in this 
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investigation. Substituting values into Equation (8.4.26) gives: 


Ak 


: 0.85(3.5)(96.00-0.24) + 50(0. 2h) 


297 kips 
The average overpressure, Pes corresponding to general yielding 
of the arch rib may now be found using Equation (8.4.25). 


(p') 


& max 


(Pees On45 


es 
men 
Oi 


_ eke 1000 
~ 0.75(16.667 )(1 


550 psi 


Since 330 psi is greater than the maximum value of P, in Table 8.6, i.e. 
P, = - 97.7 psi, the arch should not fail by general yielding. 

A qualification should be made here regarding the use of an 
“average radial overpressure in this case. Since the crown is subjected to 
an overpressure level greater than the "average" value, local yielding may 
occur which could lead to general failure of the arch. This is particularly 
true when the depth of cover over the crown is zero or very small. 

In considering the elastic stability of the arch rib, one must 
keep in mind that buckling is not instantaneous. For the higher yield 
nuclear weapons, however, the duration of the positive phase of the shock 
wave is quite long when compared with the characteristic response time of 
the structure so that buckling is a possibility which should be investigated. 
If one assumes that the arch rib may be approximated by a bar pinned at the 
crown and at the abutment and restrained by the soil, the elastic stability 


of the arch rib may be investigated using Equation (8.4.27). 629) 
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hk 
t“EI 2 Pt 
Ti = oe We +) 
L mt ET 
h 
BR 
t= 22 (n? + *—) (8.4.27) 
R 16m°EI 
where Dee = thrust corresponding to buckling of the arch rib into a 
simple sine curve 
KE = modulus of elasticity of the arch material 
I = moment of inertia of the arch rib cross section 
£ = developed length of one-half of the arch 
m = number of half sine loops in the buckled configuration 


By = soil stiffness, here assumed to act in both directions 


R radius of the arch. 


ul 


Equation (8.4.28) may be obtained from Equation (8.4.27) by equating the 
critical thrust values for which the number of half-waves changes from m 
to m+l. 
m*(m+1)° = ul (8.4.28) 
16EL e “Te 
The minimum soil stiffness corresponding to the maximum thrust 
in the arch is obtained by noting that the maximum value of Pa in Table 8.6 
occurs at t = 0.0052 seconds. While the configuration of the arch at 
t = 0.0052 seconds is not yet known, a “lower bound” for the soil stiffness 
may be obtained by using the minimum soil restraining pressure acting at 
this time. From Table 8.6 it appears that Pas has an absolute magnitude, 


i.e. 18 psi, less than all of the other soil restraining pressures at 


t = 0.0052 seconds. 


0.0052-0.00 
(550126-0.0055) (96-4) 


~ 
Paz 18 psi 


z Pia 











oe 


The soil stiffness is now determined using Equation (8.4.15). 


Ky /a = 618.9 - 178.2 Ps 
K,;/2 = 618.9 - 178.2(-18) 
= 3827 x/ft* 
B, = 38-7/0.144 
= 26,600 psi 


The number of half sine loops corresponding to these conditions is obtained 


from Eyuation (8.4.28). 


alae = 


Z 26 ,600(100)" 


16(3.82x10~ )(268.6) 


162.03 


m(m+1) 12.729 


3.103 


m 


The critical thrust may now be found from Equation (Sie) by substituting 
the integral value of m which gives the lowest value of beet The value m = 4 
gives the minimum value in this case. (Note: Differentiating Equation 
(8.4.27) with respect to m, setting av _/dm = O, and solving for m gives 
m = 5.568). 
p R’ 
T = k 

oe Ro 16m“EI 


_ 4(3,82x10°)( 268.6) (16 + 262.03) 
. 1 


(100 )©(1000 ) 





) 


10,725 kips 
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The average overpressure, Po» corresponding to buckling of the arch may now 
pe found using Equation (8.4.25). 


vad 
ex 


RPE) se 0.75B 


- 2(1072 
satan 


11,900 psi 





This approach, although approximate, indicates that the arch should not 
buckle elastically prior to general yielding of the arch rib itself. Note 
that this is an "“elastic’ analysis so that the actual buckling of the arch 

in this case is associated with general yielding of the arch rib as described 


in the first part of Section 8.4.2. 


8.4.3 Integration 


Since it has now been determined that neither buckling nor general 
yeilding of the arch rib should occur under the given loading, it is possible 
to proceed with the integration of the equations of motion. The results 
shown in Appendices C and D were obtained using a digital computer (cf. Sec- 
tion 8.3). Two of the problems were also integrated by hand in order to 
check the computer solutions. 

A cycle of integration is given below to illustrate the operations 
involved. The numerical integration formulas used are those derived by 
Professor N. M. Newmark and presented in references (4), (5), (6) and (7). 
Equations (8.4.29) are expressions for the velocity and displacement of mass 3 
at the end of time interval At in terms of the acceleration at the end of the 
interval and the acceleration, velocity and displacement at the beginning of 


the interval. 
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. Q e 3t oe oe 
Dj, ntl 7 Dyn ae ‘he . Sal 

. al se a (8.4.29) 

= a t 

bs nal = hy gthy ploti+s B Mh, y(t) 

ee eo 

U] 
ee By nay *) 
where h. ’ h, mele = acceleration, velocity and displacement, 
JN J, ’ 


respectively, of mass j at time t = to 


= acceleration, velocity and displacement, 


h wae ,» A 
jyntl® “J, ntl® ~j,ntl respectively, of mass j at time t=t iy 


At time interval = t a 


n+l “n 


BE parameter introduced to indicate how much of the acceleration 


at the end of the time interval enters into the relations for 
the displacement at the end of the time interval. 
Since it is desired to use a linear approximation to the acceleration curve 


within each time interval, the value B' = 1/6 is used in Equations (8.4.29). 


° e At ee on 
een . ea ' (ha ) bs nar? 


(8.4.30) 


by ntl =h, + hy (ot)+(2/3)h, (at) +(1/6 dh, (At ) 


J,n 


ynt+l 


In the following illustrative example a time interval, At, of 
0.0001 seconds is used. This time interval was actually used for the computer 
solutions in order to detect the rapid changes in the loading as the shock 
wave traverses the structure. Solutions were also obtained using time 
intervals of 0.001 seconds for several of the problems which had more slowly 
varying load-time functions in order to check the significance of round-off 
error. These solutions were nearly identical with the solutions obtained 
using At = 0.0001 seconds. 

For illustrative purposes it is preferable to consider a time 
interval other than the first. Therefore, for this example three time 


intervals have already been integrated, so that at time ve = 0.0005 seconds, 
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the displacements, velocities and accelerations of the equivalent masses 


are 
h, = 40.397x10~ ft h. _ = 40,0GNMIO © Ft 
jipgal : 5,0 : 
hy = +3.927x10 ft/sec hs = +0.835x10~? ft/sec 
ae = +26.474 ft/sec” bs = +5.568 ety aac. 


Assume values of the accelerations at the end of the interval, i.e. at 


time t_447 0.0004 seconds, as follows: 


a 


h = +35,279 3) oe 


2 
ene) = +7.420 ft/sec 


ee em 


The velocities at t.., = 0.0004 seconds may be found using the first of 


ag 
Equations (8.4.30). 


by ad = (+5.972x1077) + 2 (426.471435.279)(10"") 


+7.060x10~? ft/sec 


Hie uesxior ft/sec 


eonnel 


The displacements at ca = 0.0004 seconds may be calculated using the 


a 
second of Equations (8.4.30). 


h (40.397x10 © )+(43.972x107 )(1074) 


1,n+1 
+(1/3 (26.474) (107?) +(1/6)(35.279)(107°) 
= +0. 941x107© ce 
o/s 
hy al = +0.198x10- ft 


The rotations of each joint are found using Equations (8.4.5), (8.4.6) and 
mo. 4.7 ). 


<— 
il 


0. 4098 h, - 0.1698 hs 


[0.4098(0.941) - 0.1698(0.198)]110 


6) 


40,5520 ant ane 
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-O, 405x10© radians 


<_ 
It 


+0.079x10 “6 radians 


¥5 


The resisting moments at each hinge are found from Equations (G2 43). 
me YTB 
= (40. 352x107) (1117.2) 


= +0.393x10~" ft-kips 
M. = -0.340x107? ft kips 


M, = +0.088x107 ft kips 


The equivalent lateral resistance of the arch is now found from Equations 


(3.479). 


M, (0. 2049+0. 0246 h, -0.02h6 h,)- M,,(0.0600+0. 0548 hy 


eal 


-0. 0246 hs) 


= [ 40. 395x107- }[0. 204940. 02h6(+0. 941x107) 
yl 


“6) 


~0.0246(+0.198x10 -~ [-0.340x107” } [0.0600 


+0.0348(+0.941x10 -0.0246(+0.198x10©) } 


= 40.101x10- kipes 


Q, = -0.038x10~" kips 


-6 


The displacement h = +0.941x10 ~ ft indicates that the windward haunch, 


1,n+1 
joint 1, is displaced inward into the arch. From Table 8.6 the average over- 


pressure acting on the "active" soil failure wedge adjacent to joint 1, p 
& s 





17 
at Oe 0.0004 seconds is 77.2 psi. 
_ 0.0004 =z : 
ss eh me 67.9 a 0.0013 (98.3 sf 67.9) 


T7-e psi 
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The soil stiffness corresponding to this condition is given by Equation 


fee. 12). 
Ke = 3647 - 1050 Ps j 
Kr = 3647 - 1050(-77.2) 
= 84,710 kips/radian 
The displacement h = +0,198x10°? ft indicates that the leeward haunch, 


3,nt+l 
joint 3, is displaced outward against the soil. From Table 8.6 the average 


overpressure acting on the "passive soil failure wedge adjacent to joint 3, 


» at ti45 0.0004 seconds is zero because the shock front hes not yet 


3 
Poa AG 
arrived at point 4a (cf. Figure 6.1). The soil stiffness corresponding to 


this condition is given by Equation (8.4.13). 


Koj 3647 - 1050 Psi 


K 33 


il 


3647 - 1050(0) 


3647 kips/radian 


HT 


The rotations of the windward and leeward pressure surfaces are obtained 


using Equations (4.4.13). 


I 


By h,/a 


‘ +0. 941x107 
5-095 
we 
+0.160x10 ~ radians 


3 = hs/a 
i 


_ +0.198x10 
; 5.893 


+0.034x107° radians 


un 
| 


Since these rotations are less than the yield values calculated previously, 
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Mae. , iy = 0,001 radians and Shy = 0.005342 radians, and no "permanent set" 
has yet occurred in the soil, the equivalent soil resistances may be calcu- 


lated as follows: 


Sea, Beal 
= 84, 710(+0.160x107°) 
Or ine 
853 = Kaz Ss 
= 3, 647(40.034x107°) 


40.124x107" kips 


The equivalent loads on the system may be determined from 


Equations (8.4.20). The average overpressures loading the structure at 


ia 0.0004 seconds are given in Table 8.6. 
-, = (20004)(99.4) = 28.4 psd 
Py = 0 
Pz = 0 
Po, : 
Poy : 


2 
I 


7] -0.3985 P) = 0.0182(p,,_-Po») + 0.8491 De 


= -0,3985(-28.4) 


i; 


The fact that both of the velocities, ected and eae at time 


0 


i 


(ie 0.0004 seconds are positive indicates that both haunches are moving 








1h6 


toward the right. Thus, the motion of the structure is "“Asymmetric-Right’ 
(cf. Figure 5.3) and the equations of motion for the system at this time 
are given by Equations (8.4.21). By substituting the effective load and 
resistance parameters calculated in this section into Equations (8.4.21), 
it is now possible to check the assumed values of the accelerations at time 


tel = 0.0004 seconds. 
hy na 7 35-1213 (H, ~@, -@., )+0.6565(H, -@, -@., ) 


4.1213 [11.32-(+0.101x10 )-(415. 55x10 -)] 


+0.6565[0-(-0.038x10~- )-(+0.124x107 )] 


+55.29 et Bees 


It 


AG o 
rs = +7.42 ft/sec 


Since these values are essentially the same as the assumed values, i.e. 


eo os 
ntl = +55.279 ft/sec and al 


not required and the velocities and displacements at t = 0.0004 seconds 


= +7.420 ft/sec, additional trials are 


which were calculated above are correct. Thus, the integration may now 
proceed to the next time interval. The results of the solution of this 


problem up to t = 0.09 seconds are given in Figure D.3. 


8.4.4 Rib Shortening 
The effect of rib shortening is illustrated in Figures D.24, D.25 


and D.26. Expressions for the "rib-shortening’ corrections are obtained using 
Equations (3.5.8). In the following derivations an average uniform radial 
pressure on the arch equal to 75% of the vertical overpressure in the soil 


is assumed (cf. Section 8.4.2). 





Sh, = 


uJ 


Sh 


oh 


Ooh, = 


Oh), = 


Sh. = 


Oo 
a 
i) 


Fond 


= ee 
2 Po 


~idtr 
2 Pa 
-11.29x1 


-8.47x10 


ay 7 es 
2 Pa AE 


AE(ad-bc ) 


16.667(6.578)" (2, 440 )(12) | 


8(3.82x10 )(5.893°-2, 4h0°) 
-6 | 
Pa 


BL“c 
@) 


cae 
Pa 


BL“a 
ad-be 


l=, 16.667(6.578)" (5.895 )(12) 
ee! ) 


8(3.82x10°)(5.893°-2, 440° 


27. 26x10 


20. 45x10 


a 


MIke 


54.55x10 


40.92x10 


20.45x10 


-6 rat 
Po 


6 — 
Pa 


=e (BL 
C cAE 


= ieee 
Po 


2.440(8)(4.82x10° ) 


Gas 
Pa 


- — 
Pa 


67— 
P 


a 
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(8.8.345 


(8.4.32) 


(8.4.33) 


(8.4.34) 


(8.4.35) 


(8.4.36) 
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The parameters required to plot Figure D.25 are calculated below 
in order to illustrate the operations involved. The time selected for this 
computation corresponds to the first maximum downward displacement of the 
crown. From Figure D.5, this occurs at t = 0.0198 seconds. The displace- 
ment components at this time may be calculated using the data plotted in 


Figure D.3 and Equations (8.4.1) through (8.4.4). 


h Terao me ft 


1 


eae +0.4145(h, +h, ) 
= +0. 4243 (-5.845.0)(1077) 


= Soule aon Lt 


he = +5-0x10- ft 


hy, = -0.4242 h, 
meeoui (-5.0r10- } 


eon et 


ih @ ile) seater ee 


he = Oe al h,, 
= +0, 4141(+5.0x107? ) 


= 42,1x107 ft 


The average overpressure acting on the arch at t = 0.0198 seconds 
may be obtained either from Tcble 8.6 or from a plot of p(t) vs. t. For 
this example a value of P, = -88.1 psi was obtained from a plot of p(t) vs. te 
This value is now substituted into Eyuations (8.4.31) through (8.4.36) in 


order to obtain the "rib-shortening’ corrections at t = 0.0198 seconds. 













dh, = -8.47x107° Pe 
= ~8.47x107°( -88.1) 


+0.7x107 ft 











= oh, 


-0.7x10° ft 


6h) = +20.45x10© P, 


= 420. 45x107°( -88.1) 


2 aGeaie 2 


Sh. = +40.92x107° Pe 
= +40.92x107°( -88.2) 


Semon ioe at 


She = oh), 


= aa 2 


The net displacement components, h,, of the joints of the model 
may now be obtained by adding algebraically the “rib-shortening’ corrections 


to the displacement components calculated above. 


h. =h. + 8h. 
J 


j j (8.4.37) 


= (-5.8+0.7)(107) 
= ~5,1x10 ft 
1S (-0.440)(1077) 


i 


~0.4x107 ft 
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(45.0-0.7)(1077) 


rp 
| 


= +4.3x107 ft 


ayo (42.4-1.8)(1077) 


Salo 8s 


h. = (-10.8-3.6)(107) 


eo et 


he = (+2.1-1.8)(1077) 


Bees) ass 


The resultant displacements of each HEROES Ey may be obtained as 
the square root of the sum of the squares of the net horizontal and vertical 


displacement components at that joint. 
= oe = 2 
pena NAN o> ake Ah 8.4.38 
j= va,” + G5) (8.4.38) 


Soreane Be 


= 


Ap 


* 


8.4.5 Calculation of Natural Periods of Vibration 


th.4x107 ft 


4. 3x10 7? ae 


As an illustration of the calculation of the natural periods of 
vibration of the system, the natural period of vibration is calculated below 
for the structure as it approaches its first maximum response at the leeward 
haunch. Because the "stiffness" of the soil varies with the overpressure 
restraining the soil wedges (cf. Section 4.4), the time selected for this 
example is the time corresponding to maximum soil "stiffness" up to this 


instant. 





IL 


From Figure D.4 one determines that the motion of the structure 
is “Symmetric-Out" (cf. Figure 5.3) as the system approaches its first maxi- 
mum displacement at the leeward haunch and that this maximum occurs at 
t = 0.0173 seconds. Thus, the haunches are displaced outward and a "passive" 
state of stress is developing in the soil on both sides of the arch. From 


Table 8.6, or from plots of the overpressure-time functions, the maximum 





soil restraining pressures acting on the "passive" soil failure wedges during 
the interval between t = O and t = 0.0173 seconds are found to occur at 
t = 0.0126 seconds when Ply = -89.0 psi (obtained from plot) and Pos = -96.4 
psi. 

The soil stiffnesses corresponding to these restraining pressures 


are obtained from Equations (8.4.15). 


Ky /8 = 618.9 - 178.2 Pe 

K!,/a = 618.9 - 178.2(-89.0) 
= 16,480 k/ft 

Ki,/a = 618.9 - 178.2(-96.4) 


17,800 k/ft 


The arch stiffness parameters required for substitution into the 


frequency equation may be found by referring to Equation (7.2.3). 


t. = — (5+a/c+b/d+2c/atbc/ad) 


1 
OL 
2 
Bay ee E A 2-695 , 2a 4bO 2(2.4h0) . (2.440) 
2(6.378)° a 5-893 " 5.895” (5. g93)2 


121.2 k/ft 
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KB 


p 2 (3 + a/c - b/d - be/ad) 


7 
i 





abe 893 2.40  (2.440)* 
= 212 * o.4ho ~ 5.803 7 2 
2(6.378) (5.893) 


= 66.3 k/ft 


The equivalent mass parameters for the '‘Symmetric-Out" motion are 


those appearing in Equation (8.4.18). 


UL = 0.4916 age" / PE 


1,) 54) 


it 


»: Oe Gr1GsW eeccee er 


1,3. “3,3 


The parameters obtained above may now be substituted into Equation 


(7.2.9) and the resulting equation solved for we 


} h ip 
haa | wo fy al taKis/e) 
i Hy mit *Key/2) - 5 (4 ats m) | ue 

r a A ; 
7 | (+ HRS /0) (14K /8)-(75) 
= 0 


((0.4916)°-(0.1034)*] wt - [0.4916(121+17, 800) 
+ 0.4916(121416, 480 )-66.3(0.1034+0.1034) } we 
+ [ (121416, 480) (121417, 800) -(66.3)*] 


OMe le, wt -16 , 960 we +279. 4x10° = 0 


Bee 22 29 ,000 rad”/sec* 
Dai 
we = 4h, LOO peo eee 





15) 


The value of w” corresponding to "Symmetric-Out" motion may be found by 


nm Bm 


and w’ |, in turn, into either of Equations (7.2.8). 


substituting 3,0? 
3 


Ww 
5,1 


- 2 vw eZ a 
C2 Om 1) -Ki,/a)h, + cua wt tg dh, = 0 


For ig 


al = 29,000 mele sce 
3 


[0.4916(29, 000 )-121-16, 460 }h, +[ -0.1034(29, 000 )+66.3 Jb, = 0 


h, = -0.50 h, 


since hy and ae are of opposite signs, the haunches are displaced in opposite 
directions and the configuration is "symmetric". (Substitution of ce > = 
P 


44, 4OO raa/sec” into Equations (7.2.8) gives h, = +1.16 hi» indicating that 


2 
both haunches are displaced in the same direction and that the configuration 
is therefore "asymmetric"). 

The correct value of we is thus 29,000 rad“/sec” from which 
w = 170.3 rad/sec. This value may now be substituted into Equation (7,2.16% 
in order to find the natural period of vibration of the system at this 


particular time in this particular motion. 


Tn = ax/w 
T, = on/170.3 


0.047 seconds 
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TABLE 8.1 


SHOCK WAVE PROPERTIES 











Peak Side-on Overpressure psi 1000 100 50 
Distance from 1 MI Contact yds 480 1100 1500 
| Surface Burst 
Duration of Positive Phase sec 1.48 1.48 1.60 
Velocity of Shock Front in Air fps 8700 2950 220 
TABLE 8.2 


OVERPRESSURE-TIME COORDINATES FOR 50-psi SHOCK WAVE 


ee sec O 0.065 0.20 0.40 0.80 
ect) psi 50 LO 26.5 15.5 5.0 
TABLE 8.3 


OVERPRESSURE-TIME COORDINATES FOR 100-psi SHOCK WAVE 


t sec O 0.070 0.205 0.395 0.702 
pit ) psi 100 60 31.5 16 6 
TABLE 8.4 


OVERPRESSURE-TIME COORDINATES FOR 1000-psi SHOCK WAVE 


om sec O 0.01 0.04 0.12 0.40 
Tt!) psi 1000 600 250 100 25 
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TABLE 8.5 


SUMMARY OF NUMERICAL EXAMPLES 


Depth of Reinforced Concrete Arch Corrugated Metal Arch 
Cover Over —- ee See ee ee 
Crown 100 psi 1000 psi 50 psi 100 psi 1000 psi 

0 25; (1) (1) (1) (1) 

= (Span) (1) (1) (1) (1) - 

; (Span) (1),(3) (a) GQ) Cie) Sai 

= (Span) (1) (1) (1) (1) : 


Notes: (1) K > i/h, C, = 3000 fps 
2) a Wy) dy) Syn 3) teed 
(3) C. = 1000, 2400, 3000 fps 








TABLE 8.6 


OVERPRESSURE-TIME RELATIONSHIPS 
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oading Region Parameter Coordinates 
la-lb t=t'+0.0007 sec 0.0014 “0.0707 -0n2057 On5g57 
t sec - 0.0007 0.0700 0.2050 0.4950 
-P, psi O 99.4 60 Sl 16 
2a-eb t=t'+0.0024 sec 0.0014 0.0034 0.0724 0.2074 0.3974 
ce sec - -- 0.0010 0.0700 0.2050 0.3950 
“Py psi 0 99.2 60 31.5 16 
3a-3b t=t'+0.0041 sec 0.0034 0.0048 0.0741 0.2091 0.3991 
t sec - -- 0.0007 0.0700 0.2050 0.3950 
-P,, psi O O 99.4 60 31.5 16 
2a-2 t=t'+0.0019 sec O 0.0014 0.0024 0.0719 0.2069 0.3969 
te sec - -- 0.0005 0.0700 0.2050 0.3950 
“Pr, psi O O 99.6 60 41.5 16 
2-2b t=t'+0.0029 sec O 0.0024 0.0034 0.0729 0.2079 0.3979 
va sec - -- 0.0005 0.0700 0.2050 0.3950 
“Poy psi O 99.6 60 31.5 16 
Oa-Ob t=t'-0.0008 sec 0 0.0013 0.0692 0.2042 0.392 
t' sec 0.0008 0.0021 0.0700 0.2050 0.3950 
“Py psi 67.9 98.5 60 Gules 16 
04-Ob t=t'-0.0032 sec O 0.0013 0.0668 0.2018 0.3918 
+' sec - 0.0045 0.0700 0.2050 0.3950 
oe psi 82.6 96.4 60 Ses 16 
ha-lb t=t'+0.0056 sec O 0.0045 0.0077 0.0756 0.2106 0.4006 
Oe sec - -- 0.0021 0.0700 0.2050 0.5950 
“Pes psi 0 0 98.4 60 nba. 16 
ha-kp' t=t'+0.0080 sec 0.0035 0.0126 0.0780 0.2130 0.4050 
ie sec - -- 0.0046 0.0700 0.2050 0.3950 
“Pos psi O 96.4 60 5y.5 16 
aa-ab t=t'4+0.0024 sec 0.0052 0.0724 0.2074 0.3974 
e sec - 0.0028 0.0700 0..050 0.3950 
-p psi 4 97.7 60 Ae 16 


a 








aT 


PP Lene 


[2-500 





Figure 6.1. Reinforced Concrete Arch, qd. = B/k 


My = 3.50 ft-kips —— —— — 


Kej- [tl7.2 ft-kips per radian 


Wy =0.003133 radians Yj 


Figure &.2. Moment-Rotation Diagram for Joints of Model 










\2° 


ae 


L A= 0.24 in q" 


. \ 





Figure 5.3. Arch Rib Cross Section 
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Figure 8.4. Equivalent Soil Resistance-Displacement Diagram 








ae 


CHAPTER IX 


SUMMARY 


The results of the numerical investigation performed in connection 


with this thesis are summarized in this chapter. Included also are a listing 





of the conclusions which have been drawn on the basis of this investigation 
and some recommendations for further study. The numerical examples are des- 
cribed in Chapter VIII and the results presented in tabular and in graphical 
form in Appendices C and D. Im general, the results give an informative 
picture of the lateral response of shallow-buried arches to nuclear blast 


loadings under a variety of conditions. 


9.1 Discussion 

In interpreting the results of the numerical investigation as pre- 
sented in the appendices, one should keep in mind that the relative values of 
the response are more significant than their absolute values. The relative 
values permit comparison of arch behavior as geometry and loading parameters 
are varied. The absolute values are somewhat misleading because of their 
dependence upon the assumptions made in the dynamic analysis. In general, 
the calculated displacements are smaller than the displacements which one 
would expect in an actual structure. This is true first because of the re- 
placement of the actual structure having continuous flexibility with a dynamic 
model having only two degrees of freedom, and second because of the idealized 
soil conditions assumed to exist before and during response. For example, 
any lack of compaction in the backfill surrounding an actual arch will be 
manifested directly as increased deformation of the arch because of the 


relatively large influence of the resistance mobilized in the soil on the 
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response of the arch. This is discussed in detail below. Further, 
displacements of an actual structure will be affected by differences be- 
tween the actual blast wave characteristics and those of the idealized 
blast wave used in these calculations. This effect could lead to either 
greater or smaller values of displacement. 

Another factor influencing the numerical results is the absence 
of any damping factor in the equations of motion. While damping is not con- 
sidered to have a significant effect upon the magnitude of the lateral re- 
sponse during the initial portion of the response history, it will reduce 
later maxima in all cases short of failure. For this reason, the comparisons 
made in Appendix C are based on "initial", or first, maximum values at each 
joint. 

A comparison of the natural periods of vibration of the structures 
is presented in Tables C.1 and C.2. Im general, as the depth of cover is 
increased, the natural period of vibration of the structure is decreased. 
This indicates that the "stiffness" of the soil tends to be more significant 
than the additional mass which the soil adds to the vibrating system. At 
first glance the natural periods corresponding to zero depth of cover at the 
crown appear to be out of sequence. This is because the change in mass from 
the “no cover" to the “d, = O" condition is much greater than the changes in 
mass between any two of the succeeding cover geometries listed. In fact, 
because of the limitations imposed on the dimensions of the soil masses by 
Equations (5.1.3) and (5.1.4), the change in mass with increasing depth of 
cover eventually becomes zero. On the other hand, the stiffness of the 
structure continues to increase with depth becoming, relatively, quite large 


for the case of a soil depth at the crown equal to one-half of the arch span. 
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It should be pointed out here that because of the nature of the 
assumptions upon which the dynamic model is based, the natural periods of 
vibration given in Tables C.1 and C.2 for the "Asymmetric’ modes are closer 
to the periods of the actual system than the periods for the "Symmetric" 
modes. The symmetric modes are influenced to a significant extent by the 
extensional as well as the flexural characteristics of the arch while the 
asymmetric modes are almost entirely flexural, or "bending", modes. 

Table C.2 contains the natural periods of the systems as they 
approach their first maximum response at the leeward haunch. Aside from a 
relatively small initial inward movement of the windward haunch and upward 
movement of the crown in certain cases, this is the first maximum reached 
by the system (cf. Appendix D). In all of the cases studied the system 
motion is 'Symmetric-Out” (cf. Figure 5.3) as the system approaches this 
maximum. Since the natural period is a function, in part, of the soil stiff- 
ness (cf. Equations (7.c.9) and (7.2.10)), and the soil stiffness in turn is 
a function, in part, of the overpressure restraining the soil failure wedges 
(cf. Equations (4.4.10) and (4.4.12)), it is necessary to specify the par- 
ticular time at which the natural period is to be calculated so that the over- 
pressure acting on the soil failure wedges may be found. The time selected 
in each case is the time corresponding to maximum restraint of the "passive’’ 
soil failure wedges as the system approaches its first maximum response at 
the leeward haunch. 

In Table C.2 it is interesting to note that, as in Table C.1, the 
natural period tends to decrease with increasing depth of cover. Also, the 
natural period decreases with increasing overpressure. This results from 
the fact that the overpressure acts to increase the soil stiffness but has 


no effect upon the mass of the system. 
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In Tables C.3, C.4, and C.5 the first maximum values of response 
are compared. The tendency of the arch to move initially inward at the 
windward haunch and upward at the crown in certain cases (cf. Appendix D) has 
not been considered in selecting "first" maxima. As can be seen from the 
response curves in Appendix D, the sequence of response is, in general: 

1. The response reaches a maximum (outward) value at the leeward 
haunch. 

2. The response reaches a maximum (downward ) value at the crown. 

3. The response reaches a maximum (outward) value at the windward 
haunch. 

In all cases as the depth of earth cover is increased, the magnitude 
of response is decreased. In general, as the overpressure is increased, the 
magnitude of response is increased. However, this increase in the magnitude 
of response is not as great as might be expected. This is caused by the 
increase in soil resistance with increasing overpressure level. An exception 
to this occurs in the case of the corrugated metal arch with no soil cover at 
the crown. Here an increase in the overpressure loading from 50 psi to 100 psi 
results in a reduction in the magnitude of the response. This is due to the 
fact that the mass remains constant with a change in overpressure while both 
the loading and resistance parameters vary with overpressure level. Increasing 
the overpressure level to 1000 psi overcomes the "mass effect" and the magni- 
tude of the response increases as expected. 

It should be noted here that comparing first maximum in those cases 
having no soil cover at the crown is somewhat misleading. Examination of the 
response curves in Appendix D reveals that the second maxima tend to become 
very large for certain combinations of structure geometry and loading (cf. 
Figures D.7, D.15 and D.22). For these cases, the overpressure is not suf- 


ficient to hold the structure down and the crown "pops up." 
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Another point of interest may be illustrated using the data 
presented in Figures C.3, C.4 and C.5. This concerns the relative contribu- 
tions of the arch and the soil to the lateral deformation resistance of the 
system. Consider, for example, the reinforced concrete arch with a depth of 
soil cover at the crown equal to one-fourth of the arch span and subjected 
to a 100-psi shock wave. Now consider the configuration of the structure 
at the time of its first maximum response at the leeward haunch, i.e. at 
t = 0.0175 seconds. From Table C.5 the horizontal displacement of the leeward 
haunch, hs at this time is +0.0054 feet. From Figure D.3, the horizontal 
displacement of the windward haunch, h,, at this time is ~-0.0048 feet. 
Proceeding as in the sample calculations of Section 8.4 gives corresponding 
values of the equivalent arch and soil resistances as; 

“a 
8s S 


This illustrates the fact that the arch contribution to the lateral response 


-0O.87 kips 


& 
iH 


5] -77.09 kips 


+0.91 kips +93.47 kips 


&) 
i 


of the system is relatively small as compared with the soil contribution. 
This also indicates that the larger displacements obtained in the case of the 
corrugated metal arch may be due more to the larger span of the corrugated 
metal arch rather than to the smaller resistance of the steel arch rib as 
compared with the concrete arch rib. 

The data contained in Table C.6 show the effect of varying the 
ratio, oa of the horizontal to vertical overpressures induced in the soil 
by the shock wave in air (cf. Figures D.17 through D.20). In general, as 
Por is increased, the magnitude of the response is decreased. Exceptions to 
this are the initial inward movement of the windward haunch and the accompany- 


ing upward movement of the crown. These movements are magnified by increasing 


K 
pr 
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The effect of the angle of inclination of the shock front in the 





soil was studied and the results of the study presented in Table C.7 (cf. 
Figures D.3, D.4 and D.5). Since the angle of inclination of the shock 
front in soil, 1, is a function of the relative values of the velocity of 
the shock wave in air, U, and the seismic velocity of the soil, C3 (cf. 
Appendix B), varying ) for a constant U has the same effect as a correspond- 
ing variation in C. (cf. Equations (B.3)). Table C.7 lists the maxima for 
the solutions in chronological order. In general, increasing the angle of 
inclination of the shock front in soil with respect to the horizontal 
accentuates the initial inward movement of the windward haunch (and accompany - 
ing upward movement of the crown) but has no significant effect on the magni- 
tudes of the other maxima. 

Response curves for all of the cases studied are presented in 
Appendix D, Figures D.1 through D.23. The figures are arranged with those 
for the reinforced concrete arch first followed by those for the corrugated 
metal arch. The figures are grouped by overpressure level and are ordered 
within each group by increasing depth of cover. The solutions for Pa 50 
psi and Bee = 100 psi are plotted for the first 90 milliseconds after the 
shock wave reaches the structure, those for Deas 1000 psi for 18 milliseconds. 
The curves include the first maximum response at each joint in all cases. 

Configuration diagrams, Figures D.24 through D.30, are presented 
for two of the solutions. These diagrams are arranged chronologically and 
illustrate the configuration of the structure at each of the first maxima 
(neglecting the small initial movements of the windward haunch and the crown 
as before). In addition, for purposes of comparison, Figure D.30 shows the 
configuration of the corrugated metal arch at the time of the second maximum 


response at the windward haunch. 
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The configuration diagrams are drawn first neglecting and then 
including rib shortening. The effect of rib shortening is seen to be a 
general inward and downward movement of the structure. While rib shortening 
INight seem to cause the arch to move away from the soil and thus invalidate 
some of the assumptions upon which the lateral response analysis is based, 
test results have indicatea(27? that the soil will deform also and that the 
deformations of the soil are of the same order as those of the arch. Further, 
the blast loading is transmitted to the arch by the soil so that any movement 
of the arch away from the soil tends to decrease the load on the arch, thus 
counteracting the tendency for the arch and the soil to move apart. Finally, 
the blast wave in air tends to push the soil down around the arch by the 


very nature of its action. 


9.2 Conclusions 

The numerical results discussed in the preceding section lead to 
the general conclusion that well compacted soil surrounding a shallow-buried, 
underground arch contributes the major portion of the resistance of the arch 
to lateral deformation when subjected to the blast load resulting from a 
nuclear explosion. 

If consideration is limited to those cases where the depth of 
cover is sufficient to prevent the crown from moving upward out of the ground, 
e.g. d, > B/8, and the arch rib cross section is sufficient to withstand the 
thrust developed in the arch without general yielding or buckling, it may be 
concluded that the arch will successfully withstand the applied loading. 
Further, if the following parameters are varied individually in turn, it may 
be concluded that: 

1. As the depth of cover is increased, the magnitude of the 


response will decrease. 
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2. As the overpressure level is increased, the magnitude of the 
response will increase. 

4. As the ratio of horizontal to vertical overpressure induced 
in the soil by the shock wave in air is increased, the magnitudes of the 
initial movements of the windward haunch inward and the crown upward will 
increase; however, the later magnitudes of the lateral response will decrease. 

4, As the angle of inclination with respect to the horizontal of 
the shock front in soil is increased, the magnitudes of the initial move- 
ments of the windward haunch inward and the crown upward will increase; 
however, the angle of inclination of the shock front in soil has no signifi- 
cant effect upon the magnitude of the later response. 

5. As the span of the arch is increased, the magnitude of the 
response will increase. 

6. The maximum response of the arch will occur at the crown, 
principally in the vertical direction. 

In general, it may be concluded that the natural periods of 
vibration of the bending modes of the arch will be decreased by: 

1. Increasing the depth of soil cover, and 


2. Increasing the overpressure level. 


9.3 Recommendations 

While the present study enables one to make certain predictions 
regarding the relative behavior of shallow-buried arches under various 
combinations of geometry and nuclear blast loading conditions, it is limited 
by the assumptions which are made in writing and solving the equations of 
motion for the system. This study might logically be extended in one or 


more of the following areas. 
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There appears to be a significant lack of knowledge regarding the 
interaction between the structure and the soil. While classical earth pres- 
sure theory gives an indication of the resistance offered by the soil, this 
aspect of the problem could well be investigated in detail. For example, 
there is a notable lack of experimental data relating to the amount of de- 
formation required to develop the passive state of plastic equilibrium. 

A soil resistance-deformation relationship somewhat more in agree- 
ment with existing experimental data might be obtained by considering a more 
rapid shift in the direction of the deformation resistance offered by the 
soil as the system "rebounds." This approach is unattractive from a computa- 
tional standpoint but may bear further investigation. 

A logical extension of this study is a consideration of more than 
two degrees of freedom for the dynamic system. This would give a closer 
approximation to the actual structure and would permit simultaneous considera- 
tion of extensional and bending modes of response. It might also permit 
footing motions to be considered. 

The determination of the Overpressure causing elastic buckling of 
the arch rib is an interesting problem. While the calculations in Chapter 
VIII indicate that this condition is not critical, they are approximate. A 
closer estimate of the critical overpressure would be helpful. 

A closer approximation to actual response values might be obtained 
by considering the rise time of the pressure pulse in the soil. This is dif- 
ficult to define at depths of interest in this problem using existing 
information. In fact, the entire problem of the transmission of the shock 
wave through the soil merits careful study. 

A closer approximation to the later maxima of the displacement-time 
relationships may be obtained by including a damping term in the equations 


Of motion. Since this term would appear to be most significant in those 


cases with little or no soil cover, it should be kept in mind that depths 


9f soil cover approaching zero at the crown may be unattractive from a 
practical standpoint. This is due to the fact that blast resistance is 
nly part of the over-all problem of providing protection against the effects 
of nuclear weapons. Thermal and nuclear radiation effects will also be 


significant at overpressure levels of interest in this problem. 
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APPENDIX A 


HYDROSTATIC PRESSURE RATIOS IN SOIL 


SUBJECTED TO NUCLEAR BLAST LOADING 


SSS 2 Te a a Se SR Se, oe 


This study is an investigation of the effect of blast overpressure 
on calculated values of the hydrostatic pressure ratios in soil corresponding 
to the "Active State of Plastic Equilibrium" and the "Passive State of 
Plastic Equilibrium." The calculations are based upon classical earth pres- 
sure theory using the Coulomb assumptions to describe conditions in the soil 
failure wedges (cf. Section 2.2). A perfectly homogeneous, cohesionless soil 
is assumed with the blast overpressure considered as a surcharge loading. 

Two solutions for each of the hydrostatic pressure ratios are 
compared, one obtained using the "Trial Wedge" method and the other obtained 
using an algebraic expression of Coulomb's method. The example chosen for 
these calculations is a semicircular arch, 16' 8" in diameter, with a depth 
of cover over the crown equal to one-fourth of the arch span. This structure 
is illustrated in Figure A.1. The blast loading is taken as the 100-psi 
shock wave resulting from a one-megaton nuclear explosion at the ground 
surface. The pressure surfaces are defined based on the dynamic model as 
described in Chapter IV. The soil is assumed to have a unit weight of 120 
pounds per cubic foot and an angle of internal friction of 30 degrees. The 
effective friction angle between the arch and the soil is taken as two-thirds 
of the angle of internal friction of the soil, or 20 degrees. 

The geometry and soil conditions of this problem are the same as 
those of the illustrative example discussed in Section 8.4. Thus, the 
average overpressure acting on the failure wedges may be obtained from 


Table 8.6. For the purposes of this study, the maximum value of the average 
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overpressure restraining the "passive’ soil failure wedge is used, i.e. 
Ps, = -96.4 psi. While the average overpressure restraining the “active” 
soil failure wedge has a somewhat higher maximum value, i.e. Py = -98.3 psi, 
it is considered at a later time when it has decayed to a value of %.4 psi 
so that both sets of calculations may be made for the same restraining 
pressure. This is slightly more convenient. Actually, since this investi- 
gation deals principally with the effect of "large" surcharge loadings on 
calculated values of the hydrostatic pressure ratios corresponding to plastic 
equilibrium in soil, any reasonable value of the restraining pressure may be 
used. 

Coulomb's formula for the hydrostatic pressure ratio corresponding 
to the "Active State of Plastic Equilibrium," K,, as derived in reference (14), 


A 
is given below as Equation (A.1). 


vay 
K. = cos (P-w (Am 
‘i 2 lsin(5+? )sin(P-i ie 
cos? weos(o4u) [2 +) Satnbeoatond 
where ? = angle of internal friction of the soil 


W = angle of inclination with respect to the vertical of the 
lower bars of the model 


8 = effective friction angle between the arch and the soil 


i = slope of the ground surface with respect to the horizontal. 


Substituting the soil properties listed above and the geometry shown in 


Figure A.1 into Equation (A.1) gives K, = 0.510. 


a cos”(30° -22°29.6') 


0,zn0 O- 2 
cos(22°29.6' )eos(20°+22°29.6") P + i 


cos(20°+22°29.6' )cos(22°29.6') 


0.510 
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Coulomb's formula for the hydrostatic pressure ratio corresponding 
to the "Passive State of Plastic Ejuilibriuwn, " K,, as derived in reference 


(14), is given below as Equation (A.2). 


. cos (Pew) (A. 2) 


2 i ~}sin(?+5 )sin Pia) | 
cos” weos(w-6 ) 2 [2 RES Gh 





Substituting values into Equation (A.2) gives Kp = 3.433. 


7 cos(30°+22°29.6') 
PP = ——4 =" 
fe) O 
cos(22°29.6' )eos(22°29.6'-20°) E - | —enasPeaeP esa 2F ——] 
cos(22°29.6'-20° )eos(22°29.6') 


5.433 


The results of the graphical solutions for the resultant forces on 
the pressure surface using the “Trial Wedge" method are given in Table A.1. 
The corresponding hydrostatic pressure ratios are calculated below. 


No overpressure acting: 





A 
mA" Ty a 
2 °s oO 
where P, = resultant force acting on the pressure surface for an "active: 


state of stress in the soil 


4 800 


5 (120)(12.500)* 


x = 42,800 
A ~ 9,375 


Ono ic 
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m= 
a: 2 
e ye 
where Po = resultant force acting on the pressure surface for a 
"passive state of stress in the soil. 
» 91,050 
9,375 
= 3.312 
With overpressure acting: 
K, = ‘A 
Al oo 
2 Ue as Ps 5% 
7 2, 500 
9,375+96.4(12. 500)( 1 
= 0.506 
ee 
= T as 
2 75% ? ate 
7 619, 000 
~ 182,900 
= 04 


The inclination of the failure surfaces with respect to the 
horizontal are presented in Table A.e. The inclinations obtained from the 
“Trial Wedge" solutions are compared with the assumed inclinations of (E+ 5) 
for the ‘active’ failure surface and (F - ) for the "passive failure 
surface for the numerical work in this thesis. 

From the comparative results presented in Table A.1, it may be 
concluded that the overpressure loadings have a negligible effect upon the 
hydrostatic pressure ratios in the soil. Further, both the "Trial Wedge" 


method and Equations (A.1) and (A.2) give essentially the same values of 





K, and K.- Finally, in Table A.2 it may be seen that the assumed inclina- 


sions of the failure surfaces in the soil are in reasonable agreement with 


the inclinations obtained from the "Trial Wedge” solutions. 
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TABLE A.1 


HYDROSTATIC PRESSURE RATIOS IN SOIL AS COMPUTED 
BY "TRIAL WEDGE” METHOD AND COULOMB FORMULAS 


State of Stress rial Wedge Formula 
in Soil Without Blast With Blast 
Resultant Force on Active, P, 4,80 92.5 
Pressure Surface, 
kips Passive, Pp Ss 619 
Hydrostatic Active, K, 0.512 0.506 0.510 
Pressure Ratio 
Passive, K, Sane Siler 5. 4ae 
TABLE A.2 


INCLINATION WITH RESPECT TO HORIZONTAL OF 
FATLURE SURFACES IN SOIL 


— Trial Wedge Solution ASeumea Vala 
Without Blast With Blast 
Pp 
Active 64° 64° (7 + a eso 


Passive 26° 26° (F - =) = 30 












Figure A.1. Semicircular Arch, B = 16'S, d= F/4 








179 





APPENDIX B 
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The object of this study is the determination of the configuration 
of the shock front of the air-induced shock wave in soil. Primary interest 
is limited to those depths below the ground surface associated with the 
shallow-buried arches studied in this thesis. Knowledge of the shock front 
configuration is needed in order to define the loading on the arch during 
its envelopment by the shock wave (cf. Chapter VI). 

The approach used in this study is outlined below. Essentially, 

a point of interest is selected and the location of the shock front of the 
air-induced shock wave in soil determined for the time after the explosion 
corresponding to the time of arrival of the shock wave in air at this point 
of interest. It is assumed that the shock wave in air travels with velocity, 
U, which is acontinuously varying function of the peak overpressure and 


ambient atmospheric conditions as expressed by Equation (p.1), (28) 


U=c V1 + 6p /TP_ (B.1) 


O Pso Oo 


where eae anbient sound velocity in air ahead of the shock front 
(assumed as 1130 feet per second) 
Pa peak air overpressure at the shock front 
a = ambient atmospheric pressure ahead of the shock front 


(assumed as 14.7 psi). 

The shock wave in soil is assumed to travel at the seismic velocity of the 
soil, Os This study considers a semi-infinite, homogeneous, soil mass with 
no surface irregularities and with a constant seismic velocity. Further, 

the seismic velocity of the soil is assumed not to vary with depth, direction 


or with the overpressure level in the shock wave. 
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The shock wave tends to travel symmetrically outward from Ground 
Zero as a roughly spherical segment of ever-increasing radius. Therefore, 
this problem may be considered in two dimensions and the shock front assumed 
to extend indefinitely in the third dimension. At distances from Ground 
Zero which are of interest in this problem, the radius of curvature of the 
shock front in air is quite large, so that this assumption does not lead to 
Significant error. 

The calculations in this appendix are based upon a one-megaton, 
contact surface burst. The pressure-distance relationship is taken from 
reference (17) The point of interest is selected at a distance from Ground 
Zero such that the peak value of the overpressure at the shock front is 
100 psi. Seismic velocities of 1,000, 2,000, 3,000 and 6,000 feet per 
second are considered. 

The procedure followed in the calculations may be divided into 
six steps as follows: 

1. Select the point of interest. 

©. Select a series of check points, 1, 2, 3,...n, between the 
point of interest and Ground Zero. Plot these points to scale along a hori- 
zontal line representing the ground surface. 

4. Calculate the time interval, At, required for the shock wave 
in air to travel between each of the check points and the point of interest. 

4, Calculate the distance traveled by the shock front in soil, 
di at the seismic velocity of the soil, Cos during each time interval 
found in Step 3. 


qd. * (vt ee (B.2) 


5. Using each check point, n, selected in Step 2 as a center, 
swing a semicircular arc in the "soil" with a radius equal to the value of 


a calculated for that check point in Step 4. 
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6. A curve tangent to the arcs drawn in Step 5 is thus an 





approximation of the shock front in soil at the instant the shock front 
reaches the point of interest selected in Step 1. 

The results of this study are presented in Figures B.1 through 
B.4. 

While this study is quite limited, it may be concluded tentatively 
that for the relatively shallow depths considered in this thesis, the shock 
front in soil may be considered as a plane perpendicular to the direction of 
travel of the shock wave and inclined at an angle, 1, with respect to the 
horizontal as given by Equations (B.3). These equations may be written with 


the aid of Figure B.5. 


C/U < I: | = arcsin c./U 
(B.3) 
C/U & 1: |= n/2 
where ee = seismic velocity of the soil 
U = velocity of the shock wave in air at the point of interest. 


While the shock front in soil at depths close to the ground surface for 
c_/uU > 1 remains essentially vertical, the fact that it tends to "run ahead" 
of the shock front in air may tend to increase the "rise time" of the 
overpressure-time function in the soil at a particular point of interest. 
Since an increase in the rise time tends to reduce the dynamic response of 
a system to a given impulse, this effect is not considered in the numerical 


work of this thesis (cf. Section 2.4). 
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APPENDIX C 


TABULAR COMPARISONS OF NUMERICAL RESULTS 


This appendix contains a series of tables in which comparisons 
are made of the various numerical solutions. In general, quantities are 
expressed as follows: 

Displacements, an in feet 

Overpressures, p, in psi 

Cover depth at crown, dw in proportion of arch span, B 
Time, t, in seconds 


Natural periods of vibration, TW in seconds. 


The sign convention used in these tables is the same as that stated 
in Section 1.2, i.e. displacements to the right and upward are considered 
positive. The shock wave is considered to enguif the structure from left 
to right. 

A summary of the numerical examples considered is contained in 


Table 8.5. 
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TABLE C.1 


NATURAL PERIODS OF VIBRATION WITH NO OVERPRESSURE 
ACTING ON SOIL FAILURE WEDGES 


(seconds ) 





Reinforced Concrete Arch Corrugated Metal Arch 


Soil Cover 
Asymmetric Symmetric Asymmetric Symmetric 
In Out In Out 
None O.147 Omsoew O- 115 0.484 0.589 0.389 
d, = 0 0.163 0.156 0.133 0.175 0.234 0.189 
= B/8 0.144 0.179 0.153 0.18% 0.255 0.214 
= B/& 0.129 Oolgem nO elue 0.162 0.233 0.199 
= B/o 0.092 0.136 0.123 0.111 0.173 0.154 





TABLE C.2 


NATURAL PERIODS OF VIBRATION AT FIRST MAXIMUM 
RESPONSE OF LEEWARD HAUNCH 
(seconds ) 


(Motion: "Symmetric-Out" in all cases) 


Depth of Reinforced Concrete Arch Corrugated Metal Arch 
Cover Over ee —— 
Crown 100 psi 1000 psi 50 psi 100 psi 1000 psi 


0 0.034 0.012 0.060 0.04% 0.015 
B/8 0.038 0.013 0.073 0.045 (1) 
B/4 0.037 0.015 0.073 0.054 0.019 


B/2 0.032 0.011 0.062 0.06 (1) 


(1) Not calculated. 
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TABLE C.3 


FIRST MAXIMUM RESPONSE AT WINDWARD HAUNCH 
(Neglect small initial inward movement) 











Cover Para- Reinforced Concrete Arch Corrugated Metal Arch 
Depth meter‘1) (Span = 16'6") (Span = 25' 1 1/4") 
at See 
Crown 100 psi 1000 psi 50 psi 100 psi 1000 psi 
O a -0.0135 -0.0137 -0,0191 -0.0138 -0.0199 
10 0.0214 0.0068 0.0383 0.0254 0.0087 
B/8 hy -0.0087 -0.0092 -0,0121 -0.0132 (2) 
12 0.0236 0.0075 0.0450 0.0336 
B/ 4 h, -0.0063 -0.0068 -0.0084 -0.0093 -0.0101 
c 0.0229 0.0074 0.0450 0.0337 0.0110 
B/2 ay -0. 0040 -0.0043 -0.0051 -0.0057 (2) 
t 0.0204 0.0065 0.0393 0.0296 
(1) h. = horizontal displacement of windward haunch, in feet. 


- 
t = time corresponding to h 


(2) Not computed. 


ie in seconds. 
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TABLE C.4 


FIRST MAXIMUM RESPONSE AT CROWN 
(Neglect small initial upward movement) 





Cover Para- Reinforced Concrete Arch Corrugated Metal Arch 
Depth (1) (Span = 16'8") (Span = 25' 1 1/4") 
x meter 
Crown 100 psi 1000 psi 50 psi 100 psi 1000 psi 
O he ~0.0237 -0.0281 -0.0344 -0.0275 -0,0441 
t 0.0186 0.0053 0.0325 0.0213 0.0069 
B/8 be ~0.0151 -0.0173 -0.0205 -0.0218 (2) 
t 0.0205 0.0060 0.0390 0.0285 
B/h h, ~0.0108 ~0.0123 -0.0141 -0.0151 -0.0171 
t 0.0198 0.0060 0.0388 0.0285 0.0085 
B/2 h. -0.0069 ~-0.0073 -0.0085 -0.0095 (2) 
ie 0.0174 0.0053 0.0339 0.0252 


(1) “ts = vertical displacement of crown, in feet. 


t = time corresponding to h._, in seconds 


oF 


(2) Not computed. 
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TABLE C.5 


FIRST MAXIMUM RESPONSE AT LEEWARD HAUNCH 





\Cover Para- Reinforced Concrete Arch Corrugated Metal Arch 
= PQ = t i 
mee. 1) (Span = 16'8") (Span = 25' 1 1/4") 
100 psi 1000 psi 50 psi 100 psi 1000 psi 

h, +0.0119 +0.0163 +0.0181 +0.0162 +0.0268 
t 0.0165 0.0050 0.0293 0.0195 0.0065 
h, +0.0075 +0.0094 +0.0102 +0.0108 (2) 
t 0.0181 0.0055 0.0341 0.0247 
hs 40.0054 +0. 0066 +0.0070 +0.0076 +0. 0092 
t 020175 0.0053 0.0337 0.02h6 0.0075 
h, +0.,0037 +0.0039 +0.0044 +0.0049 ‘e) 
t 0.0154 0.0045 0.0295 0.0220 


horizontal displacement of leeward haunch, in feet. 
= time corresponding to h 
(2) Not computed. 





in seconds. 


Bi 
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TABLE C.6 


EFFECT OF HORIZONTAL-TO-VERTICAL OVERPRESSURE 
RATIO ON LATERAL RESPONSE** 


(First maximum values) 











K 
Be oter (2) 

i/4 AE: 3/4 1 

Po) ax 0.0007 0.0014 020022 0.0034 
Time 0.0068 0.0086 0.0104 0.0727 

(+h, ) naxx 0.0076 0.0057 0.0043 0.0054 
Time 0.02h6 ©7022) 0.0192 0.0164 

(be) ave ~O.0151 -0.0106 -0.0060 -0.0016 
Time 0.0285 0.0279 0.0268 0.0238 

Be as -0.0093 -0.0076 -0.0059 -0,0044 
Time 020737) 0.04540 0.0542 0.0402 


(1) Based on Corrugated Metal Arch, p_. = 100 psi, d, = B/4. 


(2) h, = horizontal movement of windward haunch, in feet 
= = horizontal movement of leeward haunch, in feet 
ag = vertical movement of crown, in feet 


Time, in seconds 





Inclination of Shock 


seismic Velocity in Soil 


(1) 
(2) 


(3) 
(4) 


EFFECT OF INCLINATION OF SHOCK FRONT IN SOIL 
ON LATERAL RESPONSE 


Front in Soil 


(+h. ) ax 
Time 
(+hy ) ax 

Time 


(+h, ) ax 


Time 


Genes 
Time 


( aE) dig 
Time 


TABLE C.7 


(First meximum values 


90°00! 
3000 fps 
+0.0001 
0.0031 


+0,0004 
0.0045 


+0.0054 
0.0175 


-0.0108 
0.0198 


-0.0063 
0.0229 


(1) 


(2)) 


54° D0! 
2400 fps 
+0.0001 
0.0022 


+0.0002 
0.0036 


+0,0052 
0.0163 


-0.0105 
0.0189 


-0.0062 
0.0220 


we 


18° 4k 


1000 fps 


(3) 


(4) 


+0.0061 


OF OBES 7 


=O O0i2 


0.0174 


-0.005 8 


0.0210 


Based on Reinforced Concrete Arch, qd. = B/4, Sl 100 psi 


h, = horizontal movement of windward haunch, in feet 
Pe = horizontal movement of leeward haunch, in feet 
ae = vertical movement of crown, in feet. 


No initial upward movement of crown. 


No initial inward movement of windward haunch. 


Oe 
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APPENDIX D 


GRAPHICAL PRESENTATION OF NUMERICAL RESULTS 


SS SxsS]E]RE SSS EE Ee 


The results of the numerical work done in connection with this 
thesis are presented in this appendix in graphical form. Displacement-time 
"response curves are included for each of the problems solved. Three 


curves are plotted for each problem as follows: 


hy = horizontal displacement of the windward haunch 
oe = horizontal displacement of the leeward haunch 
Hs = vertical displacement of the crown 


The sign convention used in this appendix is the same as that stated in 
Section 1.2, i.e. displacements to the right and upward are considered 
positive. In addition to the response curves, the configuration of the 
dynamic model is illustrated for two of the problems at several significant 
times in their response histories. 

The problems are described in Chapter VIII and summarized in 
Figure 8.5. A set of sample calculations for one of the problems is pre- 
sented in Section 8.4. The significance of the results is discussed in 


Chapter IX. 
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